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Preface 


Beginning with Volume XX, the Deep Space Network Progress Report changed from 
the Technical Report 32- series to the Progress Report 42- series. The volume number 
continues the sequence of the preceding issues. Thus, Progress Report 42-20 is the 
twentieth volume of the Deep Space Network series, and is an uninterrupted foUow-on to 
Technical Report 32-1526, Volume XIX. 

This report presents DSN progress in flight project support, tracking and data 
acquisition (TDA) research and technology, network engineering, hardware and software 
implementation, and operations. Each issue presents material in some, but not all, of the 
following categories in the order indicated. 

Description of the DSN 

Mission Support 

Ongoing Planet ary /Interplanet ary Flight Projects 
Advanced Flight Projects 

Radio Astronomy 

Special Projects 

Supporting Research and Technology 
Tracking and Ground-Based Navigation 
Communications— Spacecraft /Ground 
Station Control and Operations Technology 
Network Control and Data Processing 

Network and Facility Engineering and Implementation 
Network 

Network Operations Control Center 
Ground Communications 
Deep Space Stations 
Quality Assurance 

Operations 

Network Operations 
Network Operations Control Center 
Ground Communications 
Deep Space Stations 

Program Planning 
TDA Planning 

In each issue, the part entitled “Description of the DSN” describes the functions and 
facilities of the DSN and may report the current configuration of one of the seven DSN 
systems (Tracking, Telemetry, Command, Monitor 8c Control, Test & Training, Radio 
Science, and Very Long Baseline Interferometry). 

The work described in this report series is either performed or managed by the 
Tracking and Data Acquisition organization of JPL for NASA. 
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Network Functions and Facilities 

N. A. Renzetti 

Office of Tracking and Data Acquisition 


The objectives, functions, and organization of the Deep Space Network are 
summarized; deep space station, ground communication, and network operations control 
capabilities are described. 


The Deep Space Network was established by the National 
Aeronautics and Space Administration (NASA) Office of 
Space Tracking and Data Systems and is under the system 
management and technical direction of the Jet Propulsion 
Laboratory (JPL). The network is designed for two-way 
communications with unmanned spacecraft traveling approxi- 
mately 16,000 km (10,000 miles) from Earth to the farthest 
planets and to the edge of our solar system. It has provided 
tracking and data acquisition support for the following NASA 
deep space exploration projects: Ranger, Surveyor, Mariner 
Venus 1962, Mariner Mars 1964, Mariner Venus 1967, Mariner 
Mars 1969, Mariner Mars 1971, and Mariner Venus'Mercury 
1973, for which JPL has been responsible for the project 
management, the development of the spacecraft, and the 
conduct of mission operations; Lunar Orbiter, for which the 
Langley Research Center carried out the project management, 
spacecraft development, and conduct of mission operations; 
Pioneer, for which Ames Research Center carried out the 
project management, spacecraft development, and conduct of 
mission operations; and Apollo, for which the Lyndon B. 
Johnson Space Center was the project center and the Deep 
Space Network supplemented the Manned Space Flight Net- 
work, which was managed by the Goddard Space Flight 
Center. The network is currently providing tracking and data 
acquisition support for Helios, a joint U.S./West German 
project; Viking, for which Langley Research Center provides 
the project management, the Lander spacecraft, and conducts 


mission operations, and for which JPL provides the Orbiter 
spacecraft; Voyager, for which JPL provides project manage- 
ment, spacecraft development, and conduct of mission 
operations; and Pioneer Venus, for which tlie Ames Research 
Center provides project management, spacecraft development, 
and conduct of mission operations. The network is adding new 
capability to meet the requirements of the Jupiter Orbiter 
Probe Mission, for which JPL provides the project manage- 
ment, spacecraft development and conduct of mission 
operations. 

The Deep Space Network (DSN) is one of two NASA 
networks. The other, the Spacefli^t Tracking and Data 
Network (STDN), is under the system management and 
techmeal direction of the Goddard Spaqe Flight Center 
(GSFC). Its function is to support manned^ and unmanned 
Earth-orbiting satellites. The Deep Space. Network supports 
lunar, planetary, and interplanetary flight projects. 

From its inception, NASA has had the objective of 
conducting scientific investigations throughout the solar sys- 
tem. It was recognized that in order to meet this objective, 
significant supporting research and advanced technology devel- 
opment must be conducted m order to provide deep space 
telecommunications for science data return in a cost effective 
manner. Therefore, the Network is continually evolved to keep 
pace with the state of the art of telecommunications and data 
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handling. It was also recognized early that close coordination 
would be needed between the requirements of the flight 
projects for data return and the capabilities needed in the 
Network. This close collaboration was effected by the appoint- 
ment of a Tracking and Data Systems Manager as part of the 
flight project team from the initiation of the project to the 
end of the mission. By this process, requirements were 
identified early enough to provide funding and implementa- 
tion in time for use by tlie fiiglit project in its fliglit phase. 

As of July 1972, NASA undertook a change in the interface 
between the Network and the flight projects. Prior to that 
time, since 1 January 1964, in addition to consisting of the 
Deep Space Stations and the Ground Communications 
Facility, the Network had also included the mission control 
and computing facilities and provided the equipment in the 
mission support areas for the conduct of mission operations. 
The latter facilities were. housed in a building at JPL known as 
the Space Flight Operations Facility (SFOF). The interface 
change was to accommodate a hardware interface between the 
support of the network operations control functions and those 
of the mission control and computing functions. This resulted 
in the flight projects assuming the cognizance of the large 
general-purpose digital computers which were used for both 
network processing and mission data processing. They also 
assumed cognizance of all of the equipment in the flight 
operations facility for display and communications necessary 
for the conduct of mission operations. The Network then 
undertook the development of hardware and computer soft- 
ware necessary to do its network operations control and 
monitor functions in separate computers. A characteristic of 
the new interface is that the Network provides direct data flow 
to and from the stations; namely, metric data, science and 
engineering telemetry, and such network monitor data as are 
useful to the flight project. This is done via appropriate ground 
communication equipment to mission operations centers, 
wherever they may be. 

The principal deliverables to the users of the Network are 
carried out by data system configurations as follows: 

(1) The DSN Tracking System generates radio metric data, 
i.e., angles, one- and two-way dopplef and range, and 
transmits raw data to Mission Control. 

(2) The DSN Telemetry System receives, decodes, records, 
and retransmits engineering and scientific data gener- 
ated in the spacecraft to Mission Control. 

(3) The DSN Command System accepts spacecraft com- 
mands from Mission Control and transmits the com- 
mands via the Ground Communication Facility to a 
Deep Space Station. The commands are then radiated 
to the spacecraft in order to initiate spacecraft func- 
tions in flight. 


(4) The DSN Radio Science System generates radio science 
data, i.e., the frequency and amplitude of spacecraft 
transmitted signals affected by passage through media 
such as the solar corona, planetary atmospheres, and 
plantetary rings, and transmits this data to Mission 
Control. 

The data system configurations supporting testing, training, 
and network operations control functions are as follows: 

(1) The DSN Monitor and Control System instruments, 
transmits, records, and displays those parameters of the 
DSN necessary to verify configuration and validate the 
Network. It provides the tools necessary for Network 
Operations personnel to control and monitor the Net- 
work and interface with fligiit project mission control 
personnel. 

(2) The DSN Test and Training System generates and 
controls simulated data to support development, test, 
training and fault isolation within the DSN. It partici- 
pates in mission simulation with flight projects. 

The capabilities needed to carry out the above functions 
have evolved in the following technical areas: 

(1) The Deep Space Stations, which are distributed around 
Earth and which, prior to 1964, formed part of the 
Deep Space Instrumentation Facility. The technology 
involved in equipping these stations is strongly related 
to the state of the art of telecommunications and 
flight-ground design considerations, and is almost com- 
pletely multimission in character. 

(2) The Ground Communications Facility provides the 
capability required for the transmission, reception, and 
monitoring of Barth-based, point-to-point communica- 
tions between the stations and the Network Operations 
Control Center at JPL, Pasadena, and to the JPL Mis- 
sion Operations Centers. Four communications dis- 
ciplines are provided: teletype, voice, high-speed, and 
wideband. The Ground Communications Facility uses 
the capabilities provided by common carriers through- 
out the world, engineered into an integrated system by 
Goddard Space Flight Center, and controlled from the 
communications Center located in the Space Flight 
Operations Facility (Building 230) at JPL. 

The Network Operations Control Center is the functional 
entity for centralized operational control of the Network and 
interfaces with the users. It has two separable functional 
elements; namely. Network Operations Control and Network 
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Data Processing. The functions of the Network Operations 
Control are: 

(1) Control and coordination of Network support to meet 
commitments to Network users. 

(2) Utilization of the Network data processing computing 
capability to generate all standards and limits required 
for Network operations. 

(3) Utilization of Network data processing computing 
capability to analyze and validate the performance of 
all Network systems. 

The personnel who carry out the above functions are located 
in the Space Flight Operations Facility, where mission opera- 
tions functions are carried out by certain flight projects. Net- 
work personnel are directed by an Operations Control Chief, 
The functions of the Network Data Processing are: 


(1) Processing of data used by Network Operations Control 
for control and analysis of the Network. 

(2) Display in the Network Operations Control Area of 
data processed in the Network Data Processing Area. 

(3) Interface with communications circuits for input to 
and output from the Network Data Processing Area. 

(4) Data- logging and production of the intermediate data 
records. 


The personnel who carry out these functions are located 
approximately 200' meters from the Space Flight Operations 
Facility. The equipment consists of minicomputers for real- 
time data system monitoring, two XDS Sigma 5s, display, 
magnetic tape recorders, and appropriate interface equipment 
with the ground data communications. 
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DSN VLBI Articles 


This issue of the DSN Progress Report contains two articles which relate to current 
DSN VLBI activities. 

The first article describes the functional requirements and system design of the DSN 
VLBI System, Mark 1-79, a new system being implemented in the Network. 

The second article provides a tutorial introduction to the technology of VLBI. The 
basic concepts are described, and equations are derived by which various parameters 
related to interferometry can be approximated. 

At the back of this issue, the interested reader will find a brief VLBI history as well as 
an extensive VLBI bibliography. 
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The DSN VLB! System, Mark 1-79 

B. D, LMuihall 

TDA Engineering Office 


The DSN VLB! System has been established as a new Network System. This article 
describes the system and discusses the system functional requirements. 


I. Introduction 

Very long baseline interferometry is a new capability being 
implemented m the Network to support flight project naviga- 
tion requirements and provide station frequency standard 
stability monitoring of sufficient accuracy to effectively vali- 
date hydrogen maser performance. 

In order to manage this major new implementation, the 
DSN VLBI System was established, and a review of the system 
functional requirements was held om|24 February 1978. A 
description of the system and the phases of its implementation 
follow. 

II. Definition 

Very long baseline interferometry is a method of measuring 
the time of arrival of a radio signal at two locations very 
distant from each other on the Earth’s surface. From the 
measurement of difference m time of arrival, the position of 
the radio source and/or several other parameters of the prob- 
lem can be determined. These other parameters mclude 
Universal Time One (UTl) (the instantaneous rotational angle 
of the Earth), polar motion, the relative position of the two 
stations, and the time offset and rate of change of the clocks 
at the two stations. 


III. System Characteristics 

The DSN VLBI System can be characterized as a new 
implementation in the 64-meter subnet and the Network Oper- 
ations and Control Center. The implementation modifies exist- 
ing station equipment; implementation of the 34-meter subnet 
is planned later (e.g., the multimission open-loop receivers, tlie 
occult ation data assembly, and the digital recorder assembly). 
New hardware and software will be implemented in the Net- 
work Data Processing Area at JPL Pasadena. 

VLBI requires at least two stations and a central facility for 
correlation and estimation processors. 

The station diagram consists of an antenna pointed at the 
signal being used to make the measurement, low-noise micro- 
wave amplifiers, and receivers for receiving the signal and 
converting it to zero offset. At this point, analog- to-digital 
converters and multiplexers convert the data to digital form 
and direct them to the appropriate device for recording. The 
diagrams show various implementation phases. These are 
Block 1 (Phase 1, 2, and 3) and Block 2. 

Block 1 handles data differently than does Block 2. As 
shown in Fig. 1, the Block 1 system records the data on 
computer-compatible tape for playback and transmission to 
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the Ground Communications Facility (GCF) wideband data 
lines. 

The DSS configuration for Block 1, Phase 1 is shown in 
Fig. 3. GCF transmission capability will be achieved by July 
1979 by employing the Block IV receiver and DSN Advanced 
System IF to video converters to provide data to the DSS 
Radio Science Subsystem. By July 1980, the multimission 
open-loop receiver will be transferred. This configuration 
(Fig. 4) is called Block 1, Phase 2. The Block 1 configuration is 
completed with Phase 3 by the addition of the Block 1 proces- 
sor (Fig. 6). 

The Block 2 interferometer (Fig. 2) requires such large 
volumes of data and is used so relatively infrequently that the 
data are recorded on a digital instrumentation recorder, the 
digital recorder assembly. These tapes are then shipped or 
transported to the correlation facility. Figure 3 shows the DSS 
configuration with the Block 3 Monitor and Control Subsys- 
tem, which centralizes the collection of configuration and 
calibration data. 

The correlation facilities are also divided by Blocks 1 and 2 
in Figs. 1 and 2. These two processor facilities employ the 
same design and equipment, special-purpose hardware for cor- 
relation, control computer, and correlation software. However, 
as before, the method of providing data to the correlator is 
significantly different for Block 1 than for Block 2. In the case 
of Block I, the computer edits the data received over wide- 
band data lines and provides two serial streams of bits repre- 
senting the information recorded at the two separate stations. 
Tliese data are accepted by the correlator and the cross- 
correlation is performed. This facility serves only DSN opera- 
tional functions, and consequently, will be located in the 
Network Data Processing Area (NDPA). 

In the case of Block 2, the same correlator design is em- 
ployed; however, the data are provided by playback from the 
digital recorder assembly. This facility will be located outside 
of the Network Operations Control Center (NO^C) at another 
location within JPL, Pasadena. 

IV. DSS Configuration 

The DSS configuration (Figs. 3, 4, and 5) includes the 
portions of the Ground Communications Facility (GCF) that 
exist at the station. The configuration at both VLBI stations is 
identical. 

The signal is received and passes through the Antenna 
Microwave Subsystem. In this subsystem, the phase calibrator 
tones generated by the Receiver-Exciter Subsystem are 
injected as high up, as near the front end, as possible. These 


tones provide a calibration of any drift between the injection 
point and the analog-to-digital converters. The Antenna Micro- 
wave Subsystem also provides water vapor radiometer data to 
the DSS Monitor and Control Subsystem for logging. 

The signal next passes to the Receiver-Exciter Subsystem, 
where the signal is converted from RF to 300 MHz IF. The 
signal is then divided into eight IF-to-video converters. These 
IF converters are in turn interfaced to the A-to-D converters in 
the DSS Radio Science Subsystem. At this point the data 
either go through a multiplexer to the occultation data assem- 
bly for recording for Block 1 functions or the digital recoiding 
assembly for Block 2 functions. After the observation, if 
Block 1 configuration is in use, the data are played back from 
the real-time record, blocked up in GCF blocks, and provided 
to the station wideband data assembly for transmission. If 
Block 2 functions are being used, the recording on the digital 
recording assembly is shipped back to the correlation facility. 

The Monitor and Control Subsystem has the functions of 
configuring the stations and collecting various calibration and 
configuration data. These calibration and configuration data 
are provided to the digital recording assembly for logging and 
transmission over wideband lines with the actual VLBI data or 
they are provided over high-speed lines for monitoring pur- 
poses. If Block 2 data are being recorded, then the computer 
recorder on the occultation data assembly is used for collect- 
ing the calibration and configuration data and providing them 
with the digital recording assembly tape. The Monitor and 
Control also receives predicts from the Network Operations 
and Control Center, which are disseminated to the appropriate 
subsystems, including Antenna Mechanical, for pointing the 
antenna. 

V. Ground Communications Facility and 
Network Operations and Control 
Center Configurations 

Figures 6 and 7 show the entire GCF, including the station 
portion and the NOCC configuration. The NOCC configura- 
tion is further divided between the equipment in the network 
data processing terminal and the network data processing area. 
The network data log provided by the GCF can be used for 
collecting the VLBI and the calibration and configuration data 
for Block 1. The VLBI real-time monitor collects configura- 
tion and calibration data to verify the performance of the 
network in real time. 

The VLBI data are replayed and metered out to the VLBI 
correlator assembly that is under the control of the VLBI 
processor assembly, a computer with the correlation program 
running in it. The two data streams are cross-correlated and 
the resulting fringes are recorded by the VLB! processor 
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assembly. After the correlation process, the post-correlation 
tape is replayed through the same computer and post-correla- 
tion and estimation processing are performed. The final prod- 
uct is clock sync, UTl, the polar motion. 

The Block 2 correlator does not employ the GCF transmis- 
sion, and consequently, the tapes of the digital recording 
assembly and occultation data assembly are brought to this 
facility, where the con elation process is quite similar using the 
same design correlator and identical software except for addi- 
tional features and functions required of Block 2 over Block 1 
capability. 

The Display Subsystem provides station configuration data 
to the network operations and, from this area, controls are 
sent to the station. In the Support Subsystem, the VLBI 
predicts program is used to generate the sequence of events for 
the station, including the antenna pointing angles, which are 
then transmitted through GCF for Blocks 1 and 2 data acquisi- 
tion. 


VI. Block 2 VLBI Future Requirements 

The future of the Block 2 VLBI processor beyond the DSN 
requirements for radio source catalog maintenance are being 
studied and defined. The evolutionary growth of this facility, 
so that* it can serve the radio astronomy community and the 
National Geodetic Survey requirements as well as the DSN 
needs, is briefly described. 

Table 1 shows the evolutionary steps for the processor for 
radio astronomy applications. The columns list the number of 
stations or radio telescopes which can be simultaneously cor- 
related, the number of BWS channels that can be correlated in 
parallel, and finally, the need dates for the expanded capabil- 
ity. The first step would be the two-station, 8-channel corre- 
lator capability. This capability employs digital instrumenta- 
tion recorders with 24 tracks, which are capable, of extension 
to 28 tracks. In order to reduce tape usage, the eiglit BWS 
channels are blocked up after A and D conversion and written 
across the 24 tracks. This requires deskewing upon playback. 

The first step in increasing capability would be to add 
Mark II compatibility for system validation and for processing 
of Mark II tapes, which will still be in use when the processor 
becomes operational. The need for this Mark II capability 
could be as early as 1979, even though the processor is not 
operational at that time. At the same time, capability should 
be provided to read eight tracks of data directly off the tape 
and provide each track to the respective correlator channel so 
that the processor would be capable of performing Mark III 
processing with multiple passes through the tape. The capabil- 


ity for parallel reading without deskewing could also be used 
as early as 1979. 

The second step would be either to expand to three stations 
and three baselines with eight channels or expand the two 
telescope processors to 28 channels, all of which would be 
read directly off independent tracks of the 28-track recorder. 
The expansion from eight to 28 tracks increases processor 
thruput by a factor of 3.5 if all 28 tracks are parallel BWS 
channel recordings. On the other hand, going from two to 
three stations increases thruput by a factor of 3 (three base- 
lines versus one). This latter step is considerably more expen- 
sive, since an additional tape drive is required. 

By 1981, a third step could be justified in terms of the 
needs of the radio astronomy community. This third step 
would be to provide three baseline correlators with 28 BWS 
channels per station. At this time, a decision point is reached. 
For the ultimate capability, the number of radio telescopes in 
the radio astronomy network needs to be identified. Presently, 
there are four, perhaps five, radio telescopes in the United 
States which will shortly be equipped with Mark III instrumen- 
tation. By 1981, we should be able to say if the number of 
radio telescopes in the Network Users Group will increase to 
eight or possibly ten. 

Based on the ultimate number of potential telescopes, the 
decision can be made to proceed to increase the processor to 
four, possibly five, stations, leaving room for ultimate expan- 
sion to the final number. At this point also, the existing 
hardware correlator assembly probably needs to be repackaged 
to make room available for cabling to interconnect all the 
other station correlators so that in Step 4, when expansion is 
made to four or five telescopes, the repackaging would leave 
enough room for additional cabling and modules to allow all 
the necessary cross-correlations to be performed in parallel. 

The expansion in BWS channels beyond 28 may also be 
required, but this decision can probably be made later. The 
increase in hardware for expansion in channels goes by n, 
while the increase for expansion in telescope goes by and 
consequently the ultimate number of stations to be processed 
simultaneously is a more critical decision than the ultimate 
number of BWS channels. 

Table 2 shows the growth of processor capability required 
primarily for geodetic but also for geodynamic (UTl -polar 
motion) applications. Some differences between these require- 
ments and the astronomy requirements are: 

(1) BWS channel widths for radio astronomy range from 
0.125 MHz to 2 MHz (possibly both narrower and 
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wider widths will be required), while geodesy requires 
2 MHz to 25 MHz channel bandwidth. 

Radio astronomy requires independent track recording, 
while geodesy requires maximum tape efficiency while 
keeping the ability to speed up or slow down the data 
rate during the observation, depending on the source 


strength. This requires blocking the data and writing 
across all tracks regardless of the tape speed. 

The immediate impact of these Block 2 processor plans is 
the design requirement to be able to expand the correlator to 
multiple baseline. The Block 1 processor is not affected since 
it will remain a two-station device. 



Table 1. Radio astronomy applications ' 


Step 

Number of 
telescopes 

Number of 
BWS channels 

Date needed 

0 

2 

8 

(blocked and wntten 
over 24 tracks) 


la 


(Mark II compat) 

1979 

lb 

2 

8 

(parallel, each track 
recorded independently) 

1979 

H 

3 

8 



2 

28 


3 

3 

28 

1981 


(Decision points for ultimate number of telescopes) 

4 

4 (5?) 

28 

1983 


Table 2. Geodetic and/or geodynamrc applications 


Step 

Number of 
antennas 

Number of 
BWS channels 

BWS 

channel width 

Date needed 

1 

2 

2 

2 -25 MHz 

1981 

2 

2 

6 

2- 8 MHz 

1981 

3 

6 

6 

2-, 8 MHz 

1982/3 
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Fig. 4. DSS-GCF block diagram, Block 1 and 2 VLBI configuration before Block 3 BMC 
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Fig. 7. GCF-NOCC block diagram, Block 2 processor 
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This article gives a tutorial presentation of the basic principles underlying Very Long 
Baseline Interferometry (VLBI) using bandwidth synthesis. Although many subtle details 
are iptored, the article presents the basic signal processing approach and summarizes 
results showing the tradeoff of measurement accuracy with spanned bandwidth, source 
strength, antenna size and efficiency, system noise temperature, and data volume. Results 
pertaining to minimization of required antenna time for a given baseline measurement 
accuracy are also discussed. 


I. Introduction 

A concise description of VLBI, given in Ref. 1, is repeated 
below. 

In very long baseline interferometer (VLBI) measure- 
ments, the radio signal produced by a distant source is 
recorded simultaneously at two radio antennas. Because 
of a difference in raypaths, the signal will be delayed in 
time at one antenna relative to the other. By cross 
correlating the- two signals, the time delay and/or its 
time derivative may be determined. In addition, cor- 
related amplitude measurements can yield source 
strength and structure. If the radio signal is generated by 
an extragalactic object, the radio source may be regarded 
as a fixed object because of its great distance. In this 
case, the time dependence of the time delay is generated 
by the Earth's motion but depends, of course, on the 
source location and the baseline vector between the two * 
antennas. In general, measurement of the time delay 


and/or its derivative for many sources can lead to a 
least-squares determination of source locations, the base- 
line vector, and Eartli motion parameters, such as UTl 
and polar motion. 

Although much has been written about VLBI, a need was 
felt for a tutorial introduction to the basic principles involved. 
It is much easier to understand the more detailed analyses of 
VLBI after the central ideas are understood. Thus this tutorial 
concentrates only on the basic principles and ignores many 
subtle details that are discussed thoroughly in the literature. 

II. Monochromatic Source (Ref. 2) 

Consider two antennas separated by a distance (baseline) B 
receiving a signal from a distant extraterrestrial source. The 
geometry is shown in Fig. 1. If the source is assumed to be 
monochromatic and noise (receiver and background) is 
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neglected, the received signals at the two antennas can be 
represented as 

V^(t)=A^ cos lirf^t 
and 

^2(0 = ^2 ~ '^p 

where 

A^,A^ represent the strength of the received signals at 
antennas 1 and 2, respectively 

= (B cos \jj)/c is the time delay (typically a slowly 
varying function of time due to the earth’s 
rotation) in the reception of a given wavefront at 
antenna 2 relative to antenna 1 

c is the velocity of liglit 

These signals are processed by the system shown schemati- 
cally in Fig. 2 to produce a cross-correlation function 
?)• Prom Fig. 2, 

Ml COS27T//] [Mj COS(27 t4^j?- 0j)] 

= —^ cos [27r(4-4^j> + 0j 
^1^1 

+ cos [27 t(/; +4oP^- (2) 

where 

fj^oi effective local oscillator mixing frequency at 

station 1 

IS the amplitude of the local oscillator at station 1 

represents phase shifts due to the receiver and local 
oscillator (Phase shifts due to other propagation ef- 
fects such as the ionosphere can also be included in 
this term.) 

Tlie low-pass filter is designed to pass the difference term and 
reject the sum. term. (Assume for simplicity that the filter gain 
is 1 for the difference frequency. Equivalently it can be 
lumped in with and ^j.) Thus the recorded signal at 
station 1 is given by 

cos[27T(r,-4^^)f + 0j (3) 


In an entirely similar fashion, 

[ 257(4 ■ ^ + ^2] ( 4 ) 

As indicated in Fig. 2, Vu^Ct) and V 2 i(t) are typically 
recorded on magnetic tape (digitally in most cases, as discussed 
later) and shipped to a central facility for correlation and 
further processing. 

Correlation consists of advancing the recorded signal V 2 i(f) 
from station 2 by a model time that approximates the 
expected geometric delay r^, taking the product of the signals 
and again low-pass filtering. This is equivalent to multiplying 
and + tJ and low-pass filtering. From Eq. (4), 

A.M^ 

^2^2 

= COS [277(4-4.02)? 

+ 2^CT,-4:o2)"^-27r4r^ + 023 (5) 

Taking the product of and expanding 

into sum and difference terms, and noting that the low-pass 
filter removes the sum (high-frequency) term yields 

where 

A A M. M 
8 - 

and 

+ 27747^ + dj- 02 

The cross-correlation function t) is processed 

further (as will be discussed in a later section) and the delay 
as well as the delay rate which are the basic VLBI data 
types, are determined by analyzing the phase of the cross- 
correlation function. Note that 0^(f) may also be written in 
the form 

= 27r(4o2 - 27 t4^^ 
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where 


( 12 ) 


Ar = r ^ T 

g m 

A. Angular Sensitivity 

The geometric delay enters into the phase of the cross- 
correlation function as the ‘term Suppose it is 

possible to detect a variation in the phase of 27 t radians (a 

one-fVinge change). This corresponds to a change in geometric 
delay (assuming the other terms of do not change) given 
by 


27rj^Ar^ = 27T 


( 8 ) 


-0.0027 

sin 


arc sec 


or 


As stated in Ref. 2, “This illustrates the potential of VLBI: a 
change in position of arc sec causes detectable changes in 
the cross correlation. At the present though, there is an ambi- 
guity in determining which fringe the source is on.” That 
ambiguity problem is brought out in the next two sections. 

B. Phase Measurement 

The cross-correlation function R(r^, , t) is multiplied by 

the cosine and sme of a model phase 0^(0 and integrated over 
an interval of T seconds (the integrator eliminates the sum 
frequency term of the product) to produce the so-called cosine 
and sine-stopped fringes 




Now the relationship between small changes in and small 
changesin the^angle to the source 0 can be- approximated by 


9^. 

At ^ A\jj 


Substituting = (B cos \p)/c yields 


, -B sm 0 , , 

Ar = A0 

Sc 


( 9 ) 


or 


A0 = 


-c 


. At 
5 sin 0 s 


( 10 ) 


and 






(13) 


The phase of the cross-correlation function can, in principle, 
be determined (within a multiple of 27 t) as 

♦ *)-<> (0 = «n- 


or 


M (^» T ) 

0X0-0 (0 + tan-^ % ", 

ixXt.rJ 


(14) 


Using ATg = V gives 


In actual practice the phase is determined by fitting sinusoids 
to the cosine and sine-stopped fringes. 


A0^ 


-c 


sin 0 7^ B sin 0 


(11) 


since c/f^ - X, where X is the wavelength of radiation at fre- 
quency/^. 

As an example, suppose B = 10,000 km and X= 13 cm 
(corresponding to an S-band signal). Then 


Ar- -1.3X10“^ ‘1.3X10"® , 

A0 = p rad 


1 X 10'^ sin ^ 


III. Polychromatic Source and 
Bandwidth Synthesis 

As noted above, the phase of the cross-correlation function 
t) due to geometric delay can be determined only 
within a multiple of 27 t (i.e., 0^0 = + Im). If the 

source consists of several spectrd lines (polychromatic) it is 
possible to determine the delay- *by combining measurements 
from the cross-correlation functions of the individual spectral 
lines. As an. example, suppose the source consists of two 
spectral lines at frequencies and / 2 * phases of the 
cross-correlation functions due to the geometric delay term are 
then measured as 
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^>j.j(0 = 27Tj^i'r^ + 27rm 


and 


ambiguities depending on the spanned bandwidth, accuracy of 
the a priori information, and uncertainty in individual phase 
measurements* 


+ 2nn 


IV. Wideband Source and Noise 


where m and n are integers (see Fig. 3)* If the a priori 
uncertainty in slope iS" of phase versus frequency is less than 
A5j 2 corresponding phase points may be 

connected (i.e., the relative phase ambiguity resolved) and the 
delay is determined by the slope of the connecting line, 
since 

or 


It is important to note that the uncertainty in S includes 
uncertainties in both and instrumental delays. The instru- 
ment delay uncertainty can be essentially eliminated by using 
a phase calibrator. In addition, There is of course some uncer- 
tainty in the measurement of the phase <^y(r) due to noise as 
indicated by the error bars in Fig. 3 . (Noise is considered in a 
later section.) If the uncertainty in the measurement of 
at a given frequency is then the uncertainty in (neglect- 
ing certain errors due to ionospheric effects) is given by 


A. Cross-Correlation Function 

In most cases the source used in VLBI (a natural radio 
source for example) is broadband and can be approximated as 
white noise over the band of interest. In this case the correla- 
tion function can be shown (Ref. 1) to be of the form 

>0 = cos 4>^e) (18) 

where 

Ar = T - T 

g m 

+ 2nfLt + 27 ^ 4 ^ 2 '^^, ^1 ■ 

Note that the correlation function differs from the mono- 
chromatic case in that the amplitude of the correlation func- 
tion Z)(Ar) is a function of rather than a 

constant. For example, if the receiver has a rectangular band- 
pass of width W centered at/^, the amplitude of the coi rela- 
tion function is given by (Ref. 1) 


V2a^ 


= 


2 .( 4 , - 4 ,) 


(17) 


where the >/2 results from incoherent differencing of phase 
errors. The point is that for a given the larger the term 
fs 2 “ fs\ (commonly called the spanned bandwidth), the* 
smaller will be the delay uncertainty This technique is 
known as bandwidth synthesis. 

In this simple example, increasing the spanned bandwidth 
would, also require more accurate a priorUnformation in order 
to resolve the relative phase ambiguity. This dilemma can be 
avoided if an additional intermediate channel is added as 
shown in Fig. 3. Assume the relative phase ambiguity is re- 
solved between channels 1 and 3 using a priori information . 
The more precise estimate of determined by using chan- 
nels 1 and 3 may now be used to resolve the ambiguity 
between the outside channels 1 and 2. Obviously several inter- 
mediate channels may be required to resolve the relative phase 


where is a constant. This dependence on model delay 
reflects the accuracy with which the two signals have been 
aligned in time. 

Note further that, in the case of a broadband source, delay 
information can be obtained by analyzing the amplitude as 
well as the phase of the cross-correlation function. To deter- 
mine from the amplitude of the cross-correlation function, 
the model delay is adjusted to maximize £)(Ar), at which 
point At = = 0 or = r,„. 

The accuracy with which delay can be determined using the 
amplitude of the cross-correlation function depends on the 
width of the peak of Dipa) versus At. In the case of a 
rectangular bandpass, the first zero of D(At) (see Fig, 4) 
occurs at 7rK^Ar= tt or At= 1/H^. Thus the accuracy with 
which delay can be determined using this technique is propor- 
tional to the reciprocal of the channel bandwidth. Typically 
this bandpass shape can be used to estimate the delay t^ with 
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an accuracy of 0 . 01 / 2 IV to 0 . 1 / 2 IV seconds ( 2.5 to 25 nsec for 
^ PV= 2 MHz). 

B. Bandwidth Synthesis 

Rather than use the amplitude of the cross-correlation 
function to determine the delay, which would require record- 
ing a wide-bandwidth signal, several narrower bandwidth 
channels are used (see Fig. 5). Cross-correlation functions are 
computed for each channel, and the delay is determined by 
analyzing the phases of the cross-correlation functions and 
using bandwidth synthesis as explained above for the poly- 
chromatic source. -However, the amplitude of the single- 
channel cross-correlation function can be used to improve a 
priori knowledge of prior to performing bandwidth 
synthesis. It can be shown that for a given delay uncertainty 
, bandwidth synthesis requires less data than the amplitude 
technique for many cases of practical interest (Ref. 3). 


J, B. Thomas has drawn up an alternate way of looking at 
the bandwidth synthesis process shown schematically in Fig. 6 . 
First of all, the amplitude of the correlation function is used 
to give an unambiguous delay measurement based on the 
channel bandwidth IV (2 MHz in this case). This delay measure- 
ment can then be used to resolve the delay ambiguity resulting 
from the most closely spaced channel pair (5 MHz in this 
example). The delay ambiguity results from the phase slopes 
spaced by 3 in Fig. 3 [Ar^ 3 = (1/27t) AS^ 3 ]. Once the 
ambiguity in the most closely spaced channel pair is resolved, 
the more accurate measured delay obtained is in turn used to 
resolve the ambiguity in the more widely spaced channel pair 
(20 MHz in this case). Note, that both the uncertainty in the 
delay measurement and the ambiguity spacing decrease as the 
channel separation increases (see Eqs. 16 and 17) as indicated 
in Fig. 6 . 


V. Effect of Noise and Digital Processing 

When noise is included, the cross-correlation function may 
be written 


where 


stopped fringe signal-to-noise ratio 
< 7 ^ - rms fringe noise 


^ai*^a 2 ~ radio source temperatures (K) at sta- 

tions 1 and 2 , respectively 


^sV^s 2 “ system noise temperature (K) at stations 1 
and 2 , respectively 

W= channel bandwidth (Hz) 

T = total integration time (sec) 


In the actual implementation of this technique each 
channel is hard-limited and sampled ( 1 -bit quantization) at a 
rate 2 IV. The effect of hard-limiting for low signal-to-noise 
ratios is to reduce the signal to noise ratio S/N by a factor of 
2 /tt (Ref. 4) and thus for the actual system 


1^2l 


( 22 ) 


Note that maximum SjN requires that the model delay 
must be accurate enough to keep Z)(At-) near its maximum 
Z)( 0 ) to provide the most accurate measurement in the 
presence of noise as described below. In addition, further 
processing requires “seeing” fringes, which in practice requires 
that at least 5 and preferably 10. 

A. Phase and Delay Measurement Accuracy 
(Refs. 5 and 6) 

The error in determining 0^(/) - the stopped 

fringes may be illustrated by noting that r^) and 
ju^(f, r^) are projections of a vector consisting of the signal 
vector plus a noise vector on two orthogonal axes as shown in 
Fig. 7, where 


'^m’ COS (0 + n{i) ( 20 ) 

where w(r) is a random function due to receiver and back- 
ground noise. A rather involved analysis (Refs. 4 and 5 ) shows 
that 


5 _ -PCO) 

^ 


Ittwt 

/ al a2 

T T 

V si s2 


( 21 ) 


and 




«,(0 


( 23 ) 


Since the noise is spherically symmetric, it follows (see Fig. 7) 
that for small rms phase deviations, 


/ noise \ _ noise _ 1 _ N 

\ signal / ^ signal S/N S 


(24) 
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A more rigorous analysis (Ref. 5) yields the same result, 
and thus (assuming the relative phase ambiguities have been 
resolved) the system noise-limited delay accuracy is given 
by ^ 


■\/2o 


(25) 


where 

- mean sampling rate in megabits/sec for each channel 
T = total integration time in seconds 

- spanned bandwidth in MHz 


Then 


where now refer to the center frequencies of 

channels 1 and 2, respectively. Defining the spanned band- 

width Rj = /^2 "4i yields 


a- (cm) = 




BJd,d, 

S 12 . 


T T 
e,e,ST 

l 2 r 


cm 


(29) 


Substituting (Ref. 7) 


p =ii^ 
2 K 


Je.^d^ 

I 4 j 


1 = 1,2 


(26) 


(27) 


where 


= spanned bandwidth in MHz 
/ = correlated flux of radio source in Janskys 
e. = antenna efficiency i- 1,2 
d, = antenna diameter in meters 
= 5.26 X 10^ 


where 

= antenna polarization loss 

J = correlated flux of radio source in Janskys (Jansky = 
10-26 ^/j^2 Hz) 

= antenna efficiency 

== antenna diameter in meters 


Again it is important to note that Eq. (29) gives the 
standard deviation in the measurement of the delay due to 
thermal noise only. Other error sources such as instrumental 
effects, propagation (ionosphere), modelling, etc., must also be 
considered. The real importance of the equation is to show the 
manner in which various parameters affect VLB! accuracy. 

As an example, using the following values for a typical 
Mark II observation, 


K = Boltzmann’s constant {K = 1.38 X 10"^ ^ Joule/ 
Kelvin) 

gives 


== 


1.24X 10^ 


(T j 
^ Sl ^ S 2 


sec 


(28) 


It is sometimes convenient to give in cm by multiplying 
by the speed of light (3 X lO^^cm/sec). If, m addition, the 
spanned bandwidth B^ is measured in MHz (10^ Hz) and if the 
sampling rate is twice the channel bandwidth so that the 
number of megabits of data A" is given by 


= 5.26 X 10"^ 

/ ^ 1 Jansky 

= 4 Mbps 
J? = 40 MHz 

S 

d^ = d^=6Am 

€j ~ 0.55 

T = 150 sec (2-1/2 min) 

yields 


2WT 

Mbits = ^r 
10^ '' 


. 5.26X10“ / (30)(30) 

% ~ (40)(1)(64)(64) y (0.55)(0.55)(4)(1S0) " 
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Note that if the spanned bandwidth and antenna param- 
eters are held fixed, varies as where 

N- total number of bits for the observation. In the above 
example, 

N = 4(150) = 600 Mbits =.6 X 10® bits 

The sampling rate S^. is proportional to the channel band- 
width (2 MHz for the Mark II VLBI system) and thus, if the 
channel bandwidth is increased, the integration tune T can be 
decreased while maintaining a given 0^^ as long as = jS'^T’is 
held constant. A wider channel bandwidth may allow a more 
accurate single-channel delay to be computed, however, which 
affects the resolution of phase ambiguities in the bandwidth 
synthesis process. 

It should also be pointed out that T refers to integration 
time, during which the same radio signal from the source is 
recorded' at two stations. If short segments of the signal are 
recorded, careful timing is required to insure signal overlap. 

vr. Measurement of Source Position and 
Earth Motion Parameters 

As stated in Ref. 1 for an extragalactic radio source, ‘'the 
time dependence of the time delay is generated by the earth’s 
motion but depends, of course, on the source location and the 
baseline vector between the two antennas.” Epr example, the 
relationship between the geometric time delay r^, baseline By 
and source position is given by (Ref. 2) 


~c 5 sin 5^ + cos 5 cos 5^ cos “ a)] (30) 

where 

B - length of the baseline 

d = declination of the source 

5^ = declination of the baseline 
o 

a = right ascension of the source 
= right ascension of the baseline 
c - velocity of light 

In addition, the instantaneous spin axis of the earth (which 
is the reference for the source and baseline coordinates in 
Eq. 30 above) changes both with respect to an inertial refep 
ence frame (precession and nutation) and the earth’s crust 
(polar motion). 


A number of delay measurements at various times and for 
various sources can be used to make a least-squares determina- 
tion of source coordinates, the baseline vector, and earth 
motion parameters. A discussion of some of the tradeoffs 
involved in this type of measurement is given in the next 
section using the baseline as an example. 


A. Baseline Accuracy 

The standard deviation of the error in baseline measure- 
ment is approximately given by (Refs. 8 and 9) ' 



where 

“ number of parameters solved for in the multi- 
parameter fit 

number of observations 

O DS 

A = constant ranging from 1 to 4 depending on a 
covariance analysis based on source locations, etc. 

and 

a^=o^ + 

X o 

where 

a = that given 5>reviously 
'^g 

G - standard deviation of other uncertainties such as 

o 

ionospheric effects, tropospheric effects, uncertain- 
ties in source positions (if taken as given), frequency 
standard deviations from linear performance in time, 
etc. 

For example, consider solving for three earth parameters 
and four clock parameters using 28 observations, with = 
0.7 cm (see previous calculation) and .4=4. Assume for sim- 
plicity that 0 ^ = 0; then 

a = (0.7)(4) 1 -4 cm 

The total number of bits required in this case is given by 

6X 10® = L68X 10*“ bits 
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B, Delay Rate Accuracy 

For a given standard deviation in phase the standard 
deviation in delay rate may be approximated by noting that, 
from Eq. (6), 




= -pp cm/sec 




^LOl^ 


(32) Dividing by the wavelength of the signal at frequency /^, 


and thus 

Gi =27r/a- 
or 


% 

Since for an integration interval of T seconds, can be 
approximated by taking the difference in phase over the inte- 
gration interval and dividing by T, 



radians/sec 


and thus 


a- 


- 2 ^ sec/sec 


(34) 


(35) 


Note that there is no ambiguity problem in determining 
delay rate and the effective spanned bandwidth for delay rate 
is the rf channel frequency Since varies 
varies as Note that although for a given spanned ban(£ 
width and antenna parameters, Oj depends only on the 
number of bits N= Sy.T (thus allowing a tradeoff between 
integration time and channel bandwidth), depends 
explicitly on the integration time T (varying as 7“^ for a fixed 
CF^ ). From Eq. (26), 


(7 


V2 \ 


Substituting into Eq. (35), 



Using the previous Mark II example yields 


40X10^ 

(3X 10’<*)(150) 


0.7 = 6.22 X 10"® Hz- 


Accurate measurement of strengthens the multiparameter 
fit for baseline measurement, especially the z-component 
(along the earth's rotation axis), as schematically indicated by 
Fig. 8. A tradeoff of r vs f accuracy for a given baseline 
component accuracy thus exists. This affects observation 
strategy and has not; been analyzed in de tailgas* yet, sLiice the 
experimental data on typical at X-band' are” only now, 
becoming available. ^ 


VII. Required Antenna Time for a Given 
Baseline Accuracy 

As discussed previously, the standard deviation of the error 
in baseline measurement is given by Eq. (31): 


a = a^A 



where 

ol = ol 

X o 

and 


'^s 


K, 

BJd^d^ 




(38) 


(7- = sec/sec 


Following an analysis by D. Rogstad, the total antenna time 
( 37 ) and number of bits Nrj. required for a given baseline 

accuracy are given by 
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and 


T =N . (T+t ) 

ant obs ^ 


(39) of = ^obs(^ ^m) subject to a constraint on a (or a^). 
The result (which can be derived using Lagrange multipliers) is 


where is the average time required to move from one 
source to another (= 4 min for a 64-m/64-m antenna pair), 
and T is the integration time per source (assumed uniform). 
For convenience, define 


Then 


/5ST 


(41) 



(42) 


and 


= (a; 


A^N 


+ a^) 

T 

S 


N 


obs 


= [a 



A^N 


N 


obs 


(43) 


Clearly, there is a tradeoff between T and in maintain- 
ing a constant o. Rogstad has determined the minimum value 


p 



This corresponds to 




Normalized curves of the results are given in Ref. 9, as well as 
a discussion of how should be handled. Rogstad has subse- 
quently obtained further results in this area. 


Vllf. Summary 

''TKis'tutorians intended to give the reader an appreciation 
of the basic principles of VLBL It should not be used as the 
sole basis of further analysis. A great deal of technical literature 
exists for those interested in further information and more 
detailed analyses. Especially recommended as a starting point 
are the articles by J. G. Williams (Ref. 2), A. E. E. Rogers 
(Ref. 10), and J. B. Thomas (Ref. 1, 4, and 1 1). 
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Fig. 3. VLSI relative phase ambiguity, resolution and bandwidth synthesis 



Fig. 4. Amplitude of cross-correlation function D(Ar) 

vs bit stream misalignment (At) pjg^ 5^ Frequency response of bandwidth synthesis VLSI system 
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Fig. 6. Schematic example of ambiguity resolution In bandwidth 
synthesis 



Fig. 7. Effect of noise on accuracy of phase Fig. 6. Standard deviation In Z-^omponent of baseline 
measurement as a function of and 
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Viking Extended Mission Support 

R. L Gillette 
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This report covers the period from 1 March through 30 April 1978, It reports on DSN 
support of Viking spacecraft activities during the period and continues reporting on the 
DSN Viking Command and J) aching support It also continues the reports on the status 
of Viking DSN Mark III Data Subsystem Implementation Project (MDS) related testing. 


I. Viking Operations 

A. Status 

The Viking Orbiter 1 spacecraft continued to make detailed 
photomaps of Mars during this reporting period and to return 
scientific data on Martian cloud patterns, temperatures, and 
water vapor in the atmosphere. 

To conserve gas, the Viking Orbiter 2 scientific observations 
of the poles of Mars were discontinued during this reporting 
period and the spacecraft commanded into a roll drift mode of 
operation. In this mode Orbiter 2 cannot observe Mars or sup- 
port relay links from Lander 2. However its radio signals were 
used to map gravity variations on Mars, caused by pockets of 
dense material called Mascons, which sliglitly change the or- 
bital speed of the spacecraft. 

The Viking Landers continued their investigation into tlie 
composition of the Martian soil. Plans were completed for 
automatic operation of both Landers with little or no instruc- 
tions from Earth starting in early June. 

B. Orbiter 2 Aft-Bioshield Jettison 

The Viking Orbiter 1 (VO-1) aft-bioshield (Fig. 1) was sep- 
arated shortly after Lander separation in 1976. The shock 
impact on the VO-1 radio subsystem was a momentary change 


in the uplink signal level, indicating a momentary out-of-lock 
even with an uplink signal present. Also, a transient occurred 
in the traveling wave tube helix current, which required a half 
hour to recover. Because of these effects, and a failure in a 
gyro power supply on Viking Orbiter 2 (VO-2) at Lander 
separation, the VO-2 aft-bioshield was not separated during 
the Viking primary mission. The science community was inter- 
ested in having this separation done during the extended mis- 
sion, so as to provide a wider field of view for instrument 
operation. 

The decision was made to separate the bioshield on March 3, 
1978. The 64-meter Deep Space Station (DSS 63) located in 
Spain was the primary tracking station supporting the bio- 
shield jettison with the 26-meter Spanish Deep Space Station 
(DSS 61) assigned as a backup if needed for commanding. 

The separation sequence went entirely without problems. 
The first indication of a successful separation was an abrupt 
change in the spacecraft two-way doppler from DSS 63. The 
VO-2 aft-bioshield separation is tlie last major spacecraft 
event planned on either Orbiter. 

C. Spacecraft Problems 

The VO-2 spacecraft developed a serious gas leak in the 
attitude control system in late March resulting in loss of half 
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of the spacecraft’s remaining gas supply. To conserve gas the 
Orbiter’s redundant attitude control system was switched on. 
To prevent depletion of the remaining gas due to leaks in the 
roll jets of the redundant attitude system, VO-2 was com- 
manded into a roll drift mode of operation while scientists 
prepared the most useful science sequences for the remaining 
lifetime of the spacecraft. 

D. Radio Science 

Radio Science activities and experimentation continued 
during March and April. These activities include near- 
simultaneous Lander/Orbiter ranging, VO-1 and VO-2 Earth 
occultation coverage, and the Gravity Field Experiment. 


II. Network Support 

Table 1 shows the Viking Extended Mission (VEM) track- 
ing support for March and April 1978. Noticeable during this 
period is the increase of 64-meter support (DSS 14, 43, 63) 
during April. Thisris* primarily due to the gas leak'prdblem on 
VO-2 requiring the additional 8-dB antenna gain from the 
64-meter network over the 26-meter network to make up for 
the 21-dB signal loss when VO-2 was commanded into a roll 
drift mode and the transmitter configured to the spacecraft 
low-gain antenna. 

Table 2 gives the total number of commands transmitted 
by the DSN for Viking Project during March and April 1978. 


III. DSN Mark III Data Subsystem (MDS) 
Implementation Testing and Status 

As indicated in the last report of this series, MDS test and 
training had been completed at all DSN stations except DSS 
11 at Goldstone, California. DSS 11 was released on 15 Jan- 
uary to begin their MDS upgrade and complete test and train- 
ing during this reporting period. 

A. DSS 11 Test Status 

Test and training activities for Viking started in March and 
consisted of three Operational Verification Tests (OVTs) and a 
DSN/MCCC System Integration Test (SIT). 

Only the third of the three Viking OVTs was considered 
successful. Problems in the Simulation Conversion Assembly 
(SCA) at DSS 11 and the Test and Training Computer in the 
Network Operations Control Area (NOCA) at JPL prevented 
testing of the telemetry system. Both the SCA and the Test 
and Training Computer are used by the DSN to simulate a 
spacecraft and generate simulated spacecraft telemetry data. 

The first System Integration Test scheduled for 29 March 
was cancelled due to SCA problems and a last minute priority 
placed on the Viking Mission Computers by the Viking Project 
to support real-time operations. A Viking SIT retest on 7 April 
was considered successful. 

After several successful demonstration passes with the Vik- 
ing spacecraft, DSS 11 was placed under DSN Viking Configu- 
ration Control, effective 26 April. 
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Table 1. V^M tracking support 1978 

DSS ^ 

Marcli 

April 

Tracks 

Hours 

Tracks 

Hops 

11 

19 ‘ 

167 

_ 

_ 

12 

5 

22 

1 ‘ 

8 

14 

37 

*272 

47 ‘ 

365 

42 

30 

201 

30 

192 

43 

15 

68 

27 

163 

44 

6 

43 

- 

- 

61 

25 

281 

15 

159 

62 

3 

22 . 

3 

‘29 

63 

30 

293 

42 

440 * 

Total 

170 

1369 

165 

1356 


Note: Number of tracks represent the summation of all Viking 
spacecrafts tracked. ^ Track time, ' in hours, represents 
actual scheduled station support. 


Table 2. Commands transmitted during March and April 1978 


DSS 

Commands 


March 

April 

11 

119 


12 

1 

374 

14 

1326 

3032 

42 

261 

1079 

43 

124 

275 

44 

3 

— 

61 

1073 

1503 

62 

461 

73 

63 

2597 

4005 

Total 

5965 

10341 
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Fig» 1. Location of aft bioshteld separation points 
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This article covers the time period from 1 January 1978 through 31 May 1978. 
Ongoing Pioneer Project operations (Pioneers 6 through 11) are briefly mentioned. The 
DSN prelaunch preparations and launch operations for the Pioneer Venus Orbiter are 
described in detail^ along with a brief update on the Pioneer Venus Multiprobe launch and 
encounter preparations. 


I. Pioneers 6, 7, 8, 9, 10, 11 

Pioneers 6- 9 have only been tracked on a few rare occa- 
sions, due to the higher priority of other DSN users, and all 
the spacecrafts’ status* appear to be unchanged. Pioneer 10 and 
11 are being tracked almost daily; Pioneer 10 mostly by 
DSS 63, and Pioneer 11 mostly by DSS 42 and DSS 62. Pie- 
cession maneuvers on Pioneer 10 and 11 are routinely being 
conducted approximately every four weeks. 

II. Pioneer 12 (Pioneer Venus Orbiter) 

A. Prelaunch Training and Testing for 
Launch and Cruise 

The DSN prelaunch test and training program was planned 
to follow the established sequence of: 

(1) Start DSN Operational Verification Tests (OVTs) (on 1 
Oct. 1977) 

(2) Start support of ARC Ground Data System (GDS) tests 
(on 1 Nov. 1977) 


(3) Start support of MOS Operational tests (in January 
1978) 

(4) Start DSN Initial Acquisition tests at launch minus 30 
days 

(5) Support prelaunch Operational Readiness Tests (ORTs). 

The initial DSN OVTs uncovered some problems such as 
the inability of the Telemetry Processor Assembly (TP A) to 
process Pioneer Venus uncoded telemetry, and several station 
and NOCC problems associated with the newly implemented 
high-speed data 22-bit error polynomial. These problems were 
discussed in detail in The DSN Progress Report 42-45 (Pioneer 
Venus 1978 Mission Support for April 1978). 

The above problems delayed the DSN readiness for GDS 
support by approximately one month; however, the principal 
obstacle in preparing for the Orbiter launch was severe prob- 
lems with the ARC simulation, telemetry, and command soft- 
ware. To expedite the development of the ARC software a 
series of single system “data flow” tests were initiated in 
January 1978, occasionally supported by DSS 11, 44, 62 and 
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^dIL-71; but the bulk of these numerous tests were conducted 
)etween ARC and CTA-21, 

'The DSN OVT schedule was also extended to keep the DSS 
:rew proficiency at an acceptable level, and these added OVTs 
vere also utilized as a data source to provide ARC with 
)pportunities to further develop -their telemetry and command 
oftware, ' 

Formal CDS tests started on 1 March 1978 and MOS tests 
n 20 March 1978; and the ORTs in May 1978, with no 
ignificant DSN problems. 


Pioneer Venus Orbiter Launch 

(1) The Pioneer Venus Orbiter Initial Acquisition study 
was published and sent to the Pioneer Acquisition 
station (DSS 44) and the back-up station (DSS 42) on 
15 April 1978 prior to the Initial Acquisition OVTs. 
The Initial Acquisition OVTs were conducted success- 
fully, with both stations confirming that station and 
JPL crews'scheduled to support the launch were fully 
trained. 

(2) All DSN launch-required facilities were green prior to 
launch. There was, however, a major antenna failure 
during a rain storm at DSS 44, with repairs completed 
only 24 hours prior to launch. 

(3) Predicts had to be regenerated at L-20 minutes due to 
. a "5 kHz change in predicted frequencies supplied by 

Ames. 

(4) After a reasonably smooth minus count the Atlas 
Centaur lifted off at 13:13:00.073 GMT. The launch 
to parking orbit was nominal with Main Engine Cutoff 
(MECO) 1 occurring 7.4 seconds before predicted 
time. The Centaur second burn was nominal with 
MECO 2 and spacecraft separation occurring 4.5 sec- 
onds before predicted time. Spacecraft ‘injection 
appeared to be very nominal with midcourse correc- 
tion calculations from preliminary data calculated to 
be about 3 meters per second. 

(5) During the minus count only two communication 
problems of any significance occurred. The spacecraft 
telemkry 7.2 kilobit/s HSD line was down for 24 
minutes which caused the Project to reschedule a 
command memory readout. The ACN 6.2 kilobit/s 
spacecraft data line to the Cape was red for 8 minutes. 

(6) During the plus count all ETR, STDN, and ARIA 
supported all data intervals essentially as planned. No 
major losses of data occurred. Mark events 1 through 
13 were read out. 


(7) ACN support for possible commanding was flawless, 
however, the PMOCC elected to not send the com- 
mand to start the science sequence because of ambigu- 
ous telemetry indications. This was the most signifi- 
cant post-launch problem and first became apparent 
when processing the post-separation ACN data by the 
PMOCC showed a tremendous number of spacecraft 
problems. At the same time the Hughes computer at 
the Cape, processing supposedly the same data, 
showed all was well on the spacecraft. This ambiguity 
in the spacecraft status caused the Project to not send 
the command from ACN to start the science se- 
quence. The first DSS 44 data processed by the 
PMOCC showed all was well and the science sequence 
was started from that station. 

This discrepancy in processing of ACN data was 
caused by an error m the high-speed data system 
Block Error Decoders (BED) located at Ames and was 
not isolated until one month after launch. The BEDs 
were relabeling high-speed data blocks found to have 
errors as error-free blocks. 

(8) The initial acquisitions by DSS 44 and 42 were text- 
book except for a momentary loss at the backup 
station (DSS 42) when the antenna ran into the DEC 
prelimits during autotrack on the acquisition aid an- 
tenna. DSS 44 acquired the uplink with only 5 Hz left 
of the planned uplink sweep at VCO level. 

(9) The only anomaly during the first pass was an indica- 
tion of the spacecraft rejecting a ground command 
that was never sent after DSS 62 AOS. This anomaly 
was later determined to have been caused by excessive 
radiated transmitter power causing the command 
detector (CMD) to interpret noise, as a command. This 
CMD was then rejected* due to checksum errors. 
DSS 44 radiated 500 watts (SCM) during the initial 
acquisition pass and DSS 62 acquired at 1 kW as 
planned. Subsequent to the above diagnosis, uplink 
power was reduced to 300 watts at all stations until 
launch plus 2 days when the power was increased 
(over DSS 11) to 1 kW. 

Magnetometer boom deployment, first precision 
maneuver, and high gain antenna spin-up were sup- 
ported without incident. Extensive science calibration 
related activities have also 'been continuing. 

(10) Two other major events after launch were the Trajec- 
tory Correction maneuver on 1 June and the first 
attempt at executing the High Gain Antenna and 
X-band calibration over DSS 11 and 14 on 8 and 9 
June. The X-band downlink signal level was 10 dB 
lower than expected. Subsequent testing showed the 
loss was due to bandwidth spreading resulting from 



the uplink command modulation being turned around 
into the X-band. 

Although this report has concentrated on problems en- 
countered, the overall launch support was considered very 
successful and the view has been generally expressed that this 
was the smoothest planetary launch to date. 


III. Pioneer Venus Multiprobe 

A. Launch Preparation 

Tlie Multiprobe launch preparation is a relatively small task 
compared to the Orbiter launch, as the stations involved are 
tracking the Orbiter daily and both spacecraft are very similar 
in the launch/cruise mode as far as the DSN is concerned. 

The only DSN tests planned are initial acquisition OVTs 
with DSS 42 and 44; and support of a single MOS test on 
12 July 1978; and two ORTs on 27 July and 2 August 
respectively. 

The only change from the Orbiter launch is that DSS 42 
will be prime, with DSS 44 as back up. The reason for this is 
that the spacecraft will rise and set in the DSS 44 West 
“Keyhole'* resulting in the DSS 42 rise earlier than DSS 44 
and also resulting in DSS 42 setting approximately one hour 
ten minutes later than DSS 44, permitting an overlap with 
DSS 62. 

B. Multiprobe Entry Plahning 

As a result of the four 8-hour Multiprobe Entry procedure 
development tests in March 1978 (discussed in The DSN 
Progress Report 42-45) conducted at DSS 14, the following 
procedures have been generated and will be used in the remain- 
ing Multiprobe Entry preparation activities: 

(1) A detailed step-by-step sequence of events covering the 
entire Multiprobe Entry portion of the mission was 
developed. 

(2) A first cut at an E-3 hour checklist (E = time of entry 
of the probes into the Venus atmosphere; nominally at 
an altitude of 200 kilometers) was developed. The start 
of this check at DSS 14 and 43 is coincidental with the 
-time at which the station configuration is finally frozen 
(i.e., no further rectification of hardware failures). The 


check consists of a switch-by-switch “preflight*’ fol- 
lowed by acquisition of the bus downlink signal by the 
five closed loop receivers and five open loop receivers. 
This signal is processed by the four (closed loop re- 
ceiver) telemetry strings and routed via HSDL to JPL 
for verification (see Fig, 1). Simultaneously, the closed 
loop receiver doppler is processed and transmitted to 
JPL and the four open loop receivers and Multi Mission 
Receiver (MMR) similarly acquire the bus signal, con- 
firmed on the Spectral Signal Indicators (SSIs), and 
record the open loop doppler and telemetry on the 
DRA and PPR recorders. The open loop signal recorded 
on the PPR recorders is replayed via the up-converter 
into the closed loop receivers and again routed back to 
JPL. This positively confirms correct configuration of 
all the closed and open loop data pathi Finally, the 
checklist requires the termination of the bus signal 
tracking and all other signal sources, the loading of new 
tapes, and the- loading of all required probe parameters 
and frequencies, etc. 

(3) The E-1 hour checklist is virtually a final recheck of 
switch positions and frequencies prior to probe acquisi- 
tion. 

(4) These procedures plus the extensive multiprobe en- 
counter prepass calibrations will be exercised in tests 
starting in July 1978 and continuing every other week 
to December 1978, specifically involving the station 
crews who will be supporting the actual encounter. 

The initial exercises at DSS 14 in March indicated that a 
16-man crew is required to operate at DSS 14 and DSS 43 
during the encounter, so the above tests are planned at a time 
when the prime crews are “back- to back,” enabling both 
crews to work an extra four hours in parallel along with the 
off shift personnel needed to produce the full complement 
of 16. 

C. ALSEP Tracks 

The ALSEP tracks, are being executed to produce VLBI 
data in order to check out the VLBI Wind Measurement 
experiment which involves station equipment, CTA-21 pro- 
cessing and MIT processing. The tracks involve various com- 
binations of DSS 14, 43, and STDN Santiago and Guam ti ach- 
ing stations. This series of test tracks was described in The 
DSN Progress Report 42-45\ and the next aiticle in this series 
will include a detailed status of the results of the tests con- 
ducted from February through July 1978. 
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Thh article reports on activities of the DSN Network Operations organization in 
support of the Helios Project from 15 April 1978 through 15 June 1978. 


I. Introduction 

This article is the twenty-second in a continuing series of 
reports that discuss Deep Space Network support of Helios 
Mission Operations. Included in this article is information on 
Helios-rs seventh perihelion, Helios-2’s fifth perihelion, and 
otlier mission-related activities. 

II. Mission Operations and Status 

The Helios- 1 spacecraft passed through its seventh peri- 
helion on 29 April 1978 at 06:13:30 Universal Time Coor- 
dinated (UTC). The perihelion pass was not covered by a DSN 
station, but rather by DSS 67/68 at Weilheim, Germany. The 
spacecraft was at a bit rate of 512 bits per second (bps) coded 
in Format 1 , configured in medium power mode and with all 
experiments on. All spacecraft subsystems and experiments 
were in excellent condition. 

Hell os-2 passed through, its fifth perihelion on 30 April 
1978 at 05:08 UTC, and as was the case witli HeIios-1, 
coverage was by DSS 67/68. The bit rate was at 256 bps coded 
Format 2. The spacecraft was in a medium power mode and aU 
experiments were on. All of Helios-2 subsystems and experi- 
ments appeared to be in excellent condition. 


On 3 June 1978, the Helios Project requested that the DSN 
64-meter stations use 20-kilowatt (kW) power for uphnking to 
Helios-2. The reason for this was due to the Sun-Eartli-Probe 
(SEP) angle decreasing below 4 degrees. Since 23 May 1978, 
the Inghly active Sun began to influence uphnk and downlink 
signals heavily, with ever-increasing frame deletions. On 3 June 
1978, two commands were lost over DSS 43 (Canberra, 
Australia) prior to the 20-kW uplink request, which prevented 
the spacecraft from configuring to the gap coverage read-in bit 
rate. This resulted in a 13-hour loss of data. In addition to the 
20-k\V mode, the Project will transmit aU critical commands 
three times, wliich has proved helpful m the past during similar 
conditions. Also, the telemetry bit rate was lowered to 128 
bps, wliich is the lowest bit rate capable of reading out the 
spacecraft memory to cover gaps between stations. The above 
configuration will exist until mid-July, 

An unexpected blackout entry occurred on 8 June 1978 
over DSS 63 (Madrid, Spain). Station 63 had trouble main- 
taining telemetry lock at 128 bps and the situation did not 
improve when the rate was lowered to 64 bps. DSS 67/68 
(Weilheim, Germany) was asked to acquire, and tlie results 
were the same as experienced at DSS 63. As a result, the 
project configured the spacecraft’s memory for a 22-day 
storage period assuming a symmetrical exit angle, and the 
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spacecraft was commanded to an 8 bps read-in mode on 
8 June 1978 at 18:07 UTC. The subsequent passes atDSS43 
and DSS 63 showed a very .degraded signal and proved “almost 
blackout” conditions at an SEP angle of approximately 3.1 
degrees. This region will be exited on 29 June 1978. 

Overall coverage of botli Helios spacecraft for tliis period is 
listed in Table 1. 

III. Special Activities 

A. DSN Mark III Data Subsystem (MDS) Update 

As reported m the last article^ DSS 11 MDS update was 
completed (Ref. 1) and the station placed under configuration 
control to support Helios flight operations on 26 April 1978. 
This completes the DSN MDS implementation. 

B. Support of On-Board and Ground Experiments 

The DSN 64-meter stations have continued practicing 
polaiimetry and Meteorological Monitor Assembly (MMA) 
procedures (Ref. 1). During the blackout period of Helios-2, 
discussed earlier in this article, Faraday rotation data (Experi- 
ment 12) will be collected during all 64meter Heiios-2 tracks. 


This utilizes the same procedures for polarimetry and MMA 
procedures mentioned above. These data will also be corre- 
lated with pulsar data collected by the 100-meter antenna at 
Effelsberg, Germany, on 12, 13 and 14 June. In addition to 
this period of high interest, the time of closest approach 
(smallest SEP angle) of 17, 18 and 1*9 June, is most critical to 
experimenters and good data is hoped for. 

Concerning Experiment 10 (micrometeoroid counter and 
analyzer), the experiment on Helios-2 was turned off. on 
5 June, after the experimenter suspected a problem, due to a 
current nse in Sensor-A of the experiment package. The date 
of Experiment 10 turn-on and checkout has not yet been 
determined at this writing. 

C. German Space Operations Center (GSOC) 
22-Bit-Error Polynomial Code (EPC) Conversion 

In October, 1978, GSOC is scheduled to convert to 22-bit 
EPC to conform to the DSN 22-bit configuration. The first 
test of the system, referred to by GSOC as the Mission Control 
and Computer Center (MCCC) bypass, was conducted on 
2 June 1978. Data- were successfully passed' from DSS 62 
(Cebreros, Spain) througli the liigh-speed switch at JPL to 
GSOC while in the 22-bit EPC. Future checkout of the system 
is planned prior to final conversion in October. 
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Table 1. Helios tracking coverage 


Month 

Spacecraft 

Station 

type 

Number 

of 

tracks 

Tracking time 
(hours, 
minutes) 

April 

Helios-1 

26-meter 

42 

252:07 



64-meter 

3 

13:56 


Helios-2 

26-meter 

39 

232:05 



64-meter 

21 

116:36 

May 

Helios-1 

26-meter 

55' 

333:41 



64-meter 

1 

5:06 


Helios-2 

26-meter 

. 22 • 

124:00 



64-meter 

29 

140:20 
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The currently deployed ranging systems are subject to a variety of waveform 
distortion errors and equipment vagaries which limit the accuracy of the range data used 
for Navigation and Radio Science, In light of this fact, this article presents arguments, 
both subjective and experimental, for increasing the accuracy of the currently deployed 
ranging system by the adoption of an approximately 1 MHz sine wave for the 
precision-defining signal for ranging. Inferences are also drawn for the desigfi of more 
precise and wider bandwidth ranging systems in the future. 


I. Introduction 

We have, in recent months, reviewed the instrumental 
limitations to range measurement accuracy. Two specific 
objectives were sought: devising cost-effective improvements 
to the current generation of ranging systems, and establishing 
guidelines for future ranging system development. Our analysis 
and supporting experimentation indicate that witiun the 
bandwidtli constraints of the current, or NASA Standard 
Transponder, a 1 MHz range code plus appropriate filtenng in 
the ground systems should improve the accuracy achievable to 
the order of 1 meter (6 ns). In the long term, further 
improvements to range accuracy can be achieved only by 
increasing the bandwidth of the entire end-to-end ranging 
system. Since many elements of a two-way ranging instrument 
are common with three-way or differential one-way instru- 
ments, the information we have obtained should also provide 
useful guidance to the design of these related systems. 

One unexpected result of our work was confirmation of the 
previously suspected sensitivity of the ranging system to errors 


caused by changes in Tracking Station equipment configura- 
tion- In our current test series, errors induced by changes in 
the uplink modulator were as large as 11 meters (70 ns). This 
type of error may be a prime contributor to the ranging 
system problems currently under investigation by both Viking 
and Voyager projects. 

II. Near-Term Improvements 

In the design of precise measurement systems, it is virtually 
imperative to utilize the entire available bandwidth with 
signals that are readily detected and processed in a way which 
directly contributes to improving the accuracy of the measure- 
ments. This is not done in the current ranging systems, in 
which botli the fundamental and third harmonic of the 
precision-defining range code fall withifi the passband of 
current transponders. While the third harmonic could be 
processed to enhance ranging accuracy, this processing is not 
done, nor is it readily implemented- With the current 
range-code detection schemes, the third harmonic contributes 
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nothing to measurement precision, and could be removed by 
filtering (Ref. 1). Its principal function in the current system 
has been to allow a numerically simple conversion from 
range-code correlations to estimated range. In actuality, the 
presence of the tliird harmonic creates the opportunity for 
severe distortion errors to arise, so that its removal would 
improve accuracy, though not precision. Phase-shift and filter 
effects vary between stations and between transponders, to 
cause configuration-dependent bias errors and waveform dis- 
tortion effects which are the limiting error sources of the 
current system. 

Let us suppose that the received ranging waveform consists 
of the fundamental and third harmonic from a square-wave. As 
noted above, we have previously shown that simply throwing 
away the third harmonic does not degrade the precision of the 
range delay estimate with respect to thermal noise. However, 
aside from ambiguities, it is easy to show that a range delay 
estimate can be constructed using this third harmonic alone 
which has exactly the same precision as the delay estimate 
derived from the fur^damehtal component. Since the funda- 
mental can resolve the three-fold ambiguities of the third 
harnionlc, and since their thermal noise errors are indepen- 
dent, tlie two delay estimates could be combined to effect a 
V^impro Yemen t in precision over either one alone. 

If, however, the transponder in question can pass the third 
harmonic of the conventional 500 KHz range code with little 
attenuation, it can also pass the fundamental of a 1 MHz range 
code with little or no attenuation. Doubling the range-code 
frequency provides a direct factor-of-two improvement in 
delay-estimate precision with respect to the conventional 500 
kHz range code as conventionally processed. It also provides a 
^/2 improvement in' delay-estimate precision over tlie best 
processing for the conventional range code. As the third 
harmonic of the 1 MHz code now does not pass through the 
transponder, further improvements require first increased 
transponder bandwidths, and subsequently increased band- 
width throughout the system. 

There is an additional advantage to the 1 MHz code which 
results because none of the code harmonics pass through the 
transponder. Professor L. Rauch has recently pointed out 
(Ref. 2) that nonlinear elements such as exist in the tran- 
sponder’s modulator and power amplifiers can induce pliase- 
sliifts in the fundamental component of a waveform which 
depend upon the phase of the harmonics entering the 
nonlinear element. With all harmonics removed from the range 
code prior to such nonlinearities, delay errors induced by the 
harmonics should be minimized. We anticipate that use of the 
1 MHz range code, coupled with filters in the ground receiving 
system which pass only the fundamental of this code, will 
greatly reduce sensitivity to DSN equipment variations which 


manifest themselves in changes m the harmonic structure of 
the waveform. 

Knowing that the third harmonic of the conventional 500 
kHz range code is passed by the transponders, it has been 
suggested that a range code of that frequency be used to gain 
an additional factor-of-1.5 increase in range precision. This 
gain is unlikely to be achieved for two reasons: First, the third 
harmonic is attenuated by 3 dB in several of the present 
physical transponders, and second the 1.5 MHz signal is now 
close to the transponder filter band-edge and thus more 
subject to phase-shifts due to environmental changes in the 
transponder. 

It is worth noting here that if a PN range code were used 
instead of square-wave, only those spectral lines below the first 
zero of the spectrum are at a phase which helps improve signal 
delay resolution, while those between the first and second 
spectral null (may) degrade resolution. This is exactly analo- 
gous to the behavior of a square-wave as discussed above, 
where the third harmonic, falling between the first and second 
spectral nulls, is nonbeneflcial in resolving range delay with a 
correlation-type estimator. When we earlier asked (Ref. 3) for 
the best unit-power waveform to be tracked by a delay-locked 
loop in terms of residual phase jitter, the answer was almost 
identical to a PN which had been filtered to include only the 
frequency region below the first spectral null. The length of 
that PN depended upon the SNR that was to be achieved by 
the delaydock tracker, which effected the band-limit con- 
straint for that system. 

The PN, however, is a waveform with many spectral 
components within the^ main lobe of its spectrum. Any 
nonlinearities in the ranging circuit will generate cross-modula- 
tion terms between these components which coincide with 
some other component. This creates the potential for wave- 
form distortion effects similar to those we have now, albeit at 
a lower magnitude. If we keep the waveform spectrally simple, 
i.e., a sine-wave which occupies nearly the full bandwidth 
available, the distortion effects induced by nonlinearities can 
be scrubbed off by filtering at appropriate places within the 
system. 

The conclusion we wish to draw from this can be 
summarized as follows: To get the best possible results for 
ranging, the available bandwidth should be completely utilized 
by those signal components which are most useful. For the 
near term, a 1 MHz range code will substantially improve the 
accuracy achievable with current transponders. Our recent test 
results, to be described next, reinforce this opinion. A future 
article will discuss other currently limiting error sources, and 
the anticipated effect of increasing the available ranging 
bandwidth. 
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III. Experimental Results 

Several experiments were recently conducted, to measure 
the so-called ranging “waveform-distortion” error and test the 
concepts previously proposed (Ref. 1) as a cure. These 
experiments utilized the Mu-II Ranging System (Ref. 4) 
installed^ at the Telecommunications Development Laboratory 
(TDL). The first of four groups of experiments determined the 
precision of the current ranging system. The remaining three 
experiments measured the improvement in precision gained by 

(1). filtering all but the fundamental frequency of the 500 kHz 
range code prior to demodulation, (2) utilizing a 1 MHz range 
code, and (3) using a 1 MHz range code filtered as in (1). 

The results of our investigations indicate that a filtered 
1 MHz code (3 above) reduces the waveform distortion to 
below 1 ns of error. This also increases the ranging system 
immunity to errors induced by equipment variation in 
receivers, exciters, and transmitters. 

A. The Setup 

Figure 1 illustrates the various experiment configurations. 
The Mu-II Ranging System range code modulated a carrier 
which was echoed either by a wideband zero-delay device or 
the proof test model of the MVM’73 transponder. Mu-II 
Channels 1 and 2 were respectively connected to 10 MHz IF 
signals from the Block IV and Block III receivers. One of a pair 
of passband filters was inserted into either the Block IV or 
Block III to Mu-II IF cable to provide a band-limited received 
range code. The filter characteristics are shown in Fig. 2 for 
the “500 kHz” filter and Fig. 3 for the “1 MHz” filter. 


B. Experiments and Results 

The general test plan consisted of utilizing both the zero 
delay device and the MVM73 transponder in conjunction with 
each of the folio wmg four test conditions: 

(1) 500 kHz code, urifiltered 

(2) 500 kHz code, with filtered IF signal 

(3) 1 MHz code, unfiltered 

(4) 1 MHz code, with filtered IF signal 

The data from tests (1) and (2) or test (3) and (4) were 
obtained simultaneously by connecting one channel of the 
Mu-II directly to the Block III (IV) IF signal and inserting a fil- 
ter in the Block IV (III) IF' signal attached to the other chan- 
nel. Comparison of Block III and Block IV data showed error 
differences on the order of tenths of nanoseconds. Because 
differences between the two Mu-II channels were of the same 
magnitude, the taking of data simultaneously is* reasonable and 
justified by minimizing the test time. It must be noted that 


Mu-II channel 1 exhibited an unexpected periodic error, which 
appears as periodic hash on the data graphs displayed below. 
From extensive testing we infer that this error is due to 
hysteresis in the channel 1 A/D converter and interaction 
between the three Mu-II coders. The hysteresis was predicted 
by Layland (Ref. 5). Coder interaction was discovered experi- 
mentally. Because these inseparable errors are small (<1 ns), 
they will be ignored in this article. 

^ Each of the tests described below involved stepping the 
local model of the received range code in 7 ns steps over a span 
of one quarter of the period of the highest frequency (^de 
(clock). At each step, the code phase measurement B is 
compared to the actual phase 0. The graphs are the error 
function 6 - d and are in units of nanoseconds. This is a direct 
display of the so-called “waveform-distortion” error. 

1. 500 kHz code, zero delay device. This configuration is 
identical to current DSN ranging (Planetary Ranging Assem- 
bly — PRA) zero-delay calibrations. Figure 4a is the error 
function without filtering. The peak-to-peak error is 18.5 ns. A 
dramatic improvement is shown in Eig..4b.. By adding the 
“500 kHz filter,” the error is reduced to 1.5 ns. 

2. 500 kHz code, transponder. The peak-to-peak error in 
Fig. 5a is about 7.4 ns. This is the error of the unfiltered 
ranging system. Adding the “500 kHz filter” (Fig. 5b) reduces 
the error to about 2,1 ns. This is a substantial improvement 
over the current ranging system. 

3. 1 MHz code, zero-delay device. Tests with a 1 MHz range 
code confirmed the expected improvements. Using the zero- 
delay device, a peak-to-peak error of 3.9 ns (Fig. 6a). Using 
the “1 MHz filter” reduced this error to 1.0 ns (Fig. 6b). 

4. 1 MHz code, transponder. This final test set further 
justifies the use of a 1 MHz code and filter. Without the filter 
(Fig. 7a), the error is 14.2 ns. When the “1 MHz” filter is used, 
the error (Fig. 7b) is only 0.6 ns. 

C. End-To-End Distortion Immunity 

The “end-to-end ranging system” includes the ranging 
machine, exciter, transmitter, receiver, spacecraft transponder, 
zero delay devices, and all cables, switches, connectors, etc,, in 
the ranging path. Each component in this path is capable of 
inducing, to varying degrees, range code distortions in both the 
frequency and' time domains. Jumps in DSN station range 
calibrations and spacecraft range residuals evidence the relative 
susceptibility of the range system to equipment vagaries. 

The exciter modulator is one major source of waveform 
distortion. Two standard DSN modulators were compared in 
terms of relative spacecraft delay versus modulation index. 
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The left side of Fig*. 8 shows an almost 70 ns worst-case 
disparity between modulators. Use of the 1 MHz code and 
filter reduced this to about 10 ns, A similar modulator 
measurement was made earlier (Ref 1) using only the 500 kHz 
code ,and “500 kHz filter.” The worst-case difference was 
10 ns without the filter and 3 ns with it. During the current 
series of tests, using the 500 kHz code with no filter produced 
the 70 ns disparity noted earlier. A 30 ns worst-case difference 
resulted when the “500 kHz filter” was used. These current 
numbers were easily repeatable to within 1“2 ns during the 
six-month period while we sought to explain the differences 
between the 500 kHz results now and in 1975. We now believe 
that these discrepancies are due to the extensive modification 
to TDL which occurred between the two test series. 

In addition, we strongly suspect that the numbers- previ- 
ously seen correspond only with the behavior of a superbly 
tended calibration laboratory which TDL was in i975. Our 
current measurements follow a period when TDL has been 
“mothballed” without continuous expert attention, and may 
in fact be more representative of the ranging system behavior 
to be expected in the operational DSN tracking stations. The 
delicate nature of the current range system is obvious. 


Utilizing a 1 MHz range code and IF filter/increases. system 
immunity to equipment variation. 


IV. Comments 

We have demonstrated that use of a 1 MHz range code and 
a narrow-passband IF filter can reduce waveform distortion 
error by a factor of as much as 20! Along witli this 
improvement in ranging precision, ranging accuracy is en- 
hanced by nullifying many errors caused by equipment 
vagaries. In our opinion, the IF filter and concomitant 
modification of ranging system software (and in the case of 
the Planetary Ranging Assembly, hardware) are the most 
cost-effective way of improving the current generation ranging 
system. 

The lessons we have learned here are alsoimportant for the 
development of the next generation ranging system. Clearly we 
need to increase the available bandwidth, and ^ we do, we 
should occupy that bandwidth almost completely witli a 
sinusoidal ranging' tone for precise definition of range. 
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Fig. 5. Phase error— 500 kHz code, transponder: (a) Unfiltered, 
(b) 500 kHz filter 


2L ]L ^ E. E 

8 4 8 8 4 

e 


Fig. 7. Phase error— 1 MHz code, transponder: (a) Unllltered, 
(b) 1 MHz filter 
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Fig. 8. Range variation caused by modulator change 



Fig. 6. Phase error — 1 MHz code, zero delay device: (a) Unfiltered, 
(b) 1 MHz filter 
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The Deep Space Network is currently engaged in establishing a new high-accuracy 
VLBI celestial reference frame. This article discusses the' present ‘status of the task of 
finding suitable celestial radio sources for construe ting this reference frame. To date, 564 
VLBI sources have been detected, with 166 of these lying within of the ecliptic- 
plane, The variation of the sky distribution of these sources with source strength is 
examined. 


I. Introduction 

The Deep Space Network (DSN) is currently engaged in 
establishing a new high-accuracy celestial reference frame com- 
posed of compact extragalactic * radio sources, principally 
quasars. These radio sources are observed By means of the 


technique of Very Long Baseline- Interferometry (VLBI), 
whereby two antennas separated'by distances of up to an earth 
diameter, simultaneously observe the same source (Ref. 1). 
This- article wilLdiscuss^the effort which has been undertaken 
to find a suitable set of celestial radio sources from which the 
VLBI reference frame may be constructed. 
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In the past, astrometric or geodetic studies which required 
the use of a celestial reference frame had to rely on an optical 
star reference frame (that is, a set of optical stars whose 
relative positions were well known). Such reference frames 
were limited to accuracies of about 0.1 arcseconds (5 X lO"*^ 
radians) in relative source positions due to the effects of 
atmospheric turbulence on optical observations and the rela- 
tive angular (“proper*') motions of the stars in the sky. Since 
VLBI sources are extragalactic and thought to be on average 
about 10^ times farther away than optical stars, proper 
motions should be negUgible. In addition, the VLBI technique 
should reach angular measurement accuracies approaching 
0.001 arcseconds (5 X lO"^ radians). Hence, VLBI radio refer- 
ence frames should possess accuracies which exceed those of 
optical reference frames by dne or two orders of magnitude 
(Refs. 2, 3). 

Of what use is such an accurate celestial reference frame? 
Some examples: 

(1) Calibrating Deep Space Probe Radio Metric Tracking 
Data — The DSN is currently developing an operational 
VLBI system to synchronize clock epochs and rates 
among the worldwide DSN spacecraft tracking stations 
and to monitor earth rotational irregularities in order 
to properly calibrate radio metric navigation data from 
interplanetary spacecraft (Ref. 4). In addition, the 
measured variations in earth spin rate and pole position 
are important clues to understanding the dynamics of 
the earth's crust and interior (Refs. 5, 6, 7). 

(2) Navigating Interplanetary Spacecraft — The DSN is cur- 
rently involved in a demonstration of the value of 
VLBI measurements of spacecraft position for improv- 
ing tile navigation of the Voyager spacecraft (Ref. 8). 
Computer simulations have shown that VLBI would 
provide significant navigational enhancements to future 
interplanetary missions as well (e.g.: Galileo, VOIR), 

(3) Measuring Planetary and Lunar Orbital Motions — 
Improved knowledge of planetary and lunar motion by 
means of observations of spacecraft in orbit about or 
on the surface of these bodies will provide higher 
navigational accuracy for deep space probes as well as 
useful information for studies of solar system 
dynamics. At present, JPL researchers are observing the 
motion of both the moon and Mars by VLBI tech- 
niques (Ref. 9). 

(4) - Monitoring Earth Crustal Motions — At JPL, project 

ARIES is currently demonstrating that VLBI can be 
used to monitor local crustal motions in order to study 
their, relationship to seismic events (Refs. 10, 11). 
Eventually, VLBI will be used to measure the relative 
motion of crustal blocks on a global scale. 


(5) Earth Surveying— The National Geodetic Survey is 
considering using VLBI techniques to assist in its effort 
to accurately determine the relative separations of 
widely spaced survey markers. 

With this justification, we can now consider what the 
logical steps in the development of a VLBI reference frame 
would be: 

(1) Conduct a sky search for celestial VLBI sources which 
might be useful. 

(2) Choose the “best” subset of these sources for forming 
the VLBI reference frame. 

(3) Determine accurate relative positions of these sources. 

Although it is conceptually possible to pursue these tasks in 
purely sequential order, -in practice we have performed these 
activities to a large extent in parallel. Nonetheless, to date, the 
most intensive effort has gone into searching for sources, and 
that is the subject of the remainder of this article. 


il. The Search Concept 

The search for VLBI sources first became an organized 
pursuit in the summer of 1974. Up to that point, perhaps 60 
celestial VLBI sources had been found by DSN efforts as well 
as astronomers at other institutions. It became evident that 
this number of sources would not suffice for future VLBI 
requirements, and moreover, that this was probably only a 
small Traction of the total number of VLBI sources in the sky, 
and hence, probably not an optimum set from which to form a 
VLBI reference frame. The search began slowly, but in the 
past two years, increased in intensity. Up to the present time, 
more than 30 individual observing sessions have been 
organized. 


III. What is a VLBI Source? 

The strength of a celestial radio source, or its total flux 
density, can be expressed as: 



over the 
source 


where 

= total flux density in Jansky 
(1 Jansky = 10"^^ watts/m^/Hz) 
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6,<p - orthogonal angular celestial coordinates mea- 
sured from a designated source center in radians 

B(6^<1>) = the angular brightness distribution of the source 
in Jansky /steradian 


The VLB! strength of a celestial radio source, or its correlated 
flux density, can be expressed as (Ref. 12): 


?; = ffsiese 


ri(2nKu$ + v0) 


d0d(l>\ 


over the 
source 


where 

= correlated flux density in Jansky 

u,v= the components of the baseline vector separating 
the two observing antermas, along $, 0 plane of sky 
coordinates, in RF wavelengths 

So, for a particular observation, the correlated flux density 
is merely the ampHtude of one component of a complex 
two-dimensional spatial Fourier transform of the angular 
brightness distribution of the source. If the source angular size 
is significantly smaller than the spatial wavelengths of the 
transform component, 1/w and 1/v, then it is easy to derive 
that the correlated flux density is equal to the total flux 
density. If most of the power emanating from a source comes 
from angular areas that are large compared to the spatial ^ 
wavelengths, then virtually none of the total flux density of 
the source appears as correlated flux density. So, to observe a 
significant fraction of the total flux density as correlated flux 
density, much of the total flux density must come from 
regions smaller than 1 /m and 1/v. 

Our baselines lengths were about 7 X 10^ RF wavelengths, 
and so the spatial wavelengths, l/ii and 1/v, w^e generally on 
the order of a few milliarcseconds. Hence, the radio sources we 
detect as VLBI sources have a significant amount of their total 
flux density contained in components that are no more than a 
few milliarcseconds in angular extent. Only a small percentage 
of celestial radio sources have such tiny strong components. 
More typically, radio sources are many seconds of arc to 
several minutes of arc in angular extent. 

What types of celestial sources have such small compo- 
nents? Within our own galaxy there are not many VLBI 
sources of interest to us. A few molecular clouds emit maser 
radiation from compact regions (Ref. 13). However, because 
the molecules radiate at only discrete spectral lines, are 
extremely time variable, and aie relatively nearby, these 
sources are not the object of our search. Since stars are usually 
very weak radio sources, few stars have been detected with 


present VLBI systems (Ref. 14). The only confirmed, reliable, 
and wideband VLBI source of reasonable strength within our 
own galaxy is the galactic center source (Ref. 15). 

Exterior to our galaxy, VLBI sources within our limits of 
detectability faU into three categories: 

(1) Quasars— The most distant class of extragalactic 
objects. Optically, they are stellar in appearance, with 
emission and absorption spectral lines shifted far to the 
red. By current cosmological thought, the large red- 
shifts are indicative of great distance, generally billions 
of light years. The extraordinary power output per unit 
volume exhibited by quasars is not explained by 
current theories. A prominent suggestion is that the 
large energies are released by matter falling into massive 
black holes, perhaps as large as 10^“ 10^^ solar 
masses. Ground based VLBI studies indicate that not 
only do some of these objects possess angular radio 
structure at the milliarcsecond level, but that this 
structure can be time-variable over a few months with 
apparent internal velocities exceeding the speed of light 
(Ref. 16). Most of the VLBl sources we find are 
identified as quasars, when optical identifications exist. 

(2) BL Lacertae Type Objects — These radio sources appear 
to be similar to quasars, but their atomic emission lines 
are nonexistent or weak. They are generally somewhat 
more compact and time-variable than quasars (Ref. 17), 
Astronomers feel they may be the bright cores of 
distant elliptical galaxies. 

(3) Galactic Cores — Many galaxies possess an active 
central region from which radio emission emanates. 
VLBI observations have shown that galactic cores can 
contain angular structure at the milliarcsecond level 
and, in one case, structural changes which indicate 
apparent internal velocities exceeding the speed of light 
(Ref. 16). The VLBI features are often aligned with 
larger scale features of the source, which suggest that 
both small and large scale radio structure may have 
originated from similar explosive events in the core at 
different epochs. Of perhaps more terrestrial interest is 
the fact that our own 'galaxy has an active VLBI core. 

The noted similarities of quasars, BL Lacertae objects, and 
galactic cores may be an indication that all three are merely 
different classes of a single type of celestial phenomenon, 
perhaps at different evolutionary stages. 


IV. Selection of Candidate Sources 

If one were to examine the entire sky for VLBI sources to 
the level of sensitivity we desire, the task might take several 
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decades. Hence, our search was relegated to observing known 
celestial radio sources to see which of these might be VLBI 
sources. 

In order to keep from searching all known celestial radio 
sources for VLB! components, we have used existing informa- 
tion on radio sources as clues to source size. The following 
criteria proved useful in determining which sources might have 
structure confined to the miUiarcseconddevel: 

(1) Spectrum — Power does not fall off with increasing 
frequency as rapidly as more typical radio sources in 
the range 1-10 GHz. An even stronger indicator is if 
the spectrum is flat, sloped upward, or peaked in this 
range. Of all the indicators of milliarcsecond structure, 
spectral properties have proved the most useful. 

(2) Variability — Identified as radio or optical variable. 

(3) Optical Identification — Identified as extragalactic. 

(4) Size Limits— A short baseline interferometer can 
determine if a source has components smaller than its 
resolution limit (resolution does not exceed a few 
arcseconds). 

(5) Interplanetary Scintillation — Radio sources that show 
scintillations when viewed throu^ the solar corona 
possess components with angular sizes <0.1 arcsecond. 

(6) Detection in High Frequency Sun^eys — A larger 
percentage of sources detected in surveys between I 
and 1 0 GHz are VLBI sources than is the case at lower 
frequencies. 

We searched the literature for information on celestial radio 
sources. This was not a simple task, as the information is 
scattered through hundreds of articles. Nor is the literature 
static. Many of the high frequency radio surveys have only 
been performed during the last few years, and some of the sky 
is stUl unsurveyed at high frequency. Fortunately, a few high 
frequency radio surveys (e.g.: Parkes, NRAO-Bonn, Ohio 
State) provide data on a large fraction of the sky. In general, 
we chose as candidates only those sources whose total flux 
densities were greater than 0.7 Jansky at 2.3 GHz. 
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Australia or Goldstone-Spain baselines as pairs of antennas on 
the Spain- Australia baseline cannot simultaneously observe a 
very large area of sky. Generally, 64 - 26-meter antenna pairs 
were used, but occasionally 64 - 64-meter combinations were 
substituted. All the observations were performed at 2290 
MHz. The receiver chain generally consisted of an S-band 
travehng wave maser followed by a special phase, stable S-band 
VLBI receiver which converted the signal to an IF of 50 MHz. 
The Mark II VLBI recording system, which was developed by 
the National Radio Astronomy Observatory, then recorded a 
1.8 MHz data bandwidth by digitally sampling at a4-Mbs rate 
(Ref. 18), Digital sampling and phase stability of the receiver 
chain were controlled by rubidium or hydrogen maser atomic 
clocks. System temperatures were measured at both antennas 
for each source so that the data could be properly calibrated. 


VL Data Processing 

Matching tapes from the two antennas were then cross- 
correlated on a special hardware/software computer at 'the 
National Radio Astronomy Observatory in Charlottesville, 
Virginia. Computer manipulation of the output of this 
correlator yielded the correlation constant, p, for each 
observation, or the fraction of bits on the two tapes that were 
correlaled. Correlation constants were then converted into 
correlated flux densities, 5^, by means of the expression: 



where 


= correlated flux density in Jansky 

T = the measured system temperature at antenna i 
inK 



the inverse sensitivity of antenna i in Jansky /K 
(i.e., how strong does a source have to be in 
total flux density (Jansky) to raise the system 
temperature IK) 


To date, our hterature searches have yielded about 1100 
celestial radio sources which might prove to be VLBI sources. 
Of these, we have performed VLBI observations on almost 900 
sources, and have detected 564 as VLBI sources. 


V. The Observations 


The uncertainty in the correlated flux density measure- 
ments due to random noise effects may be expressed as: 


a. = 6.3 X 10^ 

C 


T T 1 
^12 A 




The observations were all performed with pairs of antennas 
within the DSN which were separated by intercontinental 
distances. In practice, this meant observing on the Gcldstone- 


where 

e, = the antenna efficiency of antenna i (dimensionless) 
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D. - the diameter of antenna i in meters 

5 = the recorded bandwidth in Hz 
f- the coherent integration time in seconds 

Most sources were observed for at least three minutes, but 
for data processing the observations were broken into one 
minute segments. From this we see that the 5 a detection limit 
for this search for VLBI sources was about 0.1 Jansky. In 
addition, one would expect the random uncertainty in 
detected source strength to be about 0.02 Jansky. However, in 
practice, systematic errors at the 5 or 10% level dominate the 
random contribution for most sources. 

If a priori source positions are in error, the tapes must be 
correlated over a range of relative tape delay and delay rate 
offsets in order to detect a VLBI source. Appropriate searches 
in these parameters were performed so that the sky was 
completely searched within 0.5 arcminutes of all nominal 
source positions. Almost all a priori source position errors 
should be covered by this degree of position searching. In 
addition, once detections were achieved, the measured delay 
and delay rate offsets allowed positional errors to be measured 
to about one arcsecond, a good level from which to start 
building a VLBI reference frame. 


VII. Results 

The reduction of our data is not yet at a complete stage. 
Accurate source strengths and positions have been calculated 
for less than half of our detected_sources. For the remainder 
we have strengths accurate to about 20%. Nonetheless, it is 
useful at this stage to study the general sky distribution of 
detected sources as a function of source strength. In later 
articles we will publish catalogs of detected sources. 


VIII. The General Sky Distribution 

Figures la, lb, and Ic display the general sky distribution 
of detected VLBI sources as a function of VLBI source 
strength. The plots show the entire sky, with ±12 hours of 
right ascension, or celestial longitude, being displayed hori- 
zontally, and ±90^ of declination, or celestial latitude, being 
displayed vertically. These maps are equal area representations, 
so that local spatial densities in one area may be compared to 
local densities in another. 

It should be noted that displaying these source distributions 
as a function of source strength is fraught with-problems, since 
the strengths of the sources may be uncertain for a number of 
reasons: 


(1) Due to incomplete data reduction, our source strengths 
are only accurate to about 20% for more than half of 
the sources. 

(2) The total flux density and structure of these sources 
are often time-variable over months or years. 

(3) The correlated flux density may change as a baseline 
turns with earth rotation. This is due to the fact that 
changes in the projection of the baseline into the plane 
of the sky in the source direction can cause different 
spatial frequencies of the source’s brightness distribu- 
tion to be sampled. 

(4) Differences in correlated flux density between different 
baselines may also be apparent due to the samphng of 
different spatial frequencies. 

Hence, in displaying VLBI sky maps as a function of source 
strength, our gradations in source strength should be kept 
quite coarse. Even then, the categories are not crisp. However, 
although the plots may not show the true strength of a 
particular source aLa certain epoch and on a certain' baseline, 
they do indicate the general densitites of sources as a function 
of strength. 

Figure la shows the global distribution of all 564 detected 
sources. The lower limit of correlated flux strength is about 
0.1 Jansky. Several things are evident from this plot: 

(1) In general, the distribution is dense and rather uniform. 

(2) Pairs of DSN antennas cannot commonly view the sky 
below ^^^45^^ declination. 

(3) The north polar area is more sparsely populated. This is 
directly due to a lack of sufficient high frequency 
surveys in this region to easily identify candidate 
sources to observe. 

(4) The areas near the plane of the galaxy (not shown) are 
more sparsely populated. This is particularly evident 
between ±30° declination at right ascensions of about 
“5 and ±7 hours. This is due to the fact that high 
frequency radio surveys often skip the galactic regions 
and that interstellar charged particles near the galactic 
plane can cause intrinsically small sources to be 
scattered to larger angular sizes. 

We should point out that this is not nearly a complete map 
of the VLBI sky at the level of 0.1 Jansky. Since one of our 
criteria for choosing candidate sources was to generally only 
keep sources with tot^ flux densities greater than 0.7 Jansky, 
we have missed a large number of very compact sources with 
lower total flux densities. 
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Figure lb shows the distribution, of 184 sources with 
correlated flux densities greater than 0,5 Jansky, This is also a 
rather dense and even distribution, and is a more complete 
representation of the VLBI sky than was the 0,1 Jansky map. 

The DSN Block 1 VLBI System is being developed to 
calibrate radio metric navigation data from deep space probes 
(Ref. 4). The sensitivity of this system to VLBI source 
strength is 0,5 Jansky. In order to conserve antenna time, this 
system is intended to produce the appropriate calibrations 
with only 1 or 2 hours of observing on each of two baselines. 
During that short time span, a single pair of DSN antennas will 
see only a small segment of sky in common. In order to 
properly estimate the desired parameters, there must be a good 
distribution of sources within this small common area of sky. 
This must be true at any time of day. Hence, this system 
requires a source catalog of about 100 sources. The number of 
sources we have found with strengths greater than 0.5 Jansky 
is nicely matched to the number of required sources. 

Figure Ic displays the distribution of 50 sources with 
correlated flux densities greater than 1,0 Jansky. This is also 
quite a complete representation of the VLBI sky, but the 
distribution is no longer very dense. Note that less than 10% of 
the total number of detected VLBI sources have correlated 
flux strengths greater than 1.0 Jansky. 

The falloff in number of sources with increasing source 
strength is shown more clearly in the histogram of Fig. 2. The 
number of sources at source strengths less than 0.5 Jansky is 
smaller than the true sky distribution, due to our candidate 
selection criteria and our approximate sensitivity limit of 0.1 
Jansky, At the high end of the source strength scale, we see 
that only 10 sources are stronger than 2.0 Jansky and none arc 
stronger than 4.0 Jansky. The middle range of this histogram, 
say, 0.5 to 2.0 Jansky, can be approximately described by the 
expression: 

N=50S 

c 

where 

N = the number of sources stronger than a correlated flux 
strength of Jansky 

If this relationship holds at lower -source strengths as well, then 
we might ‘expect to find perhaps 4000 sources with correlated 
flux strengths greater than 0.1 Jansky, or seven times our 
present number of detections. 


IX. The Ecliptic Distribution 

Figures '3a, 3b, and 3c display the distribution* of VLBL 
sources within 10° of the ecliptic plane. These maps are 
plotted on rectangular grids rather than equal area representa- 
tions, The ecliptic plane appears as a sinusoidal-like trace 
across these plots. 

The ecliptic plane is important in solar system astrometric 
VLBI work, for it is here that planets and spacecraft mostly 
travel. The motions of these* objects may be accurately 
monitored against the background of VLBI sources by means 
of the technique of Differential VLBI (AVLBI). With AVLBI, 
a spacecraft and an extragalactic VLBI source can be simultan- 
eously observed on the same VLBI baseline. When the VLBI 
measurements of the two sources are subsequently differ- 
enced, most error contributions cancel to a high degree, 
providing an accurate measure of the angular separation of the 
two sources. Since the remaming sizes of many of these errors 
after cancellation are proportional to the angular separation of 
the two sources, a high density of extragalactic VLBI sources 
along the ecliptic plane is desirable, so that a source close to 
the spacecraft is always available. It also follows that opera- 
tional AVLBI systems will probably be required to work at 
more sensitive levels than is required for the operational 
system which calibrates radio metric navigation data, so that 
spacecraft /extragalactic VLBI source separations can be 
reduced. Ideally, one might desire as many as 100 sources in* 
an ecliptic reference frame, so that a spacecraft would never be 
more than a few degrees from a useful source. 

Figure 3a shows all 166 detected sources within 10° of the 
ecliptic. The distribution is rather dense and uniform except 
for the two areas where the galactic plane crosses the ecliptic 
plane at about 6 and 18 hours of right ascension. It may be 
possible to increase the density of sources in the galactic 
regions by observing at X-band rather than S-band, as the 
scattering effect is thought to be inversely proportional to 
approximately frequency squared. Also shown are the posi- 
tions along the ecliptic of the critical events in the Voyager 
and Galileo missions where AVLBI might provide significant 
navigational advantages. 

Figure 3b shows the distribution of 51 sources within 10° 
of the ecliptic plane which have correlated flux densities 
greater than 0.5 Jansky, The distribution now becomes 
relatively sparse in some areas. Figure 3c shows the distribu- 
tion of the 14 sources with strengths greater than 1.0 Jansky. 
The distribution is now very sparse as demonstrated by the 
fact that no sources appear anywhere near the Voyager 
encounters of Jupiter and the whole Galileo^ Jupiter orbit 
phase is almost quasar-free. 
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X. Summary 

We are engaged in a task of searching for a suitable set of 
celestial radio sources from which a high precision VLBI 
reference frame may be constructed. The observational por- 
tion of this program is nearly completed. We have shown here 
sky distributions of the 564 detected VLBI sources greater 
than about 0.1 Jansky and the subset of 166 of these sources 
which lie within 10® of the ecliptic. These distributions are 
reasonably dense at the level of 0.5 Jansky and below. The 
DSN operational VLBI system for calibrating radio metric 


navigation data will utilize sources stronger than 0,5 Jansky. 
An operational AVLBI system would probably have to employ 
weaker sources. 

Significant data reduction and some observing remains to 
complete the task of searching for VLBI sources. Published 
catalogs of VLBI sources, source positions, and correlated flux 
densities will follow this article. The tasks of choosing the 
“best” of the available VLBI sources and determining accurate 
relative positions of the chosen sources are now receiving 
mcreased emphasis. 
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Fig. 3. VLBI sources within ±10° of ecliptic; (a) 166 sources >0.1 Jansky, (b) 51 sources 
>0.5 Jansky, (c) 14 sources >1.0 Jansky 
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A million ) charineU 3Q0 MHz bandwidth, distal spectrum analyzer is being 
constructed at the Jet Propulsion Laboratory. The desigfi, fabrication, and maintenance 
philosophy of the modular, pipelined. Fast Fourier Transform (FFT) hardware are 
described. The spectrum analyzer will be used to examine the region from I A GHz to 
26 GHz for Radio Frequency Interference (RFI) which may be harmful to present and 
future tracking missions of the Deep Space Network The design will have application to 
the Search for Extraterrestrial Intelligence (SETIj signals and radio science phenomena. 


I. Introduction 

A million channel digital spectrum analyzer with 300 MHz 
of bandwidth is under construction at the Jet Propulsion 
Laboratory. The 300 MHz bandwidth was chosen to take 
advantage of a 300 MHz bandwidth K-band traveling wave 
maser which is used with parametric up-converters to cover the 
frequency range from 1.4 GHz to 26 GHz. The purpose of the 
spectrum analyzer is to detect and identify radio frequency 
interference which may be harmful to present and future 
spacecraft tracking missions of the Deep Space Network. 

The block diagram of the spectrum analyzer is shown in 
Fig. 1. The RF system consists of an antenna, parametric 


up-conyerter, rnaser, and receiver with 300 MHz bandwidth IF 
output. An analog filter bank is used to separate the 300 MHz 
into 32 complex channels of 10 MHz bandwidth each. A 
pipelined, decimation in frequency FFT is used to process two 
of these 10 MHz channels simultaneously. The power spec- 
trum is obtained by squaring tire real and imaginary parts of 
the complex spectrum. The power spectrum is accumulated 
for a number of spectra to reduce the data bandwidth to a 
manageable value for input to a general purpose computer. 

An extensive computer simulation was performed to 
determine the optimum hardware implementation to support 
the 60 dB dynamic range required. As a result of this 
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simulation, the hardware is implemented using 8 bit A/D 
converters, 12 bit memories in the first four stages, 16 bit 
memories in the remaining 11 stages, and 16 bit, fixed point, 
hard scaled calculations in all stages. 


II. Analog Filter Bank 

The bank of analog filters, shown in Fig. 2, is used to break 
the 300 MHz total bandwidth into 32 channels of 10 MHz 
bandwidth each. The input IF signal is distributed to the 32 
complex mixers togetlier with the output from a comb 
generator. This comb generator supplies the 32 local oscillator 
frequencies separated by 10 MHz. The complex mixer is made 
of two mixers and a 90® phase shifter between the local 
oscillator input and one of the mixers. The imphase (real) and 
quadrature (imaginary) outputs of the mixer each pass through 
a 5 MHz low pass filter to an 8 bit A/D converter. These 10 
MHz A/D converters are now readily available at relatively low 
cost because of their use in digital conversion of TV signals. 

The low pass filters are 3 dB down at the band edge, 
resulting in aliasing between adjacent 10 MHz channels. This 
aliasing is compensated for by a postprocessing algorithm in 
the general purpose computer. 


III. FFT Block Diagram 

The FFT, shown in Fig, 3, consists of 15 pipelined stages 
(Ref. 1), each composed of a memory unit and a “butterfly” 
arithmetic unit. Only three types of modules are used in the 
entire FFT. The memory modules used for the first four stages 
have a maximum capacity of 16K complex words of 12 bits. 
The other 1 1 memory modules have a maximum capacity of 
1024 complex words of 16 bits. The same 16 -bit arithmetic 
module type is used in all stages. 

The memory modules are programmed with a dual-in-line 
header to provide the appropriate delay and trig coefficients 
for each stage in the pipeline. 

The input to the FFT uses the “Biplex” method (Ref. 2) to 
simultaneously process two independent 10 MHz channels in 
the pipelined architecture FFT. This method results in the 
full-time utilization of all memory and arithmetic elements. 
The complete spectrum analyzer contains 16 of these dual 
15-stage pipelined FFT’s. 


IV. Memory Units 

A block diagram which is common to both types of 
memory umts is shown in Fig. 4. A memory unit is composed 
of two delay memories and multiplexers which allow straight 
through or crossed input-output connection as required in the 
pipelined algorithm. The memory unit also contains the trig 
coefficient generator. A delay equalization stage is incorpo- 
rated to compensate for different timing delays in the two 
paths through the arithmetic unit. 

The differences between the two types of memory units 
concern the size and type of delay memory and type of trig 
generator. 

The first four memory units, called 16K max on the FFT 
block diagram (Fig. 3) use CCD memories with a capacity of 
16K complex words of 12 bits. They are implemented by 
multiplexing four Fairchild CCD461A integrated circuits to 
obtain 10 MHz bandwidth. The remaining 11 units, called IK 
max, use random access memories (Intel 2125AL) with a 
capacity of 1024 words of 16 bits; 

The trig generator in the first four stages uses differences 
between successive coefficients, stored in ROM, to calculate 
the coefficients. This technique is usable because the decima- 
tion in frequency algorithm requires trig coefficients in the 
order of increasing angle. Storing differences rather than 
coefficients results in smaller ROM size and reduced total 
integrated circuit count. The trig generator in the remaining 1 1 
stages actually stores the coefficients in ROM. 

The trig coefficients used are always negative numbers so 
that -1 is included in the 2’s complement number field. 
Appropriate sign changes are incorporated in the arithmetic 
unit to allow this convention. 


V. Arithmetic Unit 

The FFT radix 2 butterfly arithmetic unit is shown m 
Fig. 5. The complex adder/subtractor butterfly is placed in 
front of the complex multiplier in the decimation in frequency 
algorithm. The adder/subtractor operates oh 16 bits of input 
data to deliver 17 bits of output. The output is scaled and 
rounded to retain the 16 most significant bits. The adder is 
implemented with the 74S283 and the subtractor with the 
74S381. 

The complex multiplier is composed of four real multipliers 
followed by an adder (74S283) and subtractor (74S381) to 
combine the partial products. The real multipliers are imple- 
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men ted with the TRW MPY-16A. Two of these multipliers are 
connected in parallel and multiplexed to obtain a 10 MHz 
multiply rate. This is siniple because of the input and output 
registers contained within the MPY-16A. The complete com- 
plex multiplier contains eight of the MPY-16A’s. A fractional 
multiply is performed, and the 16 most significant bits are 
retained. The internal circuitry of the MPY-16A is used to 
round the result. 


VI. Fabrication 

A dual 32K point FFT prototype is now under construc- 
tion at JPL. The prototype modules are constructed on wire 
wrap boards approximately 15 by 35 cm. 


VII. Maintenance Philosophy 

Testers are being designed to completely exercise the logic 
of the memory and arithmetic units of the FFT. These testers 
are used in manufactunng and depot level maintenance. 

In the field, digital test signals, instead of the A/D converter 
outputs, can be injected into the FFT. The output of any 
module can then be compared to the expected output which is 
obtained from the general purpose computer. The failed unit is 
isolated by this technique and repaired by replacement. Failed 
modules are returned to a depot for repair. Since there are 
only three module types in the FFT, only a small inventory of 
spares is required. All diagnostic tests are conducted at full 
clock rate. 
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Pig. 1. Spectrum analyzer block diagram 



Fig. 2. Analog filter bank 
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A theory is given which allows one to obtain the optimum acquisition sweep strategy 
of a PN code despreader when the a priori probability density function is not uniform. 
This theory has application to psuedo noise spread spectrum systems which could be 
utilized in the DSN to combat Radio Frequency Interference (RFIf In a sample case, 
when the a priori probability density junction is Gaussian, the acquisition time is reduced 
by about 41% compared to a ''uniform sweep'' approach. 


I. Introduction 

The acquisition circuitry of a despreader (a PN code 
acquisition and tracking system) is commonly designed so that 
complete passes are made across the code range uncertainty, as 
shown in Fig. la, during the initial search for the code epoch. 
This search, which is commonly implemented fay retarding one 
half a chip at a time, then integrating and comparing to a 
threshold (Fig. 2), continues until the signal is acquired. This 
scheme is efficient when the a priori location of the signal in 
the uncertainty region has a uniform probafaiiity density 
function. However, when the a priori density function is 
peaked, it is more likely to find the signal in the peaked region 
than elsewhere, so the full sweep approach may not be the 
best one. 


^Portions of this work were performed for theTDRSS Project at TRW 
Systems. 


This article is concerned with a method that allows one to 
determine the optimum sweep pattern to minimize tlie 
acquisition time, while achieving a required probability of 
signal detection, for a given a priori probability density of the 
signal location. The calculation is carried out for a Gaussian a 
priori signal location probability density function as illustrated 
in Fig. lb. However, the approach is general so that it can be 
applied to any given a priori signal location probability density ‘ 
function. 

The basis of this method relies on the fact that any 
meaningful statistics (Ref. 1, for example) of acquisition time, 
which is the time required to search the code until acquisition, 
depends directly upon the number of chips (code symbols) to 
be searched. Therefore searching where there is tlie greatest 
likelihood that the signal will be found first reduces the total 
number of positions to be searched. 
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II. Probability of Detection After N Sweeps 

Consider a symmetric search centered at the mean of a 
symmetric, unimodal, a priori probability density function as 
shown in Fig. Ic for the case of N= 3 sweeps. Denote as 
the probability of acquiring at the end of TV' sweeps. Typically 
Qjy would be 0.9 or 0.5, for example, LetXp L 2 , L3 . . . , 
denote the lengths (in number of cells) of the code range 
uncertainty to be searched during the N sweeps, and assume 
that ^ L^t p(x) be the a priori 

probability density of the location of the signal. Further, let 
denote the event that the signal is not detected in any one of 
the first / sweeps over regions with lengths . . . , 

Also, we shall use the notation 5 q to denote the event that the 
signal is not detected with zero sweeps, which is of course a 
sure event. It is clear that the conditional probability density 
of the signal location x, given that no sweep has yet been 
made, is equal simply to the a priori density p(x), i.e., 


Substituting (3), (4) into (2) we thus have 


p(x|5^) = 


(l-P^)p(x) 

l-PaP(Li) 

ifxdL 


(6) 


This conditional density is illustrated in Fig. 3b. It is easy to 
sliow that 



P{x\Si )dx 


(1 - Pg)PiL^) 1-P(LQ 

l-P^P(L^) l-P^P(L,)‘^ 

(7) 


p(x\S^) = p(x) 


( 1 ) 


This density is illustrated in. Fig. 3a. Suppose no signal is 
detected during the first sweep over , i.e., suppose the event 
Sj^ has occurred. The conditional density p(x|5'j) is equal to, 
by use of Baye’s rule, 


p(jc|5,) = 


( 2 ) 


The joint probability density of x and the' event. 5'^ is thus 
given by (from Eqs. 2 and 6): 


p(x,5^) = p(xl5^)?(5^) = 


(l-?^)p(x) ifxel^ 




ifx^L^ 


( 8 ) 


since P(i5j) = 1 - F^P{L^). The joint probability of and 
X eLi is given by 


In (2) the conditional (or a posteriori) probability density 
p(i5jx), is clearly given by 


P(Si\x) = 


l-P^ if xeLj 


1 


(3) 


where P^ is the probability of detection given the signal is 
there. The notations x eL^ and x denote the fact that 
the location of tlie signal is within the set or not in Zr j , and 
P(Si) is the probability of the event : 


P(S^) = 1 - P^P{L^) 


(4) 


P{L^, S^) = / p(x, S^)dx = (1 - P^)P(L^) (9) 

*7 -Xj /2 


Since ^ tile joint probability of Sj and ';ce 
(Z,2 “ 7,j) is given by (from Eq. 8) 


p(L, 


'L^/2 


p. 

-L,.S^)= P(x,S^)dx+ P(x,S^)dx 


^ P{L^~ L^) = PiL^)-PiL^) 


( 10 ) 


where P(Ly) denotes the probabUity that the signal location x 
is witliin the set £ j : 


^ij/2 

7’(£i) = I p(pc)dx 

J-LJ2 


(5) 


Further, let P,, i = 1,2, , N denote, respectively, the 

probabihties that tlie signal is acquired during the sweep, 
but not in the first, second, 1)^’’ sweeps. If we can 

compute these P,-’s, then Qjy is clearly equal to the sum 




N 


( 11 ) 
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It is clear that the probability of acquiring in the first sweep is 
given by 


= P^PiL^) ( 12 ) 

The probability P^ by definition, then, is equal to the joint 
probability 


P^ = Prob [acq. in second sweep, fail to acq. in first sweep] 


= i>^ [P{L^) - P(ii)I +P/1 - P^)P(L^) (13) 


In computing (13) we have used Eqs. (9) and (10). 
Following identical arguments one obtains the joint prob- 
ability density 

( 0 if 

p (x, S^) = \ (1 ■ if ^ e (L^ - ij) 

VPOc) ,\ix4L^ (14) 


It is then easy to compute the a posteriori probability 
density function after two sweeps: 




{1-Pdfp(?^) 

P(S^) 

{I- Pa)p{x) 

PiS^} 

P{x) 

P{S^) 


xeL^ 

x^Ly.xeL^ 

xiL^ 


(15) 


This density function is shown in Fig. 3c. From (14) one 
obtains, forij >^5 


P3 = P^P(i3 - P3. 5,) + P^P(i, - L^, 5,) + P^P(ij, 5,) 


In general, P. can be written for z = 1, 2, . . . , as. 

o-i’/-"W-np„_i)] (17) 

n-l 


withP(io) = 0. When (17) is substituted into (11) we obtain, 
after some simple algebra: 


N 


N 




1=1 


k=i 


pj^k) 


(18) 


The fact that we can write in the form of Eq. (17) is 
intuitively clear. We can argue through the case of , as given 
in (16). The region of search during the third sweep can be 
divided into three nonoverlapping regions = P3 ~ L 2 , 
and Xj . By the time the third sweep is initiated, the region-Xj^ 
has already been- searched' twice. Hence the joint probability 
that signal is detected in this region during the third sweep but 
not during the first two sweeps is P^(l - P^)^P(Xj). Similarly, 
the probability that signal is detected in the region X2 “ X^ 
during the third sweep, but not during the first two sweeps, is 
7j(l “ Pd) [F(L 2 ) - P(Xj)], since this region has only been 
searched once before the initiation of the third sweep. The 
region L^~ has not yet been searched, and thus the 
probability of detecting the signal in this region is simply 
Pj[P(X3) “ P(X2)]. P3 is the probability of detecting the 
signal in any one of these three nonoverlapping regions during 
the third sweep, and is thus equal to the sum of these three 
terms, which is precisely the result of (15). 


III. Optimum Symmetric Search Strategy 

Suppose we design a search algorithm with search lengths 
Xj, X2, . . . , X^ in the first N sweeps, and suppose the 
resulting probability of signal acquisition in these N sweeps is 
which is computed in Section II as a function of the a 
priori probability density of the signal location, the number of 
sweeps Ny and the search lengths. Denote the acquisition time 
that is required for the probability of signal acquisition to 
reach Qj^ by For example, Pq 9 is the time required to 
arrive at a of 0.9. The basis of the following optimization 
procedure is that, regardless the actual value of 2^, the 
acquisition time is proportional to the number of chips 
searched, which is proportional to the total search sweep 
length 


N 




64 


= P^ [P(L^) - PiL^y\ +F^(1 - P^) [P(L^) - P(L{)] 
+P^{l-Pj^P(L^) ■ (16) 



Therefore our problem becomes: Determine the optimum 
search lengths SO that Qj^ equals the desired 

acquisition probability and 


N 

= E h 

' i=I 

is minimized, thereby minimizing our acquisition time. Our 
method of solution is to use the LaGrange Multiplier metliod.' 
Let 


^ = -Prf E E (19) 

k=\ k^i 


where 




( 24 ) 


lr = 


L. 

l_ 

a 


(25) 


Substituting for (the normalized chip search numbers) 
back into the equation for allows us in principle to solve 
for V: 


ejT ■ E 0 - 


k=i 


where \ is the unknown LaGrange Multiplier. Up to this point 
the theory is quite general, the only requirement being that 
the a priori density be unimodal and symmetric and thatP(E;^,) 
be differentiable. Since this problem was motivated by the 
desire to improve acquisition time for the spread spectrum 
receivers and since a reasonable estimate for the a priori 
iocation of the signal is Gaussian, we shall illustrate the theory 
by assuming that the a priori density function is Gaussian. 
Now with this assumption we have 



X 


X 





(26) 


This equation appears nearly impossible to solve analy- 
tically; however it can be solved simply on a digital computer 
by trial and error, picking so that equals the desired 
probability. This must be done for all A^; of course in practice 
a few values of N will be sufficient. Before we discuss the 
results of the computer solutions let us consider the improve- 
ment over an unoptimized sweep. 


Differentiating, we obtain 



we can solve ( 21 ) to obtain 


IV. Uniform A Priori. Density Sweep Strategy 

The usual strategy for sweeping to obtain acquisition is to 
start at the end of tlie uncertainty region, where the range 
delay is minimum, and then retard the range in increments of 
typically one-half chip. By sweeping from the minimum delay 
to the maximum delay tlie chances of acquiring a multipath 
signal are diminished. If tlie probability of detection, given 
that the received code and the reference code are aligned, is 
given by and if the a priori probability density function is 
uniformly distributed, then the cumulative probability of 
acquisition is as shown in Fig-4. If, for example, a probability 
of 0.5 is chosen as the desired probability of acquisition, then 
the curve could be read off of the abscissa and the associated 
time, denoted by Tq 5 , would be the time it takes to acquire 
with a probability of 0.5. 



A measure of the improvement of the optimized scheme 
over the uniform sweep scheme can be measured as follows. 
(23) Denote as the time to acquire witli a probability of Q using 
the uniform sweep approach. Next, let 7^ denote the time to 
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acquire with the optimized sweep. Then the improvement which can be obtained numerically by the Newton’s iterative 
factor of the optimized sweep over the uniform sweep is given solution, computing successive iterations according to 
by 



(27) 




(32) 


where the derivative of/(x) can be evaluated from (31) to be 


The acquisition time is then Tq Tq = Tq. Clearly 
since unity is achieved with the uniform sweep strategy, and 
therefore the method never increases acquisition time. Our 
result is for Q = 0.5 as a representative case. 




i:=l 


in (l-Prf)l 
e 

Vn[e*-(*- 1) ln(l-/>^)J 


V. Numerical Results 

As discussed in Section III, to obtain* the optimum sweep 
lengths we have to first solve for X' from (26), from which the 
optimum sweep lengths are given, for / = 1, 2, . . . , iV, by 


(33) 

Before we illustrate the procedure by numerical example a 
discussion on the choice of iV is in order. First we note that in 
order for a given N to give finite optimum sweep lengths, N 
must satisfy 


L. 

I 


2vT ay/inY - J In (1 - P^) 


(28)' 




k=i 


(34) 


As evident from (28), the solution of X' must satisfy the 
condition: 


In^- ln(l— P^)>0 (29) 

in order for L- to be meaningful. Thus it is convenient to 
define the left-hand side of (29) to be of the form 

Inf -ln(l-P^) = e=^ (30) 

A 

so that instead of solving (26) for X' we can solve for the root 
of the following equation: 


X erfy/e^ - (k~ 1) In (1 - jP^) - (31) 


since the right-hand side of (34) is the probability of signal 
acquisition with infinite sweep lengths (i.e., for the case 
jc = ®°). Next, an upper bound on can also be given by noting 
the fact that the smallest value of x is so that N must 
satisfy, from (31), 


Q^^Pa E Ak- D In 

(35) 

Equations (34) and (35) give lower and upper bounds of N, 
for given values of and Qj^, within which optimum 
solutions of sweep lengths are feasible. The following examples 
illustrate the results obtained by this optimum search tech- 
nique as compared to the uniform search algorithm. 

Example 1 

Suppose = 0.25 and the desired = 0.5. The admis- 
sible number of sweeps N for this case is 3 and 4. The sweep 
lengths in terms of a (standard deviation of the signal location 
probability density) are shown in Table 1. Also shown in 
Table 1 is the ratio between the uniform search acquisition 
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time and the acquisition time using this optimum procedure. 
The search region in the uniform sweeps is assumed to be ±3 a. 
The four-sweep search gives an improvement ratio of 1.73 over 
the uniform sweep search algorithm, i.e., a reduction acquisi- 
tion time of 41%. 

Example 2 

Suppose = 0.5 and = 0.9. It is found that only TV' = 4 
01 N - S are admissible since solutions do not exist for all 
other N. Table 2 gives the sweep lengths and the corresponding 
improvement ratios over uniform searches. 

As illustrated by these two examples the improvement 
factor of tliis optimum search strategy over the uniform search 
varies according to Pd> and N, Nevertheless, the improve- 
ment factor is always > 1. 


VI. Conclusions 

We have presented a method that can be used to optimize 
(mmiinize) Ihe acquisition time for a PN-type spread spectrum 
system when the a priori probability density function is not 
uniform. Specifically we have calculated, for an assumed a 
priori Gaussian density function, that, when tlie 0.5 prob- 
ability acquisition'-time was used as a measure, the acquisition 
time was reduced by 40% for a cell detection probability of 

0.25 when three sweeps were used. For the same parameters 
and with four sweeps the acquisition time was reduced by 
41%. When the acquisition time probability was set to 0.9 
instead of 0.5 the reduction was 25% of the uniform sweep 
time. 

This technique Has application to the DSN if psuedo-noise 
spread spectrum systems are utilized to combat RFL 
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Table 1. Pd = 0.2S, Qn = 0.5 



Opt search with 

Opt. search with 


3 sweeps 

4 sweeps 

Ljo. 

2.54 

1.19 

1 

2.96 

1.93 


3.33 

2.46 



2.89 

Improvement 
ratio over uniform 

1.66 

1.73 

search 




Table 2. Pd = 0.5, Qfj = 0.9 


Opt. seajch with 
4 sweeps 

Opt. search with 
5 sweeps 

L fa 

2 63 

0.56 

i 

3.53 

2.42 


4.24 

3-38 


4.85 

4.12 



4.74 

Improvement 

1.34 

1.34 

ratio over uniform 



search 
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UNIBUS Monitor for POP 1l' 


M. D. Dormer 
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A UNIBUS Monitor has been designed and constructed to facilitate development of 
hardware interfaces with the PDF 11 minicomputer The Monitor provides useful displays 
of UNIBUS conditions and provides the user with a flexible diagnostic tool. It can also 
serve as a simple display and data entry device^ permitting extremely simple input! output 
for development software. 


i. Introduction 

A UNIBUS Monitor has been designed and constructed to 
facilitate development of hardware interfaces with the PDF 11 
minicomputer (Figs. 1 and 2). The Monitor provides useful 
displays of UNIBUS conditions and provides the user with a 
flexible diagnostic tool. It can also serve as a simple display 
and data entry device, permitting extremely simple input/ 
output (I/O) for development software. At this time the 
Monitor is being used with the DSN Planetary Radar System, 
which uses a PDF II. 


li. The Problem 

While the PDF ITs UNIBUS system simplifies peripheral 
and I/O interfacing and handling, it also introduces a single 
point of failure which can bring down a whole system. If a 
peripheral or I/O device fails in a way that interferes with 
UNIBUS functioning, it is possible that the processor will be 
unable to execute the code. In a situation like this it would be 


^ UNIBUS and PDF are registered trademarks of the Digital Equipment 
Corporation. 


convenient to be able to view the UNIBUS status directly 
without having to probe all 56 UNIBUS signal lines individ- 
ually. In a static situation this would be comparatively simple; 
simply hook an indicator light up to each signal line. Since the 
UNIBUS is a dynamic system involving various time-varying 
signals which typically change too fast for an engineer or 
technician to see^ some capability must be added in the 
monitor to recognize and display this dynamic activity in a 
useful way. 

Hi. The Solution 

A. Single Address, Synchronous 

In this mode (MONITOR switch - ‘‘OFFr SYNCHRO/IO 
switch- "SYNCHRO” SA/RANGE switch- “5^”) the 
UNIBUS Monitor watches for bus transactions involving the 
bus address specified by the top row of octal lever wheel 
switches. If Ae A lines contain that address the monitor 
latches the A lines on the arrival of MSYN. If the C lines 
indicate a DATO or DATOB operation, the Z) lines are latched 
at the same time. If the operation is zDATI or DATIP, the/) 
lines are latched on the arrival of SSYN. 
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B. Range, Synchronous 

This mode is similar to mode A (MONITOR — “OFF,” 
SYNCHRO jlO- ^^SYNCHROr SAjRANGE - ^^RANGE'^) 
except that instead of responding only to the address dialed 
into the A switches, the Monitor responds to all addresses less 
than or equal to the A switches and greater than or equal to 
the D switches. 

C. Single Address, 10 

In this mode (MONITOR - ^^OFFr SYNCHRO / lO - 
^^lO:^ SAjRANGE - “5y4”) the UNIBUS Monitor serves as a 
UNIBUS peripheral, responding as a peripheral or input/ 
output device would. It responds to the address contained in 
the upper row of lever wheel switches. When the A lines 
contain the correct address and MSYN is asserted by the bus 
master, the UNIBUS Monitor latches down the A lines and 
responds with SSYN, On DATO or DATOB operations the 
Monitor latches down the D lines at the same time as the A 
lines. On DATI and DATIP operations the Monitor asserts the 
contents of the lower row of lever switches onto tlie D lines 
until MSYN is released. In 10 mode the monitor observes the 
UN I BUS handshaking protocols according to Digital Equip- 
ment Corp. specifications. 

D. Range, 10 

This mode (MONITOR - ^‘OFFr SYNCHROjlO - ^^lOr 
SA jRANGE — “FA/VGF”) is essentially the same as Mode C 
except that the Monitor responds to the range of addresses 
between A and D inclusive. On DATI and DATIP instructions 
the contents of the lower row of switches are asserted onto the 
D lines as in Mode C. In this mode the lower address boundary 
and the returned data are the same, but this shouldn’t cause 
any problem. 


E. Monitor Mode 

In this mode both the D and A latches are held enabled. No 
dynamic responses are made. All addresses are displayed for as 
long as they are present. This mode is useful for taking a look 
at the UNIBUS in static situations, such as when the processor 
is stalled. 

F. Non-Modal Displays 

All the other UNIBUS signals are displayed on LEDs 
directly on the front panel with no dynamic interactions with 
the bus. In addition, the C lines are decoded into the four bus 
transactions, DATO, DATOB ^ DATI, DATIP, and displayed. 
Also, the MSYN, SSYN, SACK, and BUS GRANT lines have 
pulse stretching one shots so that the signal and an “activity” 
indicator are both displayed. In the case of the BUS GRANT 
lines the one shot is triggered by the assertion of any grant 
line; there is only one activity indicator for the grant lines. 
One particular reason for activity indicators on the GRANT 
lines and on SACK is that should a UNIBUS device fail with its 
BUS REQUEST asserted and be incapable of responding to the 
BUS GRANT, the processor will attempt over and over to 
grant the request without success, timing out and starting over 
again each time. The GRANT ACTIVITY indicator will catch 
this kind of failure. 


IV. Conclusion 

After construction and testing the UNIBUS Monitor was 
taken out to DSS 14, where it was installed with the planetary 
radar equipment. In this application the Monitor soon proved 
its usefulness in analyzing UNIBUS communication 
performance. 


71 


UNIIUS 


DATA 

LINES 


ADDRESS LINES 


H 


COMPARE 


n 


TOP/SINGLE 

ADDRESS 

SWITCHES 


U 




ADDRESS 

DISPLAY 


LATCH ENABLE 


COMPARE 


DRIVER 


y 



BOnOM DATA 
SWITCHES 




OUTPUT ENABLE 


DATA 

LINES 


) 


DATA 

DISPLAY 


LATCH ENABLE 


A= T AST A>B 


CONTROL 

LOGIC 


SSYN 


MSYN 


Cl 


CO 


i 


INTR 


INIT 

r 

BBSY 


c 

BUS 

^1 

REQUESTS " 

BUS 



POWER LO 

^ c 

PARITY 1 



LED 
DISPLAYS 


PULSE 

STRETCH 


PULSE 

STRETCH 


-O 

■O 


MONITOR 


•SA/RANGE 


MODE 

INPUT 


SYNCR/IO 


LED 

DISPLAYS 


DECODER 


(DATO 
» DATOfi 
I DATI 
) DATIP 


PULSE 

STRETCH 


PULSE 

STRETCH 


Fig. 1. UNIBUS Monitor block diagram 


UNIBUS MONITOR 



IfL BOTTOM AOORCSS/OATA OUT 


MODE 

STRCMI 

NU 




RARCC 

MONITOI 






0 — 88YM 

— 0-ArnviTY 

0 BRA 

0 

% 


0 MSYN 

# - ACTIVITY 

0 BRS 

#iCS 


0 SACK 

— 0 — Acmny 

0 BR6 

• Bce 

0 DATI 

0 INTR 

0 *MIT 

0 BRT 

#BC7 

0 OATir 
0 OATO 

0 BB8T 


0 NfR 


0 OATOB 



6RAXT 0 



ACnVfTY 


Fig. 2. UNIBUS Monitor panel 


72 




DSN Progress Report 42-46 


N78-32l4ft 





May and June 1978 


DSS 13 Antenna Subsystem Automation 

H. Phillips, I. Crane, and P. Lipslus 
DSN Engineering Section 


Unattended station operation has been implemented at DSS 13, permitting full 
operational control from NOCC, Sensors have been installed in the mechanical 
subsystem to monitor critical functions. These are arranged to permit automated 
pr emission checkout of the subsystem, automated reaction to component failure, and 
identification of failed components under control of the antenna pointing computer^ 
This monitoring installation mil serve as a prototype for monitoring equipment to be 
installed throughout the DSN. j 


I. Introduction 

Unattended station operation has been implemented at 
DSS 13, permitting full operational control from the Net- 
work Operations Control Center (NOCC). Sensors have been 
installed in the mechanical subsystem to monitor critical 
functions. These are arranged to permit automated premis- 
sion checkout of the subsystem, automated reaction to 
component failure, and identification of failed components 
under control of the antenna pointing computer. This moni- 
toring installation will serve as a prototype for monitoring 
equipment to be installed in the DSN. 

II. Automated Premission Checkout 

The function of a premission checkout is to assure that 
all of the antenna subsystem equipment required for tracking 
and for automated response to an equipment failure is 
functioning. The automated checkout, controlled by the 
Modcomp pointing computer, was planned to be carried out 
in the following steps: 


(1) Verify that the wind velocity on .site is within 
acceptable limits. 

(2) Verify that the level in the servo hydraulic fluid 
reservoir is within the prescribed limits. 

(3) Start the servo high-pressure pumps and the filter 
circulating pumps and verify that these are operating. 

(4) Verify hydraulic pressure at the azimuth and eleva- 
tion servo valves. 

(5) Start the gear drive lubrication pumps and verify 
lubrication pressure or flow. 

(6) Test the drive brakes as follows: 

^ (a) Set in a zero rate command for each axis. 

(b) Release the brakes for each axis, and verify tliat 
the brakes release. 

(c) Reset the brakes for each and verify that they 
set. 

(7) Check the drive operation as follow : 
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(a) Release the brakes and sound the horn. 

(b) In low-speed mode, command a small rate in each 
axis. 

(c) Calculate A0/Ar for each axis over a small period 
and verify that the antenna is moving at the 
commanded rate. 

(d) Read the analog rate signal from the antenna and 
verify that it is within tolerance compared to the 
calculated A0/Ar, 

(e) Command a zero rate and set the brakes. 

III. Monitoring During Operations 

Monitoring during operations is required to determine if 
the mechanical components are operating within normal 
parameters and to initiate, by computer control, defensive 
action to prevent a catastrophic failure. The sensing points 


were chosen to isolate major potential failure points, and to 
provide the computer with the information necessary to 
initiate a defensive reaction when called for. The functions 
monitored and the action to be taken on a sensed failure are 
listed in Table 1. The actual sensed points are listed in. 
Table 2; 

IV. Present Status 

All of the monitoring equipment has been installed on the 
antenna and is presently operational. The microcomputer for 
collecting data and transmitting operational data has been 
installed and is operational. A demonstration of remote 
operation controlled from NOCC has been conducted. Opera- 
tional experience is now being developed to validate the 
checkout. Monitor and failure response rationale data are 
also being developed on the normal operating range of the 
various parameters prior to implementing the final monitor- 
ing software in the antenna controller. 
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Table 2. Sensor implementation 


Table 1. Monitoring during antenna operation 


Parameter 

Action on failure 

Reservoir oU level change 
Slow leakage rate 
Large leakage rate 

Drive to stow 

Stop pumps immediately 

High-pressure pump and 
motor operation 

Motor current excessive 
Pressure output excessive 

Shut off system 
Shut off system 

Low pressure at servo valves 

Shut off high-pressure pump and 
filter circulation pumps 

Gear drive lubricant flow 
failure 

Drive to stow, set brakes 

System oil temperature 
excessive 

Drive to stow, set brakes 

Wmd velocity excessive 

Drive to stow, set brakes 

Cable wrap-up limits 

Command zero rate, set brakes 

Tracking error excessive 

Shifrto low rate command in 
high-speed mode and drive to 
stow 

DC power supplies failure 
(control room) 

Automatic servo system shut- 
down if both power supplies fail 

Remote DC power supply 

Automatic servo system 
shutdown 


In control room Clockwise cable wrapup limit 

Counterclockwise cable wrapup limit 
Azimuth 1 brake status 
Azimuth 2 brake status 
Elevation left brake status 
Elevation right brake status 
Azimuth high-speed niode select 
Elevation high-speed mode select 
Disable circuit status 
Power supply 1 status 
Power supply 2 status 
Azimuth prelimit status 
Elevation prehmit status 

At antenna Elevation low-speed servo valve pressure 

differential 

Elevation system pressure 
Elevation high-speed servo valve pressure 
differential 

Elevation high-speed system supply pressure 
Elevation low-speed system supply pressure 
Azimuth high-speed servo valve differential 
pressure 

Azimuth low-speed servo valve differential 
pressure 

Azimuth high-speed system supply pressure 
Azimuth low-speed system supply pressure 
Fluid level (hydraulic tank) 

Wind speed (southwest wind tower) 

Wind speed (southeast wind tower) 

Wind direction (southwest wind tower) 
Wind direction (southeast wind tower) 
Hydraulic oil temperature 
Lube oil pressure 
75-horsepower starter status 
125-horsepower starter status 
Pump flow 75 horsepower (left) 

Pump flow 75 horsepower (right) 

Pump flow 125 horsepower (left) 

Pump flow 125 horsepower (nglit) 
Hydromech building 28 V DC power supply 
status 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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FPLA Mechanization of Arithmetic Elements to Produce 

A + B or to Pass A Only 

D. E. Wallis, H. Taylor, and A. L. Rubin 
Communications Systems Research Section 


772/5 article describes a 4'hit and a 3-bit adder which can be implemented under 
special hardware restrictions. The chip to be used is Field-Programmable Logic Array 
(FPLA) with 12 input lines, 50 AND gates inside, and output through only 6 OR gates. 
The context in which it is being used requires an '^enable” function which can suppress 
one of the two numbers to be added. The 3-bit enabled adder is compatible with 
lookahead-carry mechanizations using the 74S182. It will be used in the accumulator for 
the RFI project. 


i. Discussion and Result 

The 4-bit adder is pictured in Fig. 1. The point of interest 
is the reduction in the number of AND gates accomplished by 
sharing six Boolean terms between and . A sequence of 
these 4-bit adders to add two 4/2 -bit numbers in ripple-carry 
fashion would be rather slow, because each next adder would 
have to wait for the previous cany. 

The 3 -bit adder is pictured in Fig. 2. Sixteen of these will 
add two 48-bit numbers faster than twelve of the 4-bit adders. 
The increase in speed is achieved by a well known “look 
ahead’’ arrangement. Each 3-bit adder puts out a P bit and a 
G bit instead of a carry bit. The G bit will have value 0 if and 
only if a carry will be generated. The P bit will have value 0 
if G has value 1 and a carry would be propagated. Using the 
Ps and Gs from all the 3-bit adders at once, the whole arrange- 
ment can get the sum of two 48-bit numbers in only 5 gate- 
delay times. 

The main point of interest in the 3-bit adder is the sharing 
of AND gates by way of the function F. This makes it possible 
to implement all the functions desired with the given chip. 


The other point of interest is that F does not depend on Cq. 
This means that F will find its final value during the time that 
Cq is being computed in the lookahead network. Thus when 
Cq arrives all the gates affected by Cq will find their values in 
just 1 gate-delay time. 

The practical problem was to mechanize a 48-bit digital 
integrator (or accumulator) whose memory could be read and 
cleared at the same time that a new datum d^ was added to 
begin a new integral. The memory type being used was not, 
itself, clearable. The form of the integrator is as shown in 
Fig. 3. 

The available 74S381 4-bit arithmetic/logical unit (ALU) 
almost met the requirements. The 74S381 can be controlled 
so as either to add its two inputs (A and B), or to transmit 
zero output. The zero output clears the memory, as desired, 
but the time-coincident incoming datum is lost and does 
not initialize the next summation. 

^Vhat was required, then, was an ALU which could add two 
numbers A and B, or transmit one of them (called “pass A”). 
No such ALU was available commercially. It seemed, however. 
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that an acceptable ALU could be created by programming a 
field-programmable logic array (FPLA) such as the 74S330, 
whose package size and pinout were compatible with the 
74S381. 

The 74S330 has 50 AND gates, each of which may be pro- 
grammed to AND any of 12 input variables using arbitrary 
assignments of which “rail” 'or true/false sense of each input 
variable would be effective on each AND gate. Further, any 
combination of the AND gate outputs can be ORed into any 
of six internal OR gates. Each OR gate is made available as an 
output from the 74S330, and each OR output can be pro- ( 
grammed so as to invert its output, if desired. 

The 48-bit adder logic was to be built on a circuit board , 
having little room (“real estate”) for extra components. It 
was, therefore,- necessary to maximize the number of bits 
which could be added by each 74S330. Further, the adder 
had to be a “fast” adder, i.e., compatible with carry lookahead 
techniques, such as those based on the 74S182 lookahead 
generator. It soon became apparent that the availability of 


only 50 AND gates in the FPLA was going to be a major 
limitation in the adder word-size or adder speed or both. Thus 
a combinatorial study was made in an attempt to obtain an 
optimal design for an adder based on the 50-gate FPLA. 

First choice would have been a 4-bit adder with the P, G, 
and E functions. We think that is impossible although we do 
not know how to prove it. 

The design chosen after several trials was the 3-bit adder of 
Fig. 4. 

Arthur L. Rubin has obtained the following result. Subject 
to the requirements of the chip 74S330, a 3 -bit adder with 
P, G, and E which makes no use of a tie around (such as the 
connection of the output F to other FPLA inputs, as seen in 
Fig. 2) from the output of the chip (and hence avoids the one 
extra gate delay) cannot be designed with fewer than 50 AND 
gates. Achieving the 3-bit adder of Fig. 4 on one chip makes 
the 48-bit adder as fast as possible. 
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Fig. 4. 3-b1t, fastest lookahead-compatible adder, with enable, using no intermediate outputs 
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A Probabilistic Version of Sperner’s Theorem, With 
Applications to the Problem of Retrieving 
Information From a Data Base 

L D. Baumert, R, J. McEliece, E, R. Rodemich, and H. Rumsey, Jr. 

Communications Systems Research Section 


We show; how the design of an optimal **merged key code** information retrieval 
system involves finding the probability distribution on n-bit binary words that minimizes 
P {X <Y 2 ^ .. . where X, Yj, , , , , Y^are selected independently according to the 
given probability distribution. We then find the minimizing probability distribution in the 
case r- 1. 


I. Introduction 

In retrieving information from a large data base, such as will 
exist in the DSN’s digital RFI Surveillance System (Ref. 1), 
the technique of merged keycodes (Ref. 2) is often useful. In 
this technique each record in the data base is assumed to have 
a certain number of attributes. Each possible attribute A is 
assigned a binary code (whose length is normally that of one 
computer word) C (A), If the record R possesses attributes , 

it is then assigned the merged keycode C{R) = 
C(^i) u C(.^ 2 ) U . . . , where the symbol u denotes the 
logical OR operation. 

If one wants to locate all records in the data base possessing 
a fixed set of attributes, say ^ 2 , . . . , one computes 
the merged keycode Z) = C (Rj^) U . . U C (R^), and then tags 
each record R such that D < C(R). Clearly the set of records 
with the desired attributes is a subset of the tagged records. 
However, some of the tagged records wiU not have the desired 
attributes; such records are called false drops. It is obviously 
desirable to minimize the numbei of false drops, other things 
being equal. 


'Under certain circumstances, it is reasonable to model the 
above situation as follows, A fixed number, say r, binary 
codewords are selected independently according to a certain 
probability distribution P, which is to some extent controlled 
by the system design. Denote these codewords by Y^, Y 2 , 
. . . , — they represent the keycodes of the attributes of a 

randomly selected record. Let X^, . . , , X^,s^r denote the 
keycodes corresponding to the attributes in a random query of 
the data base. We assume the X^ are chosen independently 
according to the same probability distribution F, The proba- 
bility of false drop is then 

p u . . . u u . . . u y^} = p {Xj < y^v ... uY^y 

Thus we are led to ask: What is the probability distribution on 
«-bit binary words that minimizes the probability 
P{Xj < Kjl U . . . U for various values of rl 

' In this article, we will solve this question for the case r = 1. 
It is hoped that the techniques developed can be brought to 
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bear for larger values of r. In the remainder of this section, we then clearly < 7} = = Y} - 1/M. Hence by 

will describe our main result. Theorem 1, 


Let be a finite set with n elements, and let V = K(^2) 
denote the collection of all subsets of SI. For xeV, Ixlwill 
denote the cardinality of .xr. The relations of set inclusion^and* 
proper set inclusion will be denoted by x< y and x< y^ 
respectively. The empty set, viewed as an element of V, will be 
denoted by 0, and itself, as an element of F, by 1. 

Let p(x) be a probability distribution on F, and let X and Y 
be elements of F chosen randomly and independently accord- 
ing to p. We denote by jP{A'< F} the probability that X will 
be a subset of F. Our main result is the following. 


Theorem 1. If n > 1, then for any probability distribution, 


P{X<Y}> 



1 


Furthermore equahty holds for the probability distribution 
defined by 



This is Sperner’s famous bound (see Ref. 3, for 
example) on the maximum number of non-comparable 
elements in F. Thus (for n > 1) Theorem 1 can be 
viewed as a generalization of Sperner’s bound. , 

(3) For even values of n, it is possible to show that the 
particular probability distribution cited in the state- 
ment of Theorem 1 is the only one for which equality 
holds. For odd n > 3,.the obvious alternate distribution 


/- 

n 



p(x) = 1 

2^ 

1 if\x\ = 

n 

1 _ 


= 0 otherwise 



= 0 otherwise 


In the remainder of this section, we make some remarks 
about Theorem 1 and its proof. The proof itself occupies 
Sections TT—V. 

(1) The restriction « > 1 is necessary, since with « = 1, 
the probability distribution p(0) = p{l) = 1/2 gives 
F{X < F} = 3/4, whereas 



(2) Let F = {yi , . . . , be a set of M pairwise non- 
comparable elements of F, i.e., y^ <y j iff i = /. If we 
define a probability distribution p(x) by 

i(xcY 
= 0 otherwise 


also achieves equality. For odd n> 5 this is the only 
other equahty-achieving distribution. However, for n = 
3 there are infinitely many extremal distributions: 

p(x) = a if lx 1 *= 1 

= 113 - a if |x I - 2 

= 0 if lx! = 0or3,for0<a<l/3 

(These facts can all be proved by making a careful 
study of our proof of Theorem 1. For the sake of 
brevity, however, we shall omit the details.) 

(4) Our proof contains two main ideas. First, by using 
more-or-Iess standard calculus techniques, one can 
obtain very strong necessary conditions satisfied by any 
extremal probability distribution. This we do in 
Section II. Second, we derive in Section IV (after some 
preliminary material in Section III) a lemma dealing 
with the expected number of maximal chains through a 
point of F which is selected according to a given 
probability distribution. This lemma was motivated in 
part by Lub ell's classic proof (Ref. 3) of Sperner's 
theorem. Finally in Section V, these two ideas are 
combined to give our proof of Theorem 1. 
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(5) *A possible alternate form of Theorem 1 would concern 
minimizing the probability that X and Y are compara- 
ble. It turns out that this problem is much easier to 
handle; the result is 


F{X<: Y 


or 


X> Y}> 



-1 


for all tiy and equality occurs only for a uniform 
probability distribution on the subsets of cardinality 


, Let 5 = min :x e G}, and let U denote the Eu- 

clidean neighborhood of the function p (restricted to G) 
defined by 


. (.qQc) - pix))^ <^1 ( 6 ) 

V xeG ) 

Clearly if q e U, then qQc) > 0 for all :x: e G. Thus by the 
assumed extremal property of p, we have 

(2g(<7) > Qq(p), foiaHq eU satisfying (5) (7) 



However, this inequality follows already from a 
theorem of Motzkin and Strauss (Ref. 4), together with 
Sperner’s original theorem. 


Thus p (restricted to G) gives a local minimum of the func- 
tion Qq, subject to the constraint (5), and so by the 
Lagrange multiplier rule there exists a constant X such that 



for X e G 


(S) 


II. Lagrange Multipliers 

The probability Q (p) = P {X < 7} is given by the sum 

Q(p) = 2 {p(x)p 0) : X < 7 } (1) 

We are asked to minimize the function Q of the 2" variables 
(p(:):), X eV), subject to the following constraints; 


Using (1), this becomes 

2pQc) + P(y) + S P(y) - X, for all :c e G (9) 

y<x y>x 

Furthermore, since F{X < Y} = F{X> l^}by symmetry, 


xev 


p{x) > 0, all X e F 


(•2) 2P{X < n = F{Z <Y}+ i>{2f > y} 


Suppose p is a probability distribution that minimizes Q (p) 
subject to (2) and (3), and let G = {x : p(x) > 0}. Temporarily 
we regard p and G as being fixed. Consider now the new 
problem of minimizing the function. 

= 2 {^(x)<3f 0) ; < y, X, 7 e G } (4) 

where ^ is a real-valued function defined on G, subject to the 
single hneaf constraint 


= 53 P(x) hpix) + Yt P(y) + Z} p0)1 

xeG \ y<x y>x ) 

= X p(x) = X, by (9) 

xeG 

Thus we 'have identified the constant X in (9), and so we 
have 


^ = 1 
XeG 


(5) Mx) + Z p(v) + E p(y) = < y>, 

y<x y>x 


(N.B, q is not required to be a probability distribution.) 


forallxeG (10) 



Equation (10) is the condition which must be satisfied by 
any extremal probability distribution that we will return to 
in Section V. 

Note: Equation (9) follows immediately from the Kuhn-Tucker 
theorem of nonlinear programming (Ref. 5), and indeed one 
can also conclude from K-T that the left side of (9) is > X if 
X We have included this elementary derivation only to 
make our exposition more self-contained. 


til. Preliminaries About Chains 

A chain of length r in F is an (r + l)-tuple c- (xp, 

. . . , of elements from V such that Xq <x^ < , , . <x^. 
If in addition we have [a:, -+ 11= 1^/1+ 1, c is said to be a 
maximal chain (of length r) from Xq to Such a chain is 
said to pass through each of the points . ,x^. 

If • 5 is denote by 

MCiy^, . . . , y^) the set of all maximal chains from to 
y^ which pass through each of the The number of 
maximal chains m MCiy^^ . . . , y^^ is denoted by /Oq? 
.. . , 7 ^). Thus 

(11) 

If (yo> • • • > ‘ j ym) ^ chain, it is clear that every 

maximal chain in MC(y^, can be decomposed 

uniquely into a chain from MCOq, . . . , y^ followed by a 
chain horn MCXy icy •••>>'/«)• Hence 

/Oo= • • ■ , • • • > j'm) = /Oo’ • • • >yk)f^k> • • ■ > J'm) 

( 12 ) 

and by mductiom it follows that 

m-l 

f(yo>--->ym)= n f(y,,yi+i) (i3) 


In view of (1 3), in order to compute / (c) for a general 
chain c, it suffices to consider the case where c= (x, y) 
consists of only two elements. This we now do. 

Let X - Xq ^ ^ ^ Xj. - y he 2 i maximal chain 

from X to y and let - x^ = {wj, i= 0,1, . . . , r ^ 1. 
Then (wq, Wj, . . . , w^_i) is a permutation of the elements 


in the set y - x. Conversely, if (wq, . . . , is any per- 
mutation of y - Xy and if we define x^= x ^ {wq, , 

}, then (Xq, . . . , Xy) will be a maximal chain from = 
X to Xy= y. Hence there is a one-to-one correspondence 
between chains from x to j and permutations of - x: 

/(pc.y')^^ \y- x{\ (14) 

m-i 

f(yo> • • • . :>'m) = n - j/l ! (15) 

£=0 

As a final bit of notation, for x e 7 let N(x) be the 
number of maximal chains from 0 to 1 passing through x. 
Then 


N(x)=f(0, X, 1) 

= /(0,x)/(x, 1) 

= 1x1-! («- Ul) !=«!/( ”j) (16) 


Note that as a function of x, N(x) achieves its minimum 
value when [x| -nl2 for n even, and [x| = (« ± l)/2 for n 
odd. Thus if we define o:(«) = [nl2] I (n - [n/2]) !, we have 

» A^(x) > a(n\ all X e 7, (17) 

with equality iffx = [n/2] oi n - \n/2 ] . 


IV. A Basic Lemma 

Let G be a subset of 7. If = (Xq, . . . , x^) is a chain in 
7, and at least one element of c lies m G, we define Lq(c) 
(“the last element of c lying in G”), as follows: 

L(n(c) = x^, where k = max {/: e G} (18) 

If no element of c lies in G, Lq(c) is undefined. 

Further, we define for each x e 7, 

' A^^(x)=|{ceMC(0,l):L^(c)=x}| , (19) 
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Thus Nq(x) is the total number of maximal chains from 0 
to 1 whose last element in G is x. If now for each jc e we 
define g (x) by 


s(x) = \{c e MC (x, 1): Lq(c) = a:}1 (20) 


i.e., ^(x) is the number of maximal chains from x to 1 
whose last element in G is itself, it follows that 


No(x)-mx)g(x) (21) 

If p is a probability distribution on K, and if X is an 
element of V chosen randomly according to p, the number 
of maximal chains from 0 to 1 passing through X, N(X) is a 
random variable whose expectation is given by E{N{X)) = 
]S1p(x)jV(x): xeV}. The following lemma gives another 
formula for E(N(Xy) which is crucial in our proof of 
Theorem 1. 


Replacing /(x, 1) in (22) by the sum (23), and interchanging 
the order of summation, we get 


p(x)N(x) = ^0) ^ pQc)f(p,x)f{x,y) 

xeV. yeG x<y 

= IJ my)g(y) I; p(x) — 

yeG x^y •' 


But by (21), /(O, y)g(y) = N^fy). And by (14), 

f(0,x-)f(x,y) \x\l\y-x\l /|y|\-^ 
f(0,y) bl! ~\\x\l 


This proves Lemma 1. 


Lemma 1 


Let p be a probability distribution on V, and let G = 
{x € V: p(x) > 0}. Then 


E 

xeV 


p(x)N(x) = Y, ^-G^y) 

yeG 



Proof 

Using the fact that p(x) = 0 if x ^ G, and (16), we have 


V, Proof of Theorem 1 

We are now in a position to give a short proof of 
Theorem L The idea of the proof is to estimate the expected 
number of maximal chains through a randomly selected point 
of V in two ways. On one hand, this expectation is certainly at 
least a{n) by (17). On the other hand, using the machinery we 
have developed m Sections 11“ IV, we will show that this expec- 
tation (at least for an extremal probability distribution) is at 
most Y'}‘ n\ The resulting bound, E}>ok(«)/ 

n!, is the bound of Theorem 1 . Let us now see how tliis proof 
goes in detail. 


Y E fWr(o. x)f{x, 1) (22) 

xev xeG - 


We now^classify the chains in MC(x^ 1) according to their 
last element in G. If y>x, then the number of maximal 
chains from x to 1 whose last element in G is y is 
f(x,y)g(y). Thus 


Let p be a probability distribution that minimizes Q(p) = 
Yl and let G= {xeV: p{x) > 0}. Then by Lemma 1, 


Y P(.x)N(x ) = Y E 1^1 ) ' (24) 

xer yeG x^y 


% (1^7), N(pc) > aQi) for all x. Hence 


f(x, l)=Y y^Siy) (23) 

y'^x 


^ p{x)N(x) > a(«) (25) 

xeV 
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Let us rewrite the inner sum in (24) in the following way: 



= p(y) 

x<y 


cr' 


(26) 


<p(y) + S p(-^) ^ 

^ x<y ^ x>y 


( 28 ) 


We now claim that if n > 1, and ^ < 7 , 


We now apply (10) to (28) and conclude that 


P(jc) 



< JP(^) 


(27) 


Ys P(^) (1^1) <P{X< Y},' ify e G (29) 

^y 


If X 0, (27) is obvious since then the binomial coefficient 
will be > 2. If X = 0, the binomial coefficient is 1, but (27) is 
true anyway because p(0) = 0. This can be seen as follows. 
Let ;c = 0 in (10); we get F{X< 7}= 1/2 + 1/2 ^(O). Hence 
if p(0) > 0, then F{X < 7} > 1/2. But since « > 1 we can 
select two non-comparable elements and and define a 
probability .distribution q by setting q(x^) = q(x 2 ) = 1/2. For 
this probability distribution we clearly have Y}= 1/2. 

This shows that no distribution that minimizes F{X<Y} for 
n > I can have p(0) >0, and this completes the proof of 
(27). 


. Combining (24), (25), and (29), we get 


a(«) < P{X < 7} X) ^c(y) (30) 

yeG 

Finally we observe that each maximal chain from 0 to 1 
is counted at most once in the sum (30). (If c is such a 
chain and if Lq(c)- it is counted by the term A^^(y).) 
Since the total number of such chains is n\ (see Eq. 14), we 
get 27V^(y) < n\ and hence, finally, 


Combining (26) and' (27), we get 


Y ( !^I ) p^^ i 2 ^.(^) 

x<y 


P{X< Y}> 


a(n) _ 



(31) 


This completes the proof of Theorem 1. 
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VLBI-Laser Intercomparison Project 

B. D. L. Mulhall 
TDA Engineering Office 


The VLBI-Laser Intercomparison Project was established at the direction of NASA to 
assess state-of-the-art geodetic measurement systems being developed by NASA, A Project 
plan describing the objectives of the Pvject, the methods for making the assessment, and 
the schedule, was reviewed. The plan was approved and published. This article describes 
the contents of the plan. 


i. Introduction 

A. General 

The primary objective of the DSN VLBI-Laser Intercompar- 
ison Project is the accurate measurement of baseline vectors, 
both length and direction, between established geodetic bench- 
marks.^ 

A second objective is the intercomparison between both 
satellite and lunar laser ranging techniques with VLBL Both 
geodetic and geodynamic measurement instruments are to be 
evaluated. This is to provide potential users with an assessment 
of the value of the different instruments. This assessment is to 
be performed in terms of accuracy of the baselines measured 
and/or the geodynamic parameters measured, the operability 
of each instrument and, finally, the cost of operations of each 
system. This leads to the demonstration of the suitability of 
the various systems to potential applications. 


surveyor’s mark made on a permanent landmark that has a known 
position altitude that can be used by various geodetic measurement 
systems. 


B. Specific Objectives 

The specific objectives of the Project are divided into im- 
mediate and future objectives. The immediate objectives are 
to: assess current VLBI and Laser System performance, 
identify potential problems 111 tlie application of VLBI or 
Laser Systems by NASA or other agencies, assist system devel- 
opment to overcome problems, and finally, demonstrate readi- 
ness for technology transfer. All these objectives are intended 
to bC' accomplished by 1979, except for the intercompanson 
between VLBI and Lunar Laser Ranging (LURE) which is not 
planned to be accomplished until early 1980. 

Beyond 1980 the objectives of the Project are to demon- 
strate performance of the five-centimeter VLBI system in an 
operational environment and intercompanson with LURE. 

C. Assessment Criteria 

'The VLBI and Laser Systems will be judged on tlieir ability^ 
to measure vectors between benchmarks in terms of relative 
accuracy, since no absolute scale can be found. By relative 
accuracy is meant the repeatability of each system’s measure- 
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ments over the measurement session and, in cases where appro- 
pnate, the ability to close figures such as triangles. 

Secondly, the operability, that is the ability to operate each 
system, will be assessed. This will be done by noting the 
differences between planned observations and actually' accom- 
plished observations of the measurement’s sources during each 
session. 

Finally, the life-cycle cost of each system will also be a 
basis for judgement. By life-cycle cost is meant the implemen- 
tation cost of an operational system and the opeiating cost 
over the life-cycle of the system. Operational costs will be 
measured in terms of man-months of effort per month, and 
shipping and setup costs rather than actual salary dollars to 
avoid evaluation dependence on salary structure of the user 
agency or contractor. 

D. Accuracy Demonstration 

The objective of demonstrating geodetic and geodynamic 
measuiement accuiacy can be satisfied by an experiment 
which will produce necessary and sufficient results. For exam- 
ple, if a more accurate measurement technique existed, then 
the accuracy of VLBI could be demonstrated by comparison 
to this more accurate technique. For short baselines where 
conventional surveys can provide centimeter accuracy, this 
comparably accurate technique does exist. However, on long 
baselines there is no demonstrated technique for producing 
comparable accuracies. 

Consequently, in the absence of a more accurate measure- 
ment technique, no known test of sufficiency has been found. 
What can be done is to satisfy an exhaustive set of necessary 
conditions to infer the accuracy. 


II. Error Source Evaluation Methods 

The VLBI error sources can be grouped foi analysis as 
follows. First, there are the VLBI subsystems, some of which 
can be isolated and examined singly. Others must be lumped 
together for evaluation. Table 1 lists the subsystems, the types 
of errors that arise, the methods for evaluating each subsystem 
or group of subsystems, and the measurement session in which 
the evaluation will occur (see paragraph III for measurement 
session schedule). 

A second set of error sources can be called natural. These 
are shown in Table 2 and include the structure and position of 
the radio source, the ionosphere, and the troposphere. In two 
of these instances, ionosphere and troposphere, related re- 
search has been conducted for some years by DSN Advanced 
Systems. Also, more work is planned in the future to evaluate 


the water vapor radiometer in an absolute sense. The Project 
will closely momtor these activities. But, they, are not part of 
this Project plan since they have been or will be carried out 
regardless of this Project’s sponsorship. 

It is the intent of the Project to conduct VLBI experiments 
in certain environmental conditions such as variable humidity 
to investigate the effect on the interferometry of water vapor, 
for example. 

System evaluations will be conducted as follows: First, 
short baseline closure and repeatability will be performed to 
investigate instrument errors. Also, comparisons with con- 
ventional surveys, such as NGS will be conducted. 

Secondly, contential baseline closure and repeatability can 
demonstrate VLBI capability and comparisons with satellite 
laser ranging can be made. 

Third, intercontential baseline comparisons with LURE are 
planned as intercomparison with the VLBI system. The 
Universal Time 1 and polar motion determined by both sys- 
tems will also be compared to establish the geo dynamic capa- 
bility of the two systems. 

III. Measurement Session Schedule 

The Project is organized and conducted as a series of 
measurement sessions. Each session is designed to achieve 
specific objectives as developed in later portions of this plan. 
Sessions are generally conducted once per year. Sufficient time 
should be allowed between measurement sessions to evaluate 
the results so that the planning for the next session can take 
full advantage of the knowledge obtained from the last session. 

The measurement sessions evolve from the simple to com- 
plex by increasing system sensitivity and performance, by 
adding more calibrations and by lengthening the baseline being 
measured. Also, there are sessions which are intended to pene- 
trate deeply into one aspect or subsystem of the system which 
may also be conducted and which may in fact be simpler 
measurements but performed under more highly controlled 
conditions to isolate error sources of the system. 

Figure 1 shows the overall VLBI-Laser Intercomparison 
Project schedule. The top two lines indicate the periods for 
data acquisition for each measurement session and the sched- 
ule for project leports. As shown, each session has a pre- 
session plan describing the objectives of the session and a post- 
session report listing the results. 

The planning for each session involves preparation of a 
pre-session report to document specific objective, measure- 
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ment strategy and technique, and estimate its estimated ac- 
curacies. The >pre-session report also identifies the hardware 
and software required and the schedule for data acquisition, 
reduction, and analysis. The measurement is then performed. 

The post-session report is published to document the re- 
sults. In the post-session report the experiences learned in 
terms of operational difficulty and the equipment require- 
ments for the next session should be identified. Most impor- 
tant of all, the accuracy actually obtained will be published. 

Also, there is an overall Project plan, which is this docu- 
ment, and there will be a. final Project report at the end of the 
Project. 

The next line item on the schedule arc flight project events 
and are included to show periods of peak activity and heavy 
loading on DSN ‘ facilities. As can be seen from the data 
acquisition measurement section schedule, sessions are for the 
most part intended to avoid these hea \7 deep space probe 
activity periods to facilitate antenna scheduling. 

The next three line items show the deployment schedule 
for the mobile satellite laser ranging equipment (MOBLAS), 
the Geoceivers, which employ satellite Doppler measurements 
to determine their position, and the ARIES antenna which 
provides a mobile arm for the VLBI system. The planned 
deployments of these three systems require coordination and 
the agreement of the respective agencies before the plan can be 
finalized. Specifically, the deployment of the Geoceivers of 
the National Geodetic Survey in Session Four at Haystack, 
Goldstone, and Ft. Davis, needs to be agreed upon. 

The next six line items, VLBI data acquisition racks, RF 
phase calibrators, water vapor radiometers, hydrogen masers, 
noise adding radiometers, and the Caltech- JPL MarkTI Cor- 
relator are major- equipment items that are required for certain 
sessions as indicated by the cross-hatched bars. Note that the 
Mark II Correlator will not be employed in Session Six since 


the DSN VLBI processor subsystem will be used for correla- 
tion and estimation for this final session. 

The next items are surveys being conducted by the NGS as 
part of certain intercomparisons mainly, the survey at DSS 14 
between the intersection of axes of this antenna, a benchmark 
or monument beside DSS 14, to the MOBLAS location pad 
and the ARIES pad 300 meters from the DSS 14 antenna. A 
second survey which is planned if the intra-complex report 
indicates that the anticipated accuracies are sufficient to pro- 
vide a highly precise intercomparison with VLBI will be con- 
ducted over the entire Goldstone complex. This survey will 
primarily be used to compare the DSS 14 to DSS 13 baseline. 
However, since survey techniques require the establishment of 
a network of measuring devices the other Goldstone antennas 
can also be surveyed without a great additional cost, it should 
be noted that DSS 13 will become a wideband VLBI station 
before DSS 14 and consequently DSS 13 will be used m cer- 
tain long baseline tests in place of DSS 14 where 14 is unavail- 
able due to other comments. Consequently, by measuring the 
13^14 baseline accurately, these DSSs can be used more or 
less interchangeably. 

The next item is the implementation schedule of the DSN 
VLBI System which will be employed in Session Six for the 
operational demonstration. DSN VLBI will produce UTl -polar 
motion routinely from July 1979. The precise DSN VLBI 
System (5 cm) will be operational in time for Session Six. The 
DSS 13 VLBI configuration is the next line item and is 
planned to be completed by October 1978. At that time 100 
MHz of spanned bandwidth at X-band would be available. 

Finally, the four participating facilities not within the DSN, 
namely, the Haystack Observatory, the Owens Valley Radio 
Observatory, the Ft. Davis Observatory, and the LURE equip- 
ment at the McDonald and the Haleakala Observatories, are 
listed as are their required participation periods as shown by 
cross-hatched bars. 


89 



Table 1. Subsystem error source evaluation 


Error sources 

Evaluation 

method 

Measurement 

session 

Instrumentation 
Frequency standard^ 



Phase stability 

Phase calibrator^ \ 

Phase center 
variations to 
intersection of j 

^ Short baselines 


axes j 

and connected 


Phase path variations i 

^ element 

D 

due to antenna dis- [ 

interferometry 


tortion or sub- ) 

reflector motion 

Off-axis pointing 
phase effect j 

Mechanical 

Variation between 

1 

Comparison 

2 

intersection of axes and 

^th^urvey 


benchmark 


- 

Software 

Comparison of 

2 


independently 
formulated 
code, compari- 
son of results 
wth identical 
input data 



^Evaluation provided by Advanced Systems. 


Table 2. Natural error source evaluation 


Error source 

Evaluation method 

Measurement 

session 

Source positions and 
structure 



Uncertainty 
source positions 

Baseline repeatability 

2&4 

Changing source 
structure 

Baseline closure 

2&4 

Ionosphere^ 

S- YS X-band solution 
(<1 cm) 


Troposphere 



Dry components® 

Surface measurements 
(<1 cm) 


Wet component® 

WYR (water vapor 
radiometer) calibra- 
tion (<2 cm) 



Measurements in 
highly variable 
humidity 



^Evaluation provided by Advanced Systems and/or ARIES 
Project. 
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A Fast Computation of Complex Convolution 
Using A Hybrid Transform^ 
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In this article, it is shown that the cyclic convolution of complex values can be per- 
formed by a hybrid transform. This transform is a combination of a Winograd transform 
and a fast complex integer transform developed previously by the authors. This new 
hybrid algorithm requires fewer multiplications than any previously known algorithm. 


I. Introduction 

Several authors (Refs. 1-9) have shown that transforms 
over finite fields or rings can be used to compute circular con- 
volutions without round-off error. Recently, Agarwal and 
Cooley (Ref. 10) used the techniques of Winograd (Refs. 11, 
12) to compute cyclic convolutions. These new algorithms for 
convolutions of a few thousand points require substantially 
fewer multiplications than the conventional FFT algorithm 
(Ref. 13). 

Previously the authors (Ref. 5) defined a class of Fourier- 
like transforms over the complex integers modulo This was 
a transform over the Galors field GF{q^), where q - 2^-1 is a 
Mersenne prime for p - 2, 3, 5, 7, 13, 17, 19,31, 61, . . .Re- 
cently these complex integer transforms were specialized to a 
transform length of d points, where d\%p (Ref, 8). The advan- 
tage of the latter transform over others is that it can be accom- 

^This work was supported in part by the U.S. Air Force Office of 
Scientific ResearclTunder Grant AFOSR-75-2798. 


plished completely by circular shifts, i.e., no multiplications 
are needed (Ref. 8). 

In this article, it is shown that Winograd ’s algorithm (Ref. 
11) can be combined with the above-mentioned complex inte- 
ger transform over GF{q'^) to yield a new algorithm for com- 
puting the discrete cyclic convolution of complex number 
points. By this means a fast method for accurately computing 
the cyclic convolution of a sequence of complex numbers for 
long convolution lengths can be obtained. This hybrid algo- 
rithm is comparable in speed to that given by Agarwal and 
Cooley (Ref. 10) and is implemented readily on a digital com- 
puter. The dynamic range requirements for this hybrid algo- 
rithm are presented here m detail. 

II. Cyclic Convolution 

The following algorithm for the cyclic convolution of two 
sequences is based on ideas given by Winograd (Ref. 11). Let 
the field of rationals be R. Also let X{u) =x^ -\-x^u -hx^u^ + 
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• • • + Y(u) + . . . +y^urt-i be 

polynomials over R. The product T(u) = X(u) • Y(u) can be 
computed by 

2n-2 

T(u) = X(u) • Y(u) mod J][ (w - (1) 

where q:^ e It is shown in Ref. 1 1 that a minimum of 2^z - 1 
multiplications are needed to compute (1). 

It is readily shown that the cyclic convolution ofX( 2 /) and 
Y(u) IS the set of coefficients of the polynommal 

r(u) = X(u) • Y{u) mod (w” - 1) 

Let the polynomial w” - 1 be factored into irreducible rela- 
tively prime factors, i.e., 

k 

Z=1 

where ” 1 for z Then T(u) mod^^(w) for i - 

1,2,...,^ can be computed using Eq. ( 1 ). Finally, the Chi- 
nese remainder theorem is used to evaluate T(u) from these 
residues. The above summarizes Winograd's method for per- 
forming a cyclic convolution. 


The following theorem is due to Wmograd (Ref. 11). 

Theorem 1: Let a and b be relatively prime positive integers 
and A be the cychc ab X ab matrix, given by 

+J^ mod a • b), 0<x,y<ab 

If 7T is a permutation of the set of integers {0, L . . . , -1} , 

let 

B(x, y) =A(nix), n(y)) 

Then there exists a permutation ti such that, if B is parti- 
tioned into b X b submatrices, each submatrix is cyclic and 
the submatrices form an zz X a cyclic matrix. 

It was shown previously (Refs. 10, 12) that the number of 
multiplications needed to perform a circular convolution of 3 , 
5, 7, and 9 points of complex numbers is 4, 10, 19, and 22 
multiplications, respectively. In order to compute the cyclic 


convolution of two longer sequences of complex integers, a 
£?-point transform over GF(q'^\ where q - - 1 andcflSj^, 

will be utilized here. Since the latter transform can be evalu- 
ated without multiplications (Ref. 8) , it can be used with con- 
siderable advantage to compute a cyclic convolution of two 
<i-point complex number sequences. The number of complex 
integer multiplications required to perform this circular convo- 
lution over GF(cp’') is precisely cf, the number of multiplica- 
tions needed to multiply together the transforms of the two 
sequences. 

For the moment, let d, the transform length, be an arbi- 
tary integer. Next let « P^ be the factorization of 

d into prune integers. If one lets * ^2 • * * 

= p^, then by Theorem 1 , a X <7 cyclic matrix can be par- 
titioned into b\ = P'1 matrices of size let = 

^2 X where a^-p^. . p^_^ and b’^ = p^„i . If is not a* 
prime, then each cyclic matrix can be partitioned into 

^2 matrices of size In general, X b^+^, where 

is a prime. If ^ 1 , then each X cyclic matrix can 
be partitioned into bf^^ matrices of size X Other- 
wise, the procedure terminates. If the number of multiplica- 
tions used to compute the cyclic convolution of p^ pomts is 
for z = 1, 2, . . . , r, then Wmograd has shown (Ref. 10) that 
the number of multiplications needed to compute a d-point 
cyclic convolution is equal to = - . . rriy. 

It is necessary to choose only certain values of d as the 
transform length in order to combme the Wmograd transform 
with the fast complex integer transform over GF(q^), where 
q =z 2P ~ 1 is a Mersenne prime. For this purpose, let the num- 
ber d have the form 

d = a*2^^p ( 2 ) 

where /zz = 0, 1, 2, 3 and zz = 3, 5, 7, or 9. For most practical 
applications, it suffices in (2) to let p = 31 or 61 . 

If d, the transform length of the cyclic convolution, is given 
(2), then by Theorem 1 , there exists a permutation of rows 
and columns so that the cyclic dXd matrix can be partitioned 
into blocks of (2"* ♦ p) X (2^^ • p) cyclic matrices in such a 
manner that the blocks form an zz X zz cyclic matrix. Now the 
cyclic convolutions of zz = 3, 5, 7, or 9 complex number points 
can be accomplished by Winograd’s algorithm. Since 2^ • 
plz?^ “ 1 for m = 0, 1, 2, 3, a transform of length 2^ • p over 
GF{q^') can be found and used to compute the cyclic convolu- 
tion of the 2"* - p complex number points. The number of 
multiplications needed to perform this convolution is 2"^ - p. 
Using this and the number of multiplications needed for Wino- 
grad’s algorithm, the total number of multiplications needed 
to perform a convolution of d complex number points can be 
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computed. The results of this computation are shown in 
Table 1. The present algorithm, and Agarwal and Cooley's 
algorithm are compared in this table by givmg the total num- 
ber of complex number multiplications required to perform 
the different algorithms. 


Example: Let d - 6. Next suppose that the input function 
defined by 

“ 1 +?0 « = 0, 1 


It was shown above that Winograd's algorithm can be com- 
bined with a transform over to yield a new rather fast 

hybrid algorithm for computing the cychc convolution of 
complex values. In this algorithm it was necessary to compute 
the cyclic convolution of 2^” * p complex number points for 
m = 0, 1, 2, or 3. This cyclic convolution of two <i-point se- 
quences of complex number points is 


•k 


d-1 

«=0 


( 3 ) 


where d\^ ^ p and {k - n) denotes the residue ofk- n mod d. 
To compute this convolution the components of the truncated 
complex numbers and/^ must be converted- first to integers 
and with dynamic ranges, say, A and 5, respectively. 
Previously (Ref. 5), a sufficient dynamic range constraint for 
A and B was shown to be 


A< 


4Bd 


If^ -E, (4) reduces to 


A< 



J 


( 4 ) 


( 5 ) 


= 0 


2<«<5 


is convolved with itself. This convolution is 


5 

n=0 


where (x) denotes the residue ^ - n of modulo 6. This convo- 
lution can be written in matrix form as 


i:\ 


\ 


I ^0 ^2 ^3 ^4 ^5 \ 

I flj Cj \ 


^2 ^3 ^4 ^0 


^4 «5 ‘^0 


“4 “s «0 “l “2 S 


*^5 ^0 ^2 ^3 ^4 




where [x] denotes the greatest integer less thanx. 

If the circular convolution of and is denoted by 
fork = 0, 1 , 2, . . . , *“ 1 , then using the procedure described 
in the example of Ref. 5, can be obtained by using fast 
transforms over GF{q'^), In (3), can be obtained by scaling 
back cj. into the scale of the original complex numbers for k = 
0, l,2,.,.<i^ 1. Evidently, the only error made in this com- 
putation of the c^s is the truncation error. 

The dynamic range constraint A of the input sequence 
given in (5) is generally very pessimistic. By an argument simi- 
lar to that used for integer convolutions (Ref. 14), one can 
lessen the' severity for this dynamic range constraint and still 
maintain in the interval ±(g - l)/2 with a small probability 
of overflowing. This assertion is justified in the Appendix. 

To illustrate this new hybrid algorithm, consider the follow- 
ing example. 


By Theorem 1, there exists a permutation K of rows and col- 
umns so that the above cyclic matrix can be partitioned into a 
2 X 2-block matrix of 3 X 3 cyclic blocks as follows: 


l’:\ I 




\ / *( 


^4 S ‘^3 


®2 ^4 **5 ^3 


«3 “l S "'O ^4 «2 


‘^1 “5 S '^ 4'='2 «0 


^5 ®3 ^1 ^2 S ^^4 j 


\ 
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This matrix equation, has the form 
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where 




. i 





II 


II 

^1 

1 

C4 

II 


' h 

1 

1 


f ' 




\V 





/ S^i"'s\ 



, ^ = 

1 


1 % «4 j 


i 




1 



r\ 


1 

1 

1 

0 

1 


' 0 



1 

■ i 


i>/ 


Ijdtx^ - = l andj^Q “ 1* 'Then the 

matrix equation defined in (7) can be obtained by computing 
the convolution of the two sequences and To do this, 
use a transform over In order to avoid overflow, one 

needs to choose ^ = 7 so that the integer components of a^, 
lie in the interval ±1. 


Since 2 is a 3rd root of unity, the transform over GF( 1 '^) of 


3-1 

^ = 2® + 2^'^ for fc = 0, 1 ,'2 

«=0 


Thus,XQ=2,Xj =5,Z^=3. 


Similarly, the transform of sequence;/^ is 


Thus 


!y\ I(A+B){X^+X^) + (A-B)(X^-X^)] 

(A +5) (Xj + X^) -(A-B) (Z, - X^) I 




= 2"^ 


= 2"^ 




/d+e\ 

[d-eJ 


( 6 ) 


3-1 

forfc = 0,l,2 

n=0 


That is, Yq = 2, = 6, = 6. But Dj, = = 4, 

= 2, 7)2 = 4. These are the only complex integer multiplica- 
tions needed to perform this transform. The inverse transform 
ofD^is 


3-1 
*=0 


-i-nk 


whfte D = (A+B) (Zj + Zj), E = (A - B) (Zj - Z^).- Now 







1 a^+a^.a^+a^.a^+a^ 

1 


1 


\a^+a^,a^^a^,a^+al 




= 5(4 + 2 • 2"^ + 4 ♦ 2“^*) mod 7 for fc = 0, 1, 2 


since 3"^= 5 mod 7. Thus finally, = 1 , dj = 2, dj- = 1: 

In a similar fashion, matrix E, given in (6), can also be ob- 
tained as = ly = -2, ^2 = 1. Thus, by (6), one obtains 
finally = 1, = 2, ^2 = 1, C3 = 0, = 0 and = 0. 
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Table 1. Complexity of New Algorithm for the Convolution of 
Complex Number Points 


d 

Factors 

No. Complex Number 
Multiplications 

No. Complex Number 
Multiplications of 
Agarwal and Cooley’s 
Algorithm 

120 



560 

124 

4X 31 

124 


210 



1520 

244 

4X 61 

244 


248 

8X 31 

248 


420 



3800 

488 

8X 61 

488 


744 

3X 248 

992 


840 



10640 

1260 



20900 

1464 

3X 488 

1952 


2520 



58520 

3720 

3 X 5 X 248 

9920 


7320 

3 X 5 X 488 

19520 




Appendix A 

A Probabilistic Dynamic Range Constraint for the Transform Over GF{q^) 


Let ~ and b„ = x„ Then the cyclic con- 

volution of and b„ is given by 


d-l 




«=0 


and 


E(xl) =E(yl) = 


12 


Substituting these values in A-2, 


d-l 




1-1 = A 




K 


(A-1) 


where = d„x„ - 0„y„, v„ = a„y„ + P„x„. In many applica- 
tions, the sequences a„, P„,x„, andj^ can be regarded as mu- 
tually independent. With tlus assumption, consider the sum 


»i=n 


O 


2 



(A-3) 


Now by the central limit theorem (Ref. 15), the probability 
6f exceeding a threshold X is 


F 



>\ 


= 2(1"V^(X)) 


where ip(X) is the standard normal distribution. This equation 
can be written as 


in (A-1). The expected value of is given by 

^ j = 2(l-^(X)) (A4) 

EQi^) = E(a„-)E(xJ-E(^^)E(y^) 

To keep from overflowings one needs the inequality 

where E denotes the expected value operator. With no loss 

of generality, the means of /?„, and can be ^ 

assumed to be zero. With these assumptions, EQx^) = 0, and — 

the variance <7^ of is given by 


< = E(fil) = E(ocl) E(xl) + E(f^) E(yl) (A-2) 


Finally, assume that and are uniformly distributed over 
the dynamic range A, and that andy„ are uniformly dis- 
tributed over the dynamic ranged. 


Hence if one sets 




g~ 1 
2 


then (A4) is the probability of overflow. 


(A-5) 




12 


For example, if X = 3, then the probability of overflowing is 
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P j 15^ l> X V^t j = 2(1 - Sff(3)) = 0.0026 


-2*i<a„,/3 <2*1 


which is very small. Substituting (A-3) into (A-5) yields 


-2*2’<x ,v <2*2 


A 


B=^kzJl 

y/2d\ 


(A-6) 


Equation (A-6) is the required relation among the parameters 
A, By q, dy and Similarly, one obtains the same result, 
defined in (A‘6) for the sum 




Let 2^1 and 2^2 be binary scale factors for and re- 
spectively. Then since 


the dynamic ranges A and B are and respec- 

tively. For most applications the twoMersenne prime 2^^ - 1 
and 2^^ ^ 1 will provide enough bit accuracy and dynamic 
range for computing two 2^ • p-point sequences of complex 
numbers. To illustrate this, Jf d is chosen to be = 2^, X = 3, 
q-2^^ ~ 1 and Atj = * 2 , then by (A-6) 


^= 2 *^ 1 *^^= / ^ 2 ^^ 

n/VMX 


Thus one needs approximately - ^2 satisfy 

(A-4) with an overflow probability equal to 0.0026. This is a 
considerably better bound than one obtains using formula (5) 
for A, In fact, the dynamic range constraint (5) yields A = 
2^1 :=; 2 ^ ^ 
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One of the most exciting challenges facing gravitational theoreticians and experi- 
menters in the remaining decades of this century will be the search for ''gravitational 
waves” as predicted by Einstein *s General Theory of Relativity, Proposals have been 
advanced to search for gravitational waves in ultraprecise two-way Doppler data. In such 
an experiment, the total measurement system includes the Deep Space Network Tracking 
System, the spacecraft, the intervening media, and the data processing system. Pre- 
liminary estimates of gravitational wave characteristics are used to define a baseline 
experiment, with a total measurement system fractional frequency fluctuation of 
1 X and a desirable experiment, with a total measurement system fractional 

frequency fluctuation of 1 X 

The experiment to detect gravitational waves in ultraprecise two-way Doppler data is 
described, as are the anticipated requirements for the Deep Space Network, the space- 
craft, and the data processing system. The article concludes by describing the steps 
necessary to provide the capability to perform this experiment. 


I. introduction 

One of the most exciting challenges facing gravitational 
theoreticians and experimenters in the remaining decades of 
the 20th century will be the attempted detection and measure- 
ment of ''gravitational waves” as predicted by Einstein's 
General Theory of Relativity. Gravitational waves are spatial 
strains propagating at the speed of light which are (expected to 
be) generated by the violent collapse of stars, star clusters, 
galaxy cores, quasars, and seyfert galaxies into supermassive 
black holes (for instance, see Thorne and Braginsky, Ref. 1). 


Proposals have been advanced to use “ultraprecise” two-way 
Doppler data while tracking spacecraft at large distances (i.e., 
> 1 AU) to detect very-Iow-frequency Hz) gravita- 

tional waves (for instance, see Davies, Ref. 2, and Estabrook 
and Wahlquist, Ref. 3). Although such a scheme appears to be 
one of the more promising suggested, substantial improve- 
ments in current Deep Space Network (DSN) and spacecraft 
performance will have to be achieved, as well as provision of 
new capabilities, to bring to fruition the utilization of ultra- 
precise Doppler data in the search for gravitational waves. It is 
thus the purpose of this article to provide a brief description 


100 



of the experiment to detect gravitational waves in ultraprecise 
Doppler data, and to establish anticipated requirements for the 
spacecraft, DSN, and subsequent data' processing. 


II. Definition of “Baseline” and “Desirable” 
Experiments 

As already noted, gravitational waves are expected to be 
produced by the violent collapse of stellar bodies into super- 
massive black holes; these waves are expected to be evidenced 
as spatial strains propagating at the speed of light. Specifically, 
one expects (Refs. 1 and 3) that the waves will: 

(1) Alter the distance between separated free masses. 

(2) Produce a fractional frequency shift in ultraprecise 
Doppler data, on the order of the wave amplitude. 

Thorne and Braginsky (Ref. 1) have estimated possible 
gravitational wave characteristics; these are presented in Fig. 1 . 
As can be seen, the wave duration (r) is proportional to the 
wave amplitude {h) and the expected ranges of durations and 
amplitudes^ are: 

40 seconds < r < 40,000 seconds 
10'^'^ <h< 

In Section III it will be shown that the scheme to utilize 
ultraprecise Doppler data requires : 


r < Round-Trip-Light-Time (RTLT) 

Since the next decade of funded and proposed deep space 
missions have maximum distances of either Jupiter (Gahleo 
Mission, Solar Polar Mission) or Saturn (Saturn Orbit er 
Mission), it is immediately apparent from Fig. 1 that one is 
constrained to search for waves in which the amplitude is 
'•-lO':^^ or less. It is also obvious that one would like to search 
for waves with amplitudes as low as (r = 40 seconds). 

This then motivates the selection of a “baseline’ ' experiment 
to look for the strongest expected waves, and a “desirable” 
experiment to search for a wide variety of waves. Using this 
rationale, and noting that there are stringent requirements 
throughout the entire measurement system, consisting of 

(1) DSN Tracking System 

(2) Intervening media 


should be emphasized that these are highly tentative estimates 


(3) Spacecraft 

(4) Data processing 

one proceeds to define the baseline experiment as follows: 

(1 ) Search for waves with amplitude ^ lO" ^ ^ . 

(2) Require total measurement system fractional frequency 
fluctuation (o(AF/i0) of 10”^ ^ or less.^ 

(3) Require each identifiable independent element of the 
total measurement system to have a fractional fre- 
quency fluctuation of 3 X 10"^^ or less.^ 

(4) Require averaging times (rj between 50 and 5000 
seconds. 

Similarly, one defines a desirable experiment as follows: 

(1) Search for waves with amplitudes /z ^ 10“^'^. 

(2) Require total measurement system fractional frequency 
fluctuation (o(AF/F)) of 10"^ or less. 

(3) Require each identifiable independent element of the 
total measurement system to have a fractional fre- 
quency fluctuation of 3 X 10“^^ or less. 

(4) Require averaging times (r^) between 5 and 5000 
seconds. 

While it is considered technologically possible to achieve a 
baseline system in the 1980s, the desirable system may not be 
achievable until the 1990s, if then. 

While fractional frequency fluctuation is most appropriate 
for the basic description of ultraprecise two-way Doppler 
(frequency) data, certain elements of the measurement system 
are more conveniently thought of rin terms of phase fluctua- 
tion (a(A0)) or path length fluctuation (a(AL)). For the 
remainder of the article, the following equivalents will be 
assumed: 

o(A^) = a(AF/i0 * d^jdt • 
a(AL) = o(AF/i0 * c • 


^The symbol a will be used to represent a “generalized” measure of 
fluctuation which ultimately may be determined to most appro- 
priately be the Allan variance, RMS, etc- 

^ Since there are numerous independent elements, it is assumed appro- 
pnate to require each independent element to be specified at a level 
one-half order of magnitude below the total measurement system 
requirement. 
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where 

0 = phase 
F = frequency 
t = time 

= averaging time 

L = length 
c = speed of light 

Some typical conversions are 

(1) a(AF/iO = lX 10-^^ 


100 1000 

Tg seconds seconds 

o(A0), X-band 0.3 deg 3.0 deg 

o{AV) 0.03 mm 0.3 mm 

(2) a(AF/i0=l X 10-1 

100 1000 

seconds seconds 

o(A0), X-band 0.003 deg 0.03 deg 

a(AL) 0.0003 mm 0.003 mm 


III. Experiment Description 

The special feature which allows the usage of ultraprecise 
Doppler data for the possible detection of gravitational waves 
is a unique “three-pulse’' signature which is a function of the 
spacecraft, Earth, and gravitational wave propagation direction 
geometry (Ref. 3). The pulses (fractional frequency shifts) 
result from effects which' are conveniently described as 
follows: 

(1) “Clock speed-up” (Earth only effect) 

(2) “Buffeting” (equal Earth and spacecraft effect) 

Consider a case where 6 defines the angle between the 
gravitational wave propagation direction and the Earth- 
spacecraft line, as seen in Fig. 2. Clock speed-up is evi- 
denced as pulses of maximum amplitude -h and +/z at the 
time the wave impinges on Earth and a RTLT later, respec- 
tively. This is shown schematically in Fig. 3. Similarly, Earth 


buffeting is seen as pulses of maximum amplitude +77 cos Q 
when the wave impinges on Earth and a RTLT later. In 
addition, spacecraft buffeting produces a pulse of - 2/2 cos 9 in 
the Doppler data a one- way-light-time (OWLT) after the wave 
strikes the spacecraft. Earth and spacecraft buffeting effects 
are schematically illustrated in Fig, 4. Finally, combining the 
clock speed-up and Earth and spacecraft buffeting effects, one 
obtains the unique three-pulse signature which is expected to 
be seen in ultraprecise two-way Doppler data. The charac- 
teristics of this signature are as follows: 

(1) Pulse amphtudes 

(a) / 2 (cos 0-1) 

(b) /z(“2cos0) 

(c) /z(cos 9 + 1)' 

(2) Pulse separation times (r =' Earth- spacecraft distance) 

(a) “(1 + cos 0) 

(b) -^(1 - cos 0) 

Figure 5 schematically demonstrates the combination of clock 
speed-up and buffeting effects for 0 = 60 degrees. 

Finally, it is noted that the above exercise is for a single 
spacecraft measurement. Should dual spacecraft measurements 
be made, one could obtain the two-dimensional propagation 
direction of the gravitational wave — certainly a most valuable 
bit of additional information. 


IV. Anticipated Gravitationai Wave 
Detection Experiment Requirements 

The major categories identified as constituents of the 
gravitational wave detection measurement system are: 

(1) DSN Tracking System 

(2) Intervening media 

(3) Spacecraft 

(4) Data processing 

In general, requirements are necessary for multiple 
independent elements in each category. As already noted, each 
independently identifiable element will be “specified” at a 
performance level 1/2 order of magnitude better than the 
entire system requirement, under the assumption that they 
will combine m a fashion similar to root-sum-square. Each 
major category is discussed below. 
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A. DSN Tracking System 

The central element of the DSN Tracking System is 
obviously the frequency standard of the DSS Frequency and 
Timing Subsystem. The performance of the new hydrogen 
masers currently being implemented is considered to be 

a(AF/F) ^ 3 X 10"^ ^ 

> 100 seconds 

so that a decrease of about an order of magnitude will be 
required to achieve the baseline experiment defined in Section 
II: 

3 X 10-1^ 

50 seconds < < 5000 seconds 

Numerous other DSN Tracking System elements will have 
to be evaluated as to their current performance, and the 
technology developed to obtain the desired frequency 
stability. Below are listed the various pertinent subsystems, 
with the appropriate performance parameter indicated in 
parentheses: 

(1) DSS Frequency and Timing Subsystem (Frequency) 

(a) Hydrogen maser (or other frequency standard) 

(b) Frequency distribution 

(2) DSS Receiver-Exciter Subsystem (Phase) 

(a) Closed-loop receiver 

(b) Doppler extractor 

(c) Exciter 

(3) DSS Tracking Subsystem (Phase) 

Metric Data Assembly 

(4) Antenna Mechanical Subsystem (Length) 

(5) Antenna Microwave Subsystem (Phase) 

(6) Transmitter Subsystem (Phase) 

(7) System Cabling (Frequency) 

Table 1 summarizes the required performance for the above 
elements for both the baseline and desirable experiment cases. 


B. Intervening Media 

In passing between the Deep Space Station (DSS) and 
spacecraft, the signal passes through the following interactive 
media: 

(1) Solar Wind (dispersive) 

(2) Troposphere (neutral) 

(3) Ionosphere (dispersive) 

In the discussion to follow, the ionosphere, consisting of 
charged particles, will be considered jointly with the. similar, 
but generally much larger effect of the Solar Wind. 

(1) Solar Wind. The Solar Wind effect on the gravitational 
wave detection experiment has been explored in detail 
in a previous article (Ref. 4). The results from Ref. 4 
are summarized as follows (for - 1000 seconds): 

(a) Current S-band (both uplink and downlink) Solar 
Wind limitation: 

o{AF!F) = 3X 10-^3 

(b) Expected optimum X-band (uplink and downlink) 

Solar Wind limitation: ' 

o^^FlF) = 1 X 10-^^ 

It is a major conclusion of Ref. 4 that to achieve 
a Solar Wind fractional frequency fluctuation of 
a(AF/F) = 3 X 10“^^, one would require simul- 
taneous' S- and X-band uplink and downlink frequency 
band capability. To achieve the desirable experiment 
level (3 X 10"^^), it is speculated that the spacecraft 
will require sophisticated new capabilities, such as an 
onboard frequency standard and a Doppler counter. 

(2) ‘"Wet” Troposphere, An extremely preliminary esti- 
mate (Ref. 5) for the wet troposphere at an elevation 
of 30 degrees (20 cm total signal retardation assumed) 
and an averaging time = 1000 seconds is: 

a(AL/Z) ^ 1 X 

It is considered that to meet the baseline experiment 
requirements, water vapor radiometer calibration of the 
wet component of the troposphere will be required. 
For the desirable experiment level, new spacecraft 
capabilities are envisioned. 
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(3) ‘‘Dry” Troposphere. Again, an extremely preliminary 
estimate (Ref. 5) for the dry troposphere at an eleva- 
tion of 30 degrees (400 cm total signal retardation 
assumed) and an averaging time = 1000 seconds is: 

a(AZ,/L) ^ 7X lO-is 

It is considered that to meet the baseline experiment 
requirements, extremely accurate surface barometric 
pressure calibration of the dry component of the 
troposphere will be required. 

Media calibration capabilities required for the baseline and 
desirable experiment are summarized in Table 2. 

C. Spacecraft 

Similarly to the DSN Tracking System, individual critical 
path elements within the Spacecraft Radio Subsystem, such as 
the receiver, transponder, transmitter, need to be specified at 
o{AFIF) levels of 3 X 10“^^ and 3 X 10“^ ® for the baseline 
and desirable experiments, respectively. 

D. Data Processing 

Data processing numerical accuracies will need to be com- 
mensurate with the a(AF/iO levels of 3X10“^^ and 
3 X 10“^^ for the baseline and desirable experiments, respec- 
tively. Translated into units of cycles, the appropriate require- 
ments for data processing numerical accuracies are as follows: 

(1) Baseline Experiment 

(a) S-band - 3 X lO"^ cycle * (rjSO) 

(b) X-band - 1 X 10“"^ cycle • (rj50) 

(2) Desirable Experiment 

(a) S-band — 3 X 10“® cycle • (rj5) 

(b) X-band — 1 X 10"'^ cycle • (rj5) 


V. Discussion and Summary 

In order to progress towards a gravitational wave detection 
capability in regard to the proposed usage of ultraprecise 
Doppler data, a series of actions will need to be undertaken, as 
follows: 

(1) Development of gravitational wave detection experi- 
ment “drivers,” such as: 

(a) o(AFlF) = 3 X 10"^^ frequency standard 

(b) X-band uplink capability 

(2) Identification of all critical path elements in the Space- 
craft and DSN Tracking System, including evaluation 
of current frequency stability performance and 
development of the appropriate technology to achieve 
a fractional frequency fluctuation of 3 X 10"^ 

(3) Additional study of the fractional length fluctuation of 
the troposphere, and development of the appropriate 
tools necessary to calibrate the troposphere to 
o{ALlL) = 3X 10-1 ^ 

Prior to obtaining a total measurement system of 
g^AF/F) = 1 X 10“i ^ , it is suggested that “demonstrations” in 
the early 1980s might be possible at about the o^AFjF) ^ 
1 X 10“!^ level. Significant features of such a demonstration 
would be: 

(1) a{AFIF) « 3 X frequency standard'^ 

(2) X-band uplink and downlink 

Possible missions for such a demonstration would be the 
Solar Polar and the Galileo missions. 
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Table 1. Required performance summary 



Baseline 

Desirable 

Element 

Experiment 

Experiment 


Requirement 

Requirement 

Total system, a(AF//0 

IX 10“^® 

1 X 10-11 

DSN Tracking System, a(AF/iO 

< 1 X 10~‘5 

< 1 X 10-11 

Averaging time, seconds 

50 < < 5000 

5 <T^< 5000 

Frequency and Timing Subsystem 



Frequency standard (H2 maser), 
GiAFjF) 

3 X 10“^® 

3 X 10"1® 

Frequency distribution, u{AFlF) 

3 X 10“!*^ 

3 X 10-1® 

Receiver-Exciter Subsystem 



Closed-loop receiver, a(A<^>), X-band 

0.05' • {ta/50} 

0.00005° • {tJS} 

Doppler extractor, a(A0), X-band 

0.05° ■ {t„/50} 

0.00005° • {tJ5} 

Exciter, tr(A0), X-band 

0.05' . {ra/50} 

0.00005° • {TglS} 

Antenna Mechanical Subsystem, o{AL) 

0.005 mm • {ra/50} 

0.000005 mm • {ra/5} 

Antenna Microwave Subsystem, a(A0), 

0.05° • {tJ50} 

0.00005° • {ra/5} 

X-band 

- 


Transmitter Subsystem, a(A0), X-band 

0.05° . {t„/50} 

0.00005° • {t<j/5} 

System cabling, a(AF/F) 

3 X 10-1® 

3 X''10“1® 

Data processing software 



Numencal accuracy, X-band 

1 X 10“^ cycle ♦ {rfl/50} 

1 X 10~1 cycle • {t^S} 

Spacecraft 



Radio Subsystem, a(AF/F) 

3 X 10~1® 

3 X 10-1® 


ORIGINAL PAGE IS 
AF POOR QUALITY 
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Table 2* Required capabilities summary 


Media Effect 

Baseline Experiment 
Requirement 

Desirable Experiment 
Requirement 

Solar Wind 



• DSN 

Receive and transmit simultaneous S- and X-band 

Multiple one-way and two-way Doppler links 

• Spacecraft 

Receive and transmit simultaneous S- and X-band 

On-board frequency standard and Doppler counter; 


- 

Multiple one-way and two-way Doppler links 

Wet troposphere 



• DSN 

Water vapor radiometer 

Multiple one-way and two-way Doppler links 

• Spacecraft 


On-board frequency standard and Doppler counter; 
Multiple one-way and two-way Doppler links 

Dry troposphere 



• DSN 

Surface barometric pressure 

Afultiple one-way and two-way Doppler links 

• Spacecraft 


On-board frequency standard and Doppler counter; 
Multiple one-way and two-way Doppler links 
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Fig. 1. Relationship between gravitational wave amplitude and 
burst duration 



Af/f V , 

t / 

CLOCK W 
SPEED-UP—^ 


CLOCK SPEED-UP 
(RETURN) - 


(DIAGRAM AFTER ESTABROOK) 
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Fig. 2. Gravitational wave geometry 


Fig. 5. Unique gravitational wave signature in ultraprecise 
Doppler data 
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This article presents the preliminary findings of a study being made at CTA-21 to 
determine whether Deep^ Space Station control room equipment power may he turned off 
to conserve energy. The results of reliability analysis indicate that there may be some 
Correlation between the observed increase in failure rate and cycling of equipment power 
in the eight-month study period. 


i. Introduction 

The DSN Energy Conservation Project has been investigat- 
ing modifications to Deep Space Stations at Goldstone, 
California, in regard to energy conservation and cost reduc- 
tion. 'One of these modifications is to selectively turn off 
certain electronic equipment during periods of nonuse or low 
activity. Significant energy savings are expected to result, 
especially when equipment control is automated in the future 
via the computer-based Utility Control System. However, a 
key question in considering this procedure is whetlier cycling 
of electrical power will change equipment reliability and thus 
impact DSN maintenance and operations activities or cost. 

The DSN Equipment Compatibility Test Area (CTA-21) at 
JPL, Pasadena, was selected as a test site. This is the CTA-21 
control room m building 125 where electronic and air condi- 
tioning equipment is presently turned off on weekends. The 
use of the CTA-21 control room provides field observations of 
equipment reliability under conditions closely representing a 
Deep Space Station’s environment. 


Tentatively the study is divided into two phases. Phase one, 
which is the substance of this article, consists of data 
collection and performance of a preliminaiy analysis. It covers 
the period of April through November 1977, examining overall 
failure trends. Phase two is a detailed correlation analysis of 
Mures with possible causes. Each phase will be discussed in 
detail later. 

II. Field Observations and Data Collection 

After examining CTA-21 past failure and maintenance 
records, it was decided that the existing records were not of 
sufficient detail or comprehensiveness to support the study. 
Therefore, in conjunction with CTA-21 managemenTand staff, 
a special equipment event log was implemented. This form, as 
shown in Fig. 1, provides both the DSN Energy Conservation 
Project and CTA-21 with an unique source of information 
consisting of: 

(1) Equipment failures, including references to equipment 
event reports (EER’s) when applicable 
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(2) Power turnoff periods 

(3) Commercial power outages 

(4) Air handler failures 

(5) High ambient control room temperatures 

(6) Equipment modifications through engineering change 
orders (ECO’s) 

(7) Other installations, removals, calibrations, and main- 
tenance 

The log utilizes two coding groups to assist in data 
processing and interpretation. The first group is the event 
code, which specifically identifies the nature of the event 
reported. The second group code provides categorization of 
the type of equipment involved. The equipment event log 
serves not only the purposes of the study but also makes 
available to CTA-21 management and engineering a useful 
maintenance record. Since outstanding events are referenced in 
subsequent order in the log, relevant actions can be traced and 
monitored for optimum performance. 

The equipment at CTA-21 represents nearly ninety percent 
of the equipment located in the control room of a typical 
Deep Space Station. At times, this consists of over one 
hundred individual cabinets of widely varying types of 
electronic equipment. Therefore, the equipment was divided 
mto ten system categories to simplify analysis: 

(1) Telemetry (6) Radio metric 

(2) Communication (7)- Simulation and test 

(3) Radio frequency (8) Timing 

' (4) Command (9) Facility 

(5) Monitor and control (10) Software 

No special turnoff procedure is used at CTA-21. However, 
individual cabinets and systems are turned off in no particular 
sequence before general power is removed. The control room 
air handlers are allowed to run an additional thirty minutes 
after equipment turnoff to remove residual heat. Tliis proce- 
dure was initiated when a significant rise in temperature- 
related failure rates was noted prior to the preliminary study 
period. Likewise, individual equipment is turned on after 
general power is restored. 

HI. Data Study and Analysis 

Initial CTA-21 equipment event data have been processed 
for the thirty-four-week period from April 4, 1977, until 
November 23, 1977, and were divided into seven subperiods 
designated as either "experimental” or “control” periods. 


Experimental periods are defined as those -time intervals during 
which equipment power was turned off on weekends. That is, 
those contiguous groups pf weeks where .equipment power is 
off for nominally two days (weekends) and equipment power 
is on for five days. In control periods, the power was on all 
seven days -of -the- week. The controf periods, listed in Table 1, 
start with a weekend during which power was left on and 
cover the period until power is next turned off. 

Only failure and power turnoff events were included in this 
preliminary phase one analysis. No attempt was made to 
determine possible causes of individual failures. Turnoff of 
power at CTA-21 is conducted on a noninterference basis with 
the operations schedule, which leads to widely varying study 
period lengths. The statistics of hourly failure frequency (or 
rate) was computed to allow comparison of control and' 
experimental study periods. 

The distribution of daUy failure rates occurring in each 
subperiod over the entire study period is shown in the 
time-line bar charts in Fig. 2. The time-line bar chart was 
produced to reveal any significant distribution patterns, and to 
display overall trend. Tables 1 and 2 compare failure rates to 
indicate any significant effect of power turnoff on equipment 
rehability. 

IV. Results of the First Phase Study 

Inspection of the above time-line bar chart does not 
directly indicate any particular repetitive pattern. However, 
the overall 3 4- week failure trend was downward as shown by 
Table 1. Furthermore, the failure rates in both period types 
decreased over the duration of the study. One explanation 
could be that turning off equipment power over long time 
periods extends equipment life and is manifest in* lowered 
failure rates. Other possibilities include reduced station 
activity and unreported failures. These hypotheses will be 
examined in more detail during the second phase. 

Table 1 indicates that failure rate increased during periods 
of weekend power turnoff, and that the overall experimental 
period failure rate^ exceeds the aggregate control period failure 
rate by 73%. Table 2 shows that the majority of failures 
occurred in the telemetry and communication systems equip- 
ment. However, the greatest sensitivity to power turnoff was 
exhibited in the monitor and control system, although all 
systems show an increase in failure rate during experimental 
periods. 

Figure 3 shows the time distribution density of failures as a 
function of elapsed time. The distribution for control periods 
exhibits a somewhat uniform density with considerable ran- 
dom variation. Since CTA-21 is normally manned during the 
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week, usually no observations are made of failures that occur 
on weekends. They are reported when- discovered during the 
five-day work week. This explains the marked decrease in 
failures on weekends. The distribution of failures after power 
turnoff shows an apparent 2.5^day cycle of unknown signifi- 
cance, whereas a simple “one cycle” exponential distribution 
would tend to suggest a strong relationship between equip- 
ment failure and power Turnoff This phenomenon will be 
further explored in detail during phase two of the study. 


V. Significance of Results 

The initial phase of this study shows an increase in faflure 
rate during weekend turnoff periods. However, the results are 
based on only 282 failures occurring over 34 weeks. The 
adequacy of the data sample was tested, and data confidence 
intervals were determined to the 95% confidence level as 
shown in Table 3. 


The resulting confidence intervals indicate tliat the com- 
puted failure rates are at best rough approximations. With an 
objective of reaching a 20% confidence interval (±10% 
confidence limit) it was determined that a possible 13 
additional months of data would be required. Since new data 
may have a major effect upon final results, conclusions at this 
stage are undecisive. 

The next phase of study not only will benefit from a larger 
data base, but also detailed correlation analysis will be 
performed. Correlation analysis will indicate relationships’ 
between failures and other events such as air handler malfunc- 
tions, liigh ambient temperatures, equipment modifications, 
life-cycle, etc. Failures thus identified may be eliminated from 
data being investigated for relationsliips to turnoff periods. 
The objective of phase two will be to identify equipment that 
is sensitive to power turnoff, thereby providing a criterion by 
which electronic equipment may be selectively chosen for 
power shutoff to conserve energy. 
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Table 1, Total CTA-21 failure rate for 34-week study period 





Total 

hours 

Pwr on 
hours 

Failures/hour 

Period type 

Weeks 

Failures 

■Control 

Expt 

Control 

8 

72 

1344 

1344 

0.0536 


Expt l'’ 

8 

91 

1344 

875 


0.1040 

Control 2 

6 

52 

984 

984 

0.0528 


Expt 2 

5 

29 

864 

563 


0.0515 

Control 3 

2 

11 

336 

336 

0.0327 


Expt 3 

1 

7 

168 

118 


0.0596 

Control 4 

4 

20 

624 

624 

0.0321 


Total percentage 

34 

2^ 

5664 

4844 

0.0471 

100% 

0.0816 

173.1% 


^Control periods (eqpt. not turned off): 1, April 2-May 27. 2. July 23-September 1; 3. October 8-October 21; 4. October 29-November 23. 
^Experimental periods (eqpt. turned off on weekends): 1. May 28-July 22; September 2-October 7; 3. October 22-October 28. 


Table 2. Failure rates for the ten system categories 




Percent 
of total 

Failures/hour 


Percent 

difference 

System 

Failures 

Control period 

Expt period 

Telemetry 

85 

30.1 

0.0143 

0.0244 

70.8 

Communication 

58 

20.6 

0.0110 

0.0141 

29.1 

RF 

34 

12.1 

0.0058 

0.0096 

66.8 

Command 

31 

11.0 

0.0055 

0.0084 

52.6 

Monitor 

31 

11.0 

0.0040 

0.0116 

192 6 

Metric 

16 

5.7 

0.0024 

0.0052 

111.3 

Simulation 

14 ‘ 

5.0 

0.0024 

0.0051 

58.5 

Timing 

11 

3.9 

0.0015 

0.0039 

153.6 

Facility 

2 

0,7 

0.0003 

0.0006 

111.3 

Software 

0 

NA 

NA 

NA 

NA 

Total 

282 

100% 

0.0471 

0.0816 

+73.1% 


Table 3. Data confidence and sample adequacy 

Period type 

Failuies/hour 

Confidence mterval 

Limit 

% of rate 


Control 1 

0.0536 

±0.0152 

±28.5% 

Expt 1 

0.1040 

0.0260 

25.0 

Control 2 

0.0528 

0.0167 

31.5 

Expt 2 

0.0515 

0.0232 

45.1 

Control 3 

0.0327 

0.0269 

82.1 

Expt 3 

0.0596 

0.0222 

37.4 

Control 4 

0.0321 

0.0161 

50.4 

Total control 

0.0471 

±0.0090 

±19.1% 

Total expt. 

0.0816 

±0.0174 

±21.4% 
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EQUIPMENT EVENT LOG 
STATION: CTA-21 

DATE: ^01/77 DAY: MON 


EVENT CODE: 

8 ^BROKEN E = ECO M=MISC 

C = CALIBRATION F = FAILURE 
D = DEFECTIVE I INTERMITTENT 


GROUP CODE: 

A = ANALOG FAC * FACIUTIES 

CA^ CABLES L = LOGIC P = PWR SUPPLY 

CO CONNECTOR M«=K€CH 0= OTHER 


LOG NO. 
REF. DATE 

GMT 




WA 

EVENT 

CODE 

GROUP 

CODE 

HISTORY 

LOG 

DESCRIPTION 
OF EVENT 

ACTION 







B 

B 


SI A #17 (UNUSED PORT) FAIIED 
WRAPAROUND TEST 


1230 

6/10 

2130 


PRACOMP 

INTERFACE 



B 

I/O 

METRIC 

COMPUTER INTERFACE FAILED WRAP- 
AROUND TEST 

UNIT FROM 14, WORKED 
ON BY DAVIES (SPARE) 

1231 

6/30 

1930 


COMM BUFF 3 



B 


COMM 

C5#3 DF Ll GHT FAI LS TO U GHT . 

UNIT APPEARS TO PASS DATA PROPERLY 

LED ON ORDER 
07/29/77 

1232 

7/17 

1400 


DIS 1 

MAG TAPE 




B 

MONITOR 

READ ERROR ON BOTH MAG TAPE UNITS 

MASTER OK 

1233 

7/21 

1800 





F 

B 


L/P INOPERATIVE 


1234 

7/21 

1100 





E 

L/0 


ECO 77.177 STARTED 


1235 

7/21 

0300 


TPA3 

SEQ DECODER 



F 

B 

TELEM 

SEQ DECODER INOPERATIVE 


1236 

08/01 

1400 


TCP PAPER 
PUNCH 

m 


B 

fm 

TELEMETRY 

UNABLE TO PUNCH TAPE, ERRATIC 

CORRECTED 

1237 

8/01 

1430 



m 



B 

COMM 

PREVENTIVE MAINT. PERFORMED ON 
PRINTER/PLOTTER 


1238 

8/01 

1500 


TCP PHOTO 
READER 

\ 

% 







1239 

8/01 




m 


B 

0 

SOFTWARE 


ERROR TO ECO 
76.042 

1240 

8/01 

2000 



m 


B 

FAC/O 

RF 

CLOSEOUT OF ECO 77.027. INSTALL 
RF link JPLTO HUGHES 



Fig. 1. CTA-21 equipment event log 
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FAILURES FAILURES FAILURES 





Rg. 2. CTA-21 equipment failure time-line chart 
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The DSS Radio Science Subsystem — Data Handiing of 
Very Long Baseiine Interferometry (VLBI) Data 

A. L. Berman 

TDA Engineering Office — - 

The DSS Radio Science Subsystem, originally implemented- to provide the data 
handling capabilities for the DSN Radio Science System, will be modified and augmented 
to provide similar capabilities fox the newly created DSN Very Long Baseline 
Interferometry (VLBI) System, This article describes the key characteristics, functional 
requirements, and operation of the DSS Radio Science Subsystem (DRS) as they pertain 
to usage of the DRS by the DSN VLBI System, 

\ 

I. Introduction 

The DSN Very Long Baseline Interferometry System was 
recently instituted as a DSN data system, following a 
successful review of DSN VLBI System requirements on 
February 28, 1978. Very Long Baseline Interferometry is a 
method of measuring the difference in the time of arrival of a 
signal produced by a distant natural radio source at two or 
more mutually distant antennas via simultaneous signal recep- 
tion and recording. Because of the difference in signal ray 
paths between the antennas, the signal will evidence small time 
delays in reception between any pair of antennas. By cross 
correlating the recorded signals, the time delay and its 
derivative may be established for each pair of antennas. 

Because the natural radio source is extragalactic, and hence 
can be considered a fixed object, the measured time delay and 
Its derivative provide information on the earth’s motion and 
the baseline vectors between the various antennas. The DSN 
intends to use VLBI to determine Universal Time one (UTl; 
the instantaneous rotational angle of the earth), polar motion, 
and the relative positions of the Deep Space Stations (DSS), as 


w.:: 2 S the time offsets and rates of change between the 
various DSS clocks. 

The DSN VLBI System will be a phased implementation 
with the first two phases identified as “Block 1” and 
“Block 2.” Block 1 VLBI will provide for clock synchroniza- 
tion between the DSS and is characterized by usage of strong 
radio sources and a relatively low data rate of 500 kbits/s. 
Block 2 VLBI will provide the radio source catalog for Block 1 
operations, as well as the earth motion and baseline vector 
parameters required for the Block 1 data processing. Block 2 
VLBI can operate with weak radio sources as well as strong 
sources, and is characterized by a relatively high data rate of 
32 Mbits/s. 

The DSS Radio Science Subsystem provides the data 
handling capabilities for the recently created DSN Radio 
Science System, and was described in detail in a previous 
article (Ref. 1). VLBI data is quite similar to radio science data 
in that both are generated via open loop receivers, and hence 


115 



the DSS Radio Science Subsystem will be modified and 
augmented to additionally provide the data handling capabil- 
ities for the. DSN VLSI system. This article provides a 
description of the key charactenstics, functional requirements, 
and operation of The DSS Radio Science Subsystem in specific 
regard to DRS usage by the DSN VLBI System. 


II. Functional Description of the DSS Radio 
Science Subsystem 

A. Definition 

The DSS Radio Science Subsystem, an integral element of 
the pSN VLBI System, performs the following functions: 

(1) Receives, digitizes and records narrowband open-loop 
receiver data, and transmits via wideband data line to 
the Network Operations Control Center (NOCC). 

(2) Receives, digitizes and records wideband open-loop 
receiver data, and transmits via magnetic tape to the 
NOCC. 

(3) Receives calibration and ancillary data and transmits in 
real-time via High-Speed Data Line (HSDL) to the 
NOCC. 

DSS Radio Science Subsystem functions and interfaces are 
shown in Fig. 1, while Fig. 2 presents functions and data flow. 

B. Key Characteristics 

tern are listed below for functions ox me uuL.uitat.iw.. — 
Assembly (ODA) and the Digital Recording Assembly (DRA). 

(1) VLBI Block 1 (ODA) Key Characteristics 

(a) Hardware and software compatibility with the DSN 
Mark III configuration. 

(b) A to D conversion of narrowband open-loop 
^ receiver data with 1-bit quantization. 

(c) Total data rate of 500 kbits/s. 

(d) Merging of open-loop receiver data with calibration 
and ancillary data (“VLBI data”) in real-time. 

(e) Reconstruction of analog signals for verification of 
calibration tones. 

(f) Temporary storage of VLBI data on magnetic tape 
at the DSS with subsequent transmission via 
Wideband Data Line (WBDL) to the NOCC. 

(g) Real-time transmission of VLBI calibration and 
status data to NOCC via High-Speed Data Line. 


(2) VLBI Block 2 (DRA) Key Characteristics 

(a) Hardware compatibility with the DSN Mark III 
configuration. 

(b) Provision for recording eight wideband open-loop 
receiver channels of 4 Mbit/s data rate each. 

(c) Provision for time tagging of recorded data to one 
microsecond or better resolution. 

(d) Reconstruction ‘of analog signal for verification of 
calibration tones. 

(e) Automated control of the DRA by the ODA. 

C. Functional Operation 

1. VLBI Block 1 Phase 1. Natural radio source signals are 

received by the Block IV closed-loop receivers and processed 
by the Advanced System Equipment (ASE). Multiplexed 
signals from the ASE are received by the VLBI Converter 
Subassembly and digitized. The ODA Modcomp Computer 
receives the data from the VLBI Converter Subassembly and 
records it on magnetic tape. Concurrently, the ODA receives 
calibration and ancillary data from the DSS Tracking Subsys- 
tem (DTK). In real-time the ODA formats the calibration and 
ancillary data for HSDL transmission to the NC Radio Science 
Subsystem (NRS) via the Communications Monitor and 
Formatter (CMF). Post pass, the ODA formats the digitized 
receiver data, calibration data, and ancillary data for Wideband 
Data line transmission to the GCF Data Records Subsystem. 
During VLBI operations, a reconstructed analog signal is 
passed to the Spectrum Signal Indicator ciificailon 

j The DSS 

woiiiigured and controlled by the 
DSS Monitor and Control Subsystem (DMC), while the DRS 
status is routed to both the DMC and NRS. Figure 3 presents a 
functional block diagram of the ODA while performing VLBI 
Block 1 Phase 1 operations. 

2. VLBI Block I Phase 2. Natural Radio source signals are 
received by the Wideband Multi-mission Receiver (MMR) and 
pass through VLBI Intermediate Frequency (IF) Converters. 
The ODA configures both the VLBI IF Converters (Bandwidth 
selection) and the Wideband MMR synthesizer (fixed fre- 
quency information). From the VLBI IF converters, the 
signals are received by the VLBI Converter Subassembly, 
where 'they are multiplexed and digitized. The Occultation 
Data Assembly Modcomp computer receives the data from the 
VLBI Converter Subassembly and records it on magnetic tape. 
Concurrently, the ODA receives calibration and ancillary data 
from the DSS Tracking Subsystem (DTK). In real-time the 
ODA formats the calibration and ancillary data for High-Speed 
Data transmission to the NC Radio Science Subsystem via the 
Communications Monitor and Formatter. Post pass, the ODA 
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formats the digitized receiver data, calibration data, and 
ancillary data for Wideband IJata line transmission to the GCF 
Data Records Subsystem. During VLBI operations, a recon- 
structed analog signal is passed to the SSI for verification of 
the calibration tones in the recorded receiver data. The DSS 
Radio Science Subsystem is configured and controlled by the 
DSS Monitor and Control Subsystem, while DRS status is 
routed to both the DMC and NRS. Figure 3 presents a 
functional block diagram of the ODA while performing VLBI 
Block 1 Phase 2 operations. 

3. VLBI Block 2. Natural radio source signals are received 
by the Wideband Multi-mission Receiver and pass through 
VLBI Intermediate Frequency Converters. The ODA config- 
ures both the VLBI IF Converters (Bandwidth selection) and 
the Wideband MMR synthesizer (fixed frequency informa- 
tion). From the VLBI IF converters, the signals are received by 
the VLBI Converter Subassembly, where they are multiplexed 
and digitized. The Digital Recording Assembly receives the 
data from the VLBI Converter Subassembly and records it on 
magnetic tape. The DRA is controlled by the ODA; for 
example, recording is halted during slew between radio 
sources. Concurrently, the ODA receives and records calibra- 
tion and ancillary data from the DMC. These data are 
formatted for HSDL transmission to the NRS via the CMF, 
During VLBI operations, a reconstructed analog signal is 
passed to the SSI for verification of the calibration tones in 
the recorded receiver data. The DRS is configured and 
controlled by the DMC, while DRS status is routed to both the 
DMC and NRS. Post pass, the DRA and ODA recorded tapes 
are shipped via Network Information Control (NIC) to the 
Network (NWK) VLBI Processor Subsy stem . Figure 4 presents 
a functional block diagram of the DRA while performing 
VLBI Block 2 operations. 

III. Functional Requirements of the DSS „ 
Radio Science Subsystem 

A. Functional Requirements for Block 1 VLBI (ODA) 

Block 1 VLBI system requirements on the ODA are as 
follows: 

, (1) 250 kHz bandwidth 

(2) 1 bit/sample 

(3) 500 kbit/s data rate 

(4) 10^ bits data volume 

(5) 1X10“^ bit error rate 


(6) 2-nanosecond maximum sampling jitter of any bit with 
respect to station reference 

(7) 10-microsecond time tag accuracy 

B. Functional Requirements for Block 2 VLBI (DRA) 

Block 2 VLBI system requirements on the DRA are as 
follows: 

(1) 8 parallel channels 

(2) 2 MHz bandv/idth/channel 

(3) 1 bit/sample 

(4) 32 Mbit/s data rate 

(5) 5 X 10^^ bits data volume 

(6) 1X10“^ bit error rate 

(7) 2-nanosecond maximum sampling jitter of any bit with 
respect to station reference 

(8) 5-microsecond time tag accuracy 

(9) Control of the DRA by the ODA 

IV. DRS Planned Implementation Schedule 

The planned implementation dates for the various DRS 
capabilities are presented below: 

VLBI Block 1 Phase 1 

(Includes VLBI Converter Subassembly, modification of 
ODA hardware and software, DTK interface, wideband 
interface.) 

DSS 14 July 1, 1979 

DSS 43 July 1,1979 

DSS 63 July 1, 1979 

VLBI Block 1 Phase 2 

(Includes MMR interface and ODA software modifications.) 

DSS 14 July 1, 1980 

DSS 43 July 1, 1980 

DSS 63 July 1,1980 

VLBI Block 2 

(Includes modification of DRA, ODA, VLBI Converter 
Subassembly, and DMC interface.) 

DSS 14 April l,i981 

DSS 43 April 1,1981 

DSS 63 : April 1,1981 
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Absolute Flux Density Calibrations: 
Receiver Saturation Effects 

A. J. Freiley, J. E. Ohlson^ and B. L. Seidel 

Radio Frequency and Microwave Subsystems Section 


The effect of receiver saturation is examined for a total power radiometer which uses 
an ambient load for calibration. Extension to- other calibration schemes is indicated. The 
analysis shows that a monotonic receiver saturation characteristic could cause either 
positive or negative measurement errors, with polarity depending upon operating 
conditions, A realistic model of the receiver is made using a linear-cubic voltage transfer 
characteristic, Th^yaluation of measurement error for this model then provides a means 
for correcting radio source measurements. It also provides the means for assuring that this 
source of error is small in a particular situation. 


i. Introduction 

The Jet Propulsion Laboratory’s Antenna Gain Calibration 
Program has demonstrated the feasibility of accurately deter- 
mining the gain performance of large aperture antennas using 
radio metric measurements of natural radio sources (Ref. 1). 
One phase of this program examined systematic errors and 
their effects on the radio metric measurements. The measure- 
ment of low system operating temperature and ambient noise 
reference temperature demands operation of the radiometer 
receiving system over a fairly large dynamic range. This 
necessitated examination of receiver saturation. The saturation 
may be small, but it can have impact on radio source 
temperature measurements at the 1% error level. 


Consider a classical total power radiometer which goes 
through the sequence shown in Table 1. The parameters shown 
are defined as 

T - System operating noise temperature (cold 

s^). 

T^ = Increase in system operating noise tempera- 
ture due to radio source. 

^AMB “ System operating* noise temperature when 
terminated in an ambient load. 

G = System power gain through the receiver. 
This value is assumed to only be known 
nominally prior to use of ^ ^ 

reference. 


^Dr. Ohlson, a consultant to Section 333, is a Professor at the U.S. 
Naval Postgraduate School, Monterey, CA. 


P^ = Output powers of an ideal linear receiver as 
indicated by an ideal detector. 
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Observed Output powers of a real receiver as From ( 1 ) with P 3 substituted for the measurement of 

indicated by an ideal detector. gain for the real receiver is 


First, consider an ideal linear receiver. Step 1 (Table 1) 
gives us when the antenna is pointed. slightly away from 
the source. Step 2 gives us G(Xo^ + when on source. 
Step 3 gives us a calibration value of The value of 

^AMB is known from physical temperature measurement of 
the load and other RF hardware calculations. Using a prime to 
denote estimated values, the estimated gain factor is 


G'=—^ = C^G ( 6 ) 

Then from ( 2 ) and ( 6 ), with G ^ we 

have the measurement of source ternperature as 





T 

AMB 


( 1 ) 




r = 

s 


G' 



(7) 


However, since we are starting by assuming a linear receiver, 
we have P^ = GTj^y^^ and thus G* - G. Then the estimated 
value of is measured in the classical manner as 


t 

s 


( 2 ) 


The error in measurement is 


Ar = 


s s 



( 8 ) 


and we clearly obtain 


Observe that C 2 and C 3 are specific values of the 
function C(GTj^) where Ti^ is total system temperature, i.e., 


li. Real Receiver Analysis 

Let the real receiver be described by the compression 
factors Cj, C 2 and C 3 which have values between 0 and 1 . 
These represent the reduction of receiver output due to 
saturation so that (see Table 1) 


Cx = c(Gr^) 

Cj = c[G(r^p + 7;)] (9) 

^3 = 


p[ = qp, <p^ 

(3) 

K = <P^ 

(4) 

K = ^3^3 <^3 

(5) 


We assume that the saturation is not in the detector, but in the 
intermediate frequency (IF) amplifier prior to the detector. To 
insure that this is the case, most accurate radiometers use a 
very accurate adjustable calibrated IF attenuator prior to the 
detector. The detector is used simply as an indicator and the 
IF attenuator is adjusted so the indicator is returned to a 
standard reading. This technique eliminates nonlinearity error 
in the detector because the detector always, operates at the 
same level for each measurement step. The data is then the 
difference of the attenuation readings from the attenuator 
dial, for Steps 1, 2 and 3 in Table 1. 


We now define a monotonic compression characteristic for a 
receiver as one for which C(x) is monotonically decreasing 
with:x:. Thus, if for example 


op op s AMB 


( 10 ) 


then if a receiver had monotonic compression, it would have 


\> C^> C^> C^>Q ( 11 ) 

i.e., has less compression and is closer to unity. 

The first conclusion will now be made. The error in ( 8 ) can 
be positive or negative. To show tHis consider two special 
cases, satisfying ( 10 ) and ( 1 1 ). 

(i) Let T^p + 7 ; = Then = C 3 and from ( 8 ) AT 

is negative since Cj > - 
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( 2 ) Assume compression occurs only for ?\mb* Then 1 = 
Cj = > C 3 > 0 and from ( 8 ) AT is positive. 

For any radiometer we will usually have < ^AMB- 
Assuming is less than ^AMB> two operating conditions^ 
exist. One condition is in (10). It and the other one possible 
are tabulated in Table 2 . The error of the first was shown to 
possibly be positive or negative. For the second in Table 2 , the 
error is always negative, under the monotonic assumption in 
(11). The breakpoint between positive and negative errors is 
difficult to find by simply setting. ( 8 ) equal to zero, since we 
have not yet developed an analytic function for (9). The main 
result of this section is in ( 8 ). This is useful when the 
compression factors are known because we can then calculate 
the measurement error AT. 


An amplifier is usually characterized for sine wave inputs 
and outputs. Let - A cos ojt, where to is the IF center 
frequency. Substitution of this in ( 12 ), and manipulation 
gives 



(13) 


The second term is rejected by the bandpass filter so the 
envelope voltage amplitude of the output sine wave is 


e 


^ A- 


m! 

4 


(14) 


ill. Modeling of Compression 

The result above is useful only when and C 3 are 

each known. This is usually not the case. Usually a radi- 
ometer user knows only that he is operating some amount 
(backoff) below a reference compression point, e.g., the 
- 0.1 dB point. In this section we consider a linear-cubic 
model and can calculate the measurement error when ^ven 
only the backoff value. 

The voltage input-output relationship of an ampHfier is 
often modeled with the linear-cubic mo del 


We now choose a standard reference compression point, C^. 
Commonly, this is expressed in decibels so (in dB) = 
10 logjQ C^, e.g., “0.1 dB from = 0.97724. Since 
is a power ratio, then is the factor by which e in (14) 
is reduced below the amplitude A obtained with a linear 
receiver. We wish to find the input envelope A = where 
this compression occurs. Thus 


A 


o 



= AC^!^ 

o r 


(15) 


%ut 



and hence 

( 12 ) 


This model is widely used for intermodulation product cal- 
culation, and seems appropriate here as well. We assume a 
bandpass filter to foUow the amplifier prior to the detector. 
Tlius, second and fourth order terms could be included in 
( 12 ) but their output would be rejected by the filter and 
have no effect. The gain is unity for small Saturation 
mcreases for larger We assume that we wfil operate far 
below the point at which ( 12 ) has its maximum value and 
starts decreasing for very large 


< = fd - 

Since power comes from the envelope as ^^/2, we have the 
input power giving compression for a sine wave is 

jp, = |(1 - Cp) . (17) 


^ < 

A more general form for a real amplifier which incorporates a gain 

parameter a and a saturation parameter ^ is 
”out = “•'in - 

With no loss of generality m what is to follow, a and p can be taken 
to be unity as in (12). 


We now examine compression for a noise input. For 
Gaussian noise the envelope A has the Rayleigli probability 
density 

p(A) = ^ exp A>0 (18) 
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where = average input power of the noise. The output 
envelope after the filter is in (14) so the output power 
(observed by the detector) is 


P 


o 




j p{A) dA 


=p-eFi+^i 


(19) 


The compression factor for noise is implicitly defined as 


C =-^=l- 6 P+^p^ 

n P n 2 

n 


n 


( 20 ) 


Now let the input noise power P^ be below the sine wave 
compression point P^ by the backoff B<1: 


Solving for B gives B = 1 / 2 . Thus we conclude that a noise 
input must be backed off a factor of 2(= 3 dB) from a sinu- 
soidal input to have the same compression. This is reasonable 
since the noise input has large excursions which are com- 
pressed more than a sinusoidal input. 


IV. Application to Radiometer Error 

In ( 8 ) we have the error in measurement of in terms 
of the compression factors C^, C 2 and C 3 . In (23) we now 
have an expression for these factors under the assumption of 
a linear-cubic receiver model. 

We first characterize the receiver for a sinusoidal input by 
establishing the reference level at which a compression of Cy 
occurs, say -O.! dB. This reference level can be determined 
by use of a signal generator and precision attenuators. For 
use as a radiometer, with a noise input, this same com- 
pression occurs for a reference temperature T^. Tliis tempera- 
ture represents an mput noise power reduced by a backoff 
of 1/2 with respect to the reference sinusoid. Thus for total 
input temperature backoff is given by 


P ^ BP (21) 

ns ^ ^ 

We then have, by ( 21 ), ( 20 ) and (17) 

- I - 4^(1 - Cy^) + 6B\\ - cy^)^ (22) 


We have already presumed smalL saturation, so is very 
close to unity. Thus (1 ^ C^^^) is «1 and since ^ < 1 , the 
third term in (22) is second order in B{\ - C^^^) and is 
negligible with respect to the second term. Thus 

“ 1 - 45(1 - (23) 

This is the main result of this section. 

An important question may now be answered. At some 
reference level the compression for a sine wave input is C^. 
At what backoff from this reference level will a noise input 
have the same compression? We start with (23) and set = 
Cy-C, Since ^(1 - C^/^) «1, we use (1 +;c)W 2 \ +^/2 

to get from (23) 

- 1 - 25(1 (24) 


. = W, 

From (23) we thus have 

c„ = 1- 2^<l-cy") (26) 

ly 

This gives the compression vs input temperature under 
the condition that a compression of Cy occurs at reference 
temperature Ty. Square-rooting (26) and using the argument 
prior to (24) easily shows that when T.^ = 

For simplicity, define 

5 = 2 (1 - C\I^)IT^ (27) 

Then the compression factors needed for (8) are of the form 
in (9): 

^2 = 1 - 3 ;) 

^3 = 1 - «^AMB (28) 
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Substitution into (8) gives 


AT _ 


V. .Example 

To illustrate tlie above calculations, consider tlie following 
example parameters representative of a liigli performance 
(29) radiometer using a maser' 


Since (1 - will be small, and assuming 

much larger than we can drop the second term in the 

denominator bracket to give 


== 15 K 

op 

T = 100 K 


- 2(1 - C;/2) - 

s ^ r 


(30) 


This is the main result of tills section. In this form it can be 
used to calculate correction factors for measurements. 


^amb*300K 


(35) 


Let the receiver be set so that a temperature = 400 K gives a 
compression of -0,1 dB (C^= 0.97724). The error is, from 
(30), at the i% level: 


Several observations can be made; 

(1) The fractional error AT/T^, linearly decreases as T^ 
increases. 


AT 

- 0.0097 


(36) 


(2) The error goes to zero when 


T - T - 2T 

^AMB ^op 


The maximum error for any T^ satisfying (33) with 
<270, is, from (34) 


(31) 


(3) The error is positive or negative depending upon 
whether T^ is below or above the value in (31). 


< 0.0172 

5 


A special case is of considerable value. Let and be 
very small compared to ^am b- Then 

~ = 2(1 - (32) 

■*5 ^ r 


Further, as long as 


VI. System Constraint for Small Error 

Let us now find the constraint on the receiver com- 
pression so that the system error will be small. Assume that 
(33) holds and let us use the bound in (34) so that our 
constraint will hold for any value of T^ in (33), Let us take 
a value of AT/T^- ^MAX worst case error we can 

accept. Then for equality in (34) we have tlie tradeoff in 




(33) 


P 

MAX 


2(1 - C\!^) 


^AMB 


(37) 


the result in (30) is bounded by (32). Thus, we obtain the 
useful result which can be used to provide an upper bound 
on saturation effects (condition 33 must apply): 



2(1 - cy*) 


AMB 


(34) 


Examples of this tradeoff are shown in Table 3 for a value, 
of ^max “ 0-003. Tliis Table shows that -0.026 dB com- 
pression at a value of T, which is 3 dB above guar- 

antees tliat ATfT^ is below 0.003. Tliis is also true if there is 
-0.131 dB compression when is 10 dB above 
Clearly there is an infinite number of pairs of values satis- 
fying (34). The usefulness of (34) is that only one value of 
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compression need be specified and this can be done for a 
nearly arbitrary value of 

In many radiometers, the system calibration is done by 
use of a calibrated active noise source reference (gas dis- 
charge tube or noise diode) which adds a precise ^ount of 
noise tO’ the system input This is as opposed to the use of 
an ambient load for calibration as done above. With small 
changes, it is straightforward to apply the technique above 
to the case of an active reference. The key point is that 
Step 3 in Table 1 now has ^‘off source-noise reference on” 
and temperature Top + Te, where is the calibrated noise 
reference contribution. Also the gain estimate in (6) becomes 



The remainder of the analysis closely follows the above and 
is omitted for brevity. 


VII. Conclusions 

An analysis has been made of measurement error in a 
radiometer due to receiver saturation. The general result is in 
(8). A linear-cubic model of receiver saturation was then 
assumed and an explicit result was obtained in (30). 

The measurement error is characterized by only one mea- 
surement of compression which' can be made at nearly any 
reference temperature. The work here is useful for calculat- 
ing correction factors, or at the very least, verifying that 
error is insignificant for a particular situation. 
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Table 1. Summary of measurement 


Step 

Description 

Temperature 

Output power of 
ideal receiver 

Output power of 
real receiver 

1 

Off source 
(cold sky) 

T 

op 



2 

On source 


i>2 = G(rop + Tp 


3 

Ambient load 

^AMB 

•^3 ^ 



Table 2. Error polarity 


Condition 

Error, AT 

T +T <T 
op 

Positive/negative 

T + T > T 

^op 

Negative 


Table 3. Saturation tradeoff 


Case 


Worst case compression 

A 

0.5 (-3.0 dB) 

0.99401 (-0.026 dB) . 

B 

0.1 (-10.0 dB) 

0.97023 (-0.131 dB) 
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DSN Portable Zero Delay Assembly 

E J. Serhal, Jr., and T, Y. Otoshi 
Radio Frequency and Microwave Subsystems Section 


This article presents design and test data on portable zero delay assemblies that have 
recently been delivered to DSS 14, 43, and 63, These portable assemblies are field-use 
delay standards that will he used to periodically calibrate the Block IV Translator ranging 
paths at each 64-m antenna Deep Space Station. 

/ 

I. Introduction II. Description 


The Block IV Translator Method (Ref. 1) is currently being 
used to determine the ground station, range delays at all 64-m 
antenna Deep Space Stations (DSS)- The validity and accuracy 
of this method are dependent upon the Block IV Translator 
ranging paths being recalibrated (1) at periodic intervals (once 
or twice each year) or (2) whenever significant component or 
configuration changes have been made in the translator ranging 
path. 

The calibrations of the translator ranging path delays have 
been performed in the past years with an R & D Zero Delay 
Device (ZDD) (Ref. 2). The ZDD was shipped at various times 
to the three 64-m antenna DSS and has been used successfully 
to make translator delay measurements. Recently, the design 
features of the R & D ZDD have been improved and incorpo- 
rated into a field-worthy version called the Portable Zero 
Delay Assembly (PZDA). PZDAs have been fabricated for the 
DSN and a calibrated unit has now been supplied to DSS 14, 
DSS 43, and DSS 63. Operational test procedures for trans- 
lator measurements with the PZDA have been written by DSN 
Network Operations and are currently being verified. The 
purpose of this article is to present pertinent design informa- 
tion as well as calibration test data on the PZDA. 


Figures 1 and 2 respectively show the original R & D ZDD 
and the new PZDA. For the PZDA, pads and step attenuators 
were made to be internal components of the assembly to 
simplify test procedures for field use. In addition, the PZDA 
includes some broadband filters and isolators (see Figs. 3 and 
4) to suppress potential mismatch and higher order harmonic 
distortion errors. Suppressing tliese errors is especially impor- 
tant when calibrating the PZDA delays by the back-to-back 
method to be described later in this article. The tradeoff for a 
more accurate unit for field use is a longer delay through the 
PZDA because of the additional built-in components. AU com- 
ponents including filters and isolators have at least ±40 MHz 
bandwidths at their approximate uplink-downlink center fre- 
quencies of 2113 MHz, 2295 MHz, and 8415 MHz. 

Figure 5 shows the S- and X-band mixer subassembly for 
the PZDA. This subassembly is a standard mixer assembly 
being used in the Block IV Doppler Translator Assembly and is 
also very similar to the mixer assembly in the R & D ZDD. 

The test cables supplied witli the PZDA as shown in Fig. 2 
are phase-stabilized F282 cables made by the Flexco Co. in 
DenvUle, New Jersey. These cables were selected (instead of 
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the smaller diameter FI 82 cable used for the R&D ZDD) 
because of their lower loss ^ and superior mechanical durability 
with flexing, especially at the type N-to-cable transition joints. 
Three 7.62-m (25 ft) F282 cables are supplied with each 
PZDA for the S-X signal paths and two 3.05-m (10 ft) F282 
cables for the L.O. paths. The F282 cable has a delay charac- 
teristic of about 4:72 ns/m and attenuations of 0.66 dB/m 
in the 2.1 to 2,3-GHz range and 1.25 dB/m at 8.4 GHz. 


III. Test Method 

The PZDAs were calibrated by use of the back-to-back 
method wliicli is depicted in Fig. 6. The PZDA on the right- 
hand side is the unit under test while the unit on the left is the 
comparison standard. The comparison standard is also a PZDA 
but differs in tliat tlie S-X isolators are reversed to permit 
signal to flow in the reverse direction (J6 and J7 ports to J3 
port). When making S-S delay calibrations, the 182-MHz L.O. 
drive is turned on and simultaneously feeds port J1 of each 
unit. The uplink test signal (of approximately 2113 MHz) is 
fed into port 13 of the unit being tested. This uplink signal 
becomes up-converted to approximately 2295 MHz in the test 
unit and then down-converted back to 2113 MHz when it 
passes through the standard comparison unit. When measure- 
ments are to be made for S-X band delays, the 182-MHz L.O. 
drive is turned off and the 6.3-GHz L.O. drive is turned on to 
simultaneously feed port J5 of each unit. The 2113-MHz test 
signal becomes up-converted to approximately 8415 MHz in 
the test unit and down-converted back to 2113 MHz when 
passing through the standard comparison unit. 

Insertion phase versus frequency measurements are made 
for an input test frequency range of 2113.5 ±8.5 MHz using an 
HP 841 OA Network Analyzer and a microwave source driven 
by a frequency synthesizer. Group delay is calculated from the 
slope of the insertion phase versus frequency data. It is neces- 
sary to (1) measure the delay of the total insertion delay of 
the back-to-back assemblies as shown in Fig. 6 and (2) measure 
separately the delays of any test cables (including pads) that 
were used to interconnect the test unit to the standard unit in 
the back-to-back configuration. The delay of the unit under 
test is then determined by subtracting the delay of the test 
cables and delay of the standard comparison unit from the 
total measured insertion delay. If the delay of the standard 
comparison unit is not known, then the delay of the unit 
under test can be determined by subtracting out the test cable 
delays from the total measured insertion delay and then divid- 
ing the net result by two. This latter method is valid only if 
both the test and standard units were built to the same 
specifications and are therefore electrically identical to within 
an acceptable tolerance. 


For the test results of the article, the delay of the standard 
comparison PZDA was known and determined previously by 
connecting it to an improved R & D ZDD and calibrating it in 
the back-to-back configuration. The R&D ZDD (improved 
version) in this configuration was the basic primary standard 
and was used because it has S-S and S-X delays of 2.3 ns in 
each path which are known to accuracies of ±0.2 ns (la). 
These accuracies were established previously from extensive 
testing and cross-comparisons. 

After completion of the PZDA calibrations, the two 3-dB ■ 
pads in the test unit as shown in Fig. 6 were replaced by two 
20-dB pads (of nearly the same physical lengths and electrical 
delays) to restore the PZDA to its operational configuration of 
Fig. 3. The lower value (6 dB) padding was temporarily needed 
during the calibration because a relatively strong test signal 
level is required for the Network Analyzer insertion phase 
measurements. 


IV. Test Results 

Table 1 shows the final calibrated delays of the PZDAs that 
were sent to various DSSs. The applicable serial numbers and 
DSS identification are also tabulated. The accuracy of the 
delays are estimated to be better than ±1.0 ns (la). Due to the 
fact that all components in the PZDA including isolators and 
filters are wide-band (see Section II of this article), the delay 
calibrations in Table 1 are applicable over a wide-band and 
apply to all S-S and S-X frequencies for DSN channels 5 
througli 27 listed in DSN documents. 

The insertion losses of the PZDA in the final operational 
configuration were determined to be 54 ±2 dB for the (S-S) 
path going from port J3 to J6 and 55 ±2 dB for the (S-X) path 
going from port J3 to J7. The high insertion losses are mainly 
due to the two 20-dB pads (see Fig. 3) tliat are needed to 
attenuate an anticipated test signal level of ±19 dBm to a level 
of about -20 dBm so that internal S- and X-band mixers will 
be operated in* their linear regions. The strong test signal level 
of +19 dBm is derived at each 64-m DSS from the output of a 
54-dB waveguide coupler which samples the 20-kW transmitter 
power. 


V. Conclusion 

A PZDA that is field-worthy, operationally easy to use, and 
accurately calibrated has been developed. A PZDA has been 
shipped to DSS 14, DSS 43, and DSS 63. Preliminary results 
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of Translator delay measurements from DSS 43 and 63 indi- 
cate the- results obtained with the PZDA agree reasonably 
well with test data obtained in previous years with the R & D 
ZDD. 


The standard comparison PZDA as well as calibration pro- 
cedure documentation have been sent to the DSN Maintenance 
Center (DMC). Recalibrations of the PZDA in the future can 
be done by shipping the units to DMC. 
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DSS 


Table 1. Calibrated PZDA delays 


SN 


s-s 

(ns) 


s-x 

(ns) 


14 

43 

63 


2 

3 

4 


18.8 

18.1 

18.2 


14.3 

14.2 

13.9 



Ffg. 1. R&D Portable Zero Delay Device with external test cables and attenuators 



Fig. 2. DSN Portable Zero Delay Assembly 
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Rg. 4. Microwave components mounted on underneath side of front panel of 
Portable Zero Delay Assembly 



Fig. 5. S- and X-band mixer subassembly of Portable Zero Delay Assembly. Same 
as Translator Assembly 9456791-3 
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Fig. 6. Back-to>back test configuration 
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Cost Standard Practice 
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TDA Planning 


The DSN is developing a Life-Cycle Cost Standard Practice, This report compares the 
developing Practice to those of industry and the Department of Defense, Results show 
that the DSN uses the accepted concept of Life-Cycle Costing, tailoring the concept to 
DSN specific needs, but does not push the concept past the point of prevailing theory. 


I. Introduction 

The DSN is developing a Life-Cycle Cost (LCC) Standard 
Practice, At this stage of development, it is useful to compare 
the Practice to those of industry and the Department of 
Defense (DOD). Specifically, three questions are posed for 
study: 

(1) How widely used is the life-Cycle Cost concept? 

(2) Is the DSN LCC Practice, at this stage of development, 
in accordance with generally accepted standards? 

(3) How does the DSN LCC Practice compare with those 
of industry and the military; Le,, is the DSN at the 
state-of-the-art? 

These questions are addressed in the following report. 

II. Use of Life-Cycle Costing 

A. Department of Defense 

LCC is clearly used by the Department of Defense (DoD), 
An investigation of the LCC concept was made in the early 


1960’s by the Logistics Management Institute for the Assistant 
Secretary of Defense for Installations and Logistics (Ref. 1). 
Since that time, numerous DoD and service directives have 
clarified and ordered the use of Life-Cycle Cost analysis. An 
example is the DoD Directive 5000.28 which states: 

life-cycle cost estimates will be used as a basis for cost 
tradc-off analyses when considering acquisition versus 
O&S costs, comparing competing prototypes or 
comparing current versus new systems” (Ref. 2). 

Furthermore, in using LCC, 

‘‘Life-cycle cost objectives shall be established for each 
acquisition and separated into cost elements within the 
broad categories of development, production, operation, 
and support. As system definition continues, the cost 
elements are firmed into cost goals to which the system 
will be designed and its cost controlled” (Ref. 3). 

Thus, the Department of Defense uses Life-Cycle Costing to 
select a system, and then uses ‘‘design to Life-Cycle Cost 
elements” to control costs. One source reports that use of 
Life-Cycle Cost .analysis has resulted in “significant cost 
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savings,” validated by the Air Force Auditor GeneraFs Office 
(Ref. 4), 

B. Military-Related industries 

A sampling of military-related industries was conducted by 
the author to determine the extent to which LCC is used, 
Martin Marietta, Northrop, and Hughes Aircraft were con- 
tacted, and all use some form of Life-Cycle Cost analysis. One 
of the three has a documented standard LCC approach 
(Ref. 5), two have centralized LCC departments, and all three 
use approaches which are tailored to specific use. 

Numerous industrial articles have been written on LCC 
analysis, and will be referenced throughout the remainder of 
this report. The general conclusion is that mihtary-related 
industries do use some form of LCC analysis, but. no single 
approach is standard. 

C. Other Industries 

Other industries use the LCC technique for purposes such 
as comparing alternatives: 

“In comparing alternative solutions to accomplish a 
particular goal, the system showing the lowest life-cycle 
cost will usually be the first choice, assuming the 
performance requirements defined by the goal are 
adequately met” (Ref, 6). 

Industrial articles on LCC analysis were found in the following 
magazines and journals: 

Industrialization Forum 

Logistics SPECTRUM 

Professional En^neer 

Proceedings Annual Reliability and Maintainability Sym- 
posium 

Management Accounting 

The articles will be referenced, as appropriate, in the following 
pages. The same conclusion is reached through review of these 
articles; the use of the Life-Cycle Cost concept is widespread, 
but the specific approaches and models used are not stan- 
dardized. 

HI. Accepted Life-Cycle Cost Standards 

A. Definition of “Life-Cycle Cost” 

General agreement exists on the meaning of “Life-Cycle 
Cost.” The Department of Defense Directive 5000,28 states: 


“The Life Cycle Cost of a system is the total cost to the 
Government of acquisition and ownership of that system 
over its full life. It includes the cost of development, 
acquisition, operation, support, and where applicable, 
-disposal” (Ref. 7). 

Similarly, a student of the Air Force Institute of Technology 
writes that Life-Cycle Cost is: 

. a method for evaluating and considering the total 
costs of ownership, including the costs of acquiring, 
operating, supporting, and disposing of an item, less its 
residual or scrap value, during its useful life” (Ref. 8). 

A Maitm Marietta representative states: 

“Life-Cycle Costing (LCC) is a Department of Defense 
(DoD) acquisition or procurement technique that con- 
siders operating, maintenance, and other costs of owner- 
ship, as well as acquisition price, in the award of 
contracts for hardware and related support. The object 
of this technique is to insure that the hardware procured 
will result in the lowest overall ownership cost to the 
Government during the life of the hardware” (Ref. 9). 

A similar concept is referenced m Professional Engineer: 

“Life-cycle costing is an analysis of the total cost of a 
system, machine, device, building, and the like over its 
anticipated useful life” (Ref. 10). 

Thus the DSN definition — Life-Cycle Cost is implementa- 
tion cost, plus M&O costs over the life of the capability, plus 
deimplementation cost — is in accordance with the generally 
accepted definition. 

B. Life^Cycle Cost Uses and Purposes 

Although the definition of Life-Cycle Costing is standard, 
its uses vary with the background of those employing it. 
Donald Earles of Raytheon Company explains it: 

. over the past ten years, life-Cycle Costing has 
evolved into four things; a costing discipline, a procure- 
ment technique, an acquisition consideration, and a 
design trade-off tool.” 

The specific use depends on the user. The DSN employs LCC 
analysis as an acquisition consideration and a trade-off tool 
(Ref. 1). 

C. Costing Philosophies 

Just as LCC uses vary among users, costing philosophies 
vary with the specific application. Costs of a system may be 
full costs, delta costs projected from some baseline, or even 
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variable costs For example, a RAND weapon case study 
remarks : 

is important to note that this definition of weapon 
system LCCs is intended to capture only the variable 
costs associated with the weapon system. We make no 
attempt to reflect any of the ‘"fixed” costs which are 
necessary for the management of the organization (e.g., 
HQ/US AF) which uses the system, but which are not 
incurred by or for the weapon” (Ref, 1 1). 

This is the costing philosophy assumed for the DSN pathfinder 
LCC analysis (Ref. 12). Other companies use different ap- 
proaches. Hughes Aircraft Corporation, for example, uses 
absolute costs, wlule Northrop employs absolute or delta 
costs, as the situation demands. 

D. Project Life-Cycle Period 

Another Life-Cycle Cost concept which varies in use is that 
of project life-cycle period. An Air Force pubheation notes 
that the LCC period may be 10*20 years for a DoD weapon 
system (Ref. 13). A specific cost analysis done by the U.S. 
Army Electronics, Communications, and Material Commands 
uses a 10-year period (Ref. 14). This variability is summed up 
in Professional Engineer: 

“The Life-Cycle of the system or machine must be 
realistically appraised. It can be defined as the expected 
physical life or the estimated period to obsolescence, 
whichever occurs first” (Ref. 10). 

The variability of the life-cycle period is recognized hy ^ the 
DSN. 

E. Life-Cycle Cost Definitions 

Included in the DSN LCC Standard Practice is the 
requirement of definitions. Careful definitions msure fair 
comparisons and allow validation that no cost has been 
forgotten (Ref. 1 5). Although this is a generally accepted 
standard in government and industry, it is not practiced as 
thoroughly as is wise. Nevertheless, examples of LCC reports 
accompanied by careful definitions can be found. Two 
examples are A Summary and Analysis of Selected Life-Cycle 
<Costing Techniques and Models of the Air Force Institute of 
Technology (Ref. 16), and Life-Cycle Costing Emphasizing 
Energy Conservation, Guidelines for Investment Analysis of 
ERDA(Ref. 17). 

F. Cost Estimation Procedure 

It is clear from even a cursory examination of the 
Life-Cycle Cost literature, that there are a number of cost- 
estimating methods. Some of the titles of these methods are* 


per unit catalog price, cost-to-cost estimating relationships, 
non-cost-to-cost estimating relationships, specific analogy, 
expert opinion, simulation (Ref. 18), accounting, and para- 
metric. It is the last two which express the organizational 
nature of cost estimating and which concern us here. 

In general, the DSN uses the accounting approach, which 
lists all cost elements for a project over its life and, with a 
specific procedure, sums the costs. This method is explained: 

“The most frequently referenced of the life-cycle cost 
models in - use today are the additive accounting 
models.” They are useful for projects characterized by 
dollars and provide managerial visibility, companson of 
costs, identification of specific differences in bidders’ 
proposals, and visibility of hi^ support costs (Ref. 19). 

In comparison, the parametric approach is useful for esti- 
mating effects of learning curves, provisioning policies, 
management factors, repair and replacement policies, and 
system characteristics (Ref. 20). Only a few parametric models 
are sufficiently tested for DSN use. 

There are numerous cost models which mcorporate' these 
methods. For example, an Air Force Institute of Technology 
study lists twenty such models (Ref. 21). None of these can be 
used as a standard cost model. As a Honeywell employee 
states: 

“The state-of-the-art of cost modeling, however, is not 
to the point where existing cost models can be used 
universally in whole or in part for wide varieties of 
applications. They are highly tailored in terms of 
systems treated, scope of elements covered, and cate- 
gories of costs included” (Ref. 22). 

Thus, there are generally accepted cost estimation proce- 
dures, one of which is used by the DSN. There are no models 
incorporating these methods which can be used as a standard 
LCC model. 

G. Cost Elements 

Within the accounting cost estimation method are the cost 
elements themselves. Again, there is no single list of detailed 
cost elements applicable to all projects. Examples of cost 
element lists are given in Tables 1 and 2 and more can be 
found in the literature (Refs. 13, 20, 23). At this time, 
however, DSN cost elements must be determined for each cost 
study by those familiar with the project. 

H. Monetary Adjustment 

Another part of Life -Cycle Costing is the monetary 
adjustment of projected costs. Review of the literature shows a 
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lack of standardization here also. The Office of Management 
and Budget suggests the use of comparative inflation and a 
general discount rate of 10% (Ref. 24). A RAND study, 
however, uses a discount rate of 5% (Ref. 25). Northrop 
prefers to work solely in constant dollars and inflate only if 
required for bidding purposes. The Martin Marietta LCC 
process allows total discretion in the inflating and discounting 
procedures used on each cost estimation. 

Thus, the monetary adjustment used in LCC estimation will 
depend upon the project, the intended use of tlxe estimation, 
and any applicable laws. 

I. Life*Cycle Cost Presentation 

The form in which a LCC estimate is presented is 
dependent upon the project and the purpose of the estimate. 
The DSN pathfinder LCC analysis used a reporting format 
designed to reflect the specific cost elements and algorithm 
used in that project. Other formats can be found in the 
literature (Ref. 26). An example is in Table 3. 

J. Standard Practice 

As stated previously, the DSN is developing a Life-Cycle 
Cost Standard Practice. The question arises: Are LCC standard 
practices common in industry and government? The answer is 
that in reviewing sixty sources, the author encountered only 
four sets of documents which could be considered as standard 
practices. ERDA published the Life-Cycle Costing Empha- 
sizing Energy Conservation, Guidelines for Investment Analy- 
sis (Ref. 17), and Martin Marietta has the Computerized 
Costing Methodology which incor porate s Life-Cycle Costing 
(Ref. 5). The Navy has an Economic Analysis Handbook 
(Ref. 27), and the DoD has a series of documents published as 
guidelines for various weapon procurements (Ref. 28). 

The need for standard practices, however, is well docu- 
mented. A 1975 Air Force working group found that: 

“At least five actions are required in order to establish 

an adequately effective life cycle cost analysis capability 

in program offices.” 

“1. Program offices must be provided with a source of 
personnel familiar with analytical techniques.” 

“2. These engineers and analysts must be given general 
guidance on how to develop, adapt and use life 
cycle cost models for specific applications.” 

“3. Program office and supporting personnel should 
have access to a short course in the subject of 


development and application of life cycle cost 
models and methods.” 

“4, Periodic life cycle cost methods workshops should 
be held ” 

“5. Finally, program office personnel should be pro- 
vided with a central focus of expertise where lessons 
learned in each new life cycle application are 
integrated with existing life cycle cost models and 
methods . . (Ref. 29). 

Similarly, a RAND study states: 

“In order to advance the current state-of-the-art of LCC 
estimating, at least the following improvements” are 
needed: 

“1. Operational consistency of LCC procedures” 

“2. Better cost allocation rules” 

“3. A nomenclature directory” 

“4. Project LCC decision requirements” 

“5. An LCC data base” (Ref. 30) 

Thus, the need for standard practices is generally accepted, 
and some such documents already exist. The DSN, with 
others, is responding to this need. 


IV. General Conclusions 

Life-Cycle Costing is a popular technique, widely used in 
the military and industry. It is generally accepted that 
Life-Cycle Costing is a specific concept whose use should be 
accompanied by LCC definitions, a stated costing procedure, a 
stated project life-cycle period, monetary adjustment proce- 
dures suited to need, and at least a locally accepted LCC 
standard practice. The DSN is in accord with these generally 
accepted standards. The specific application of the standards is 
user and project dependent, as found in both the LCC 
literature and the DSN, 

The DSN is not at the state-of-the-art of Life-Cycle Costing. 
The state-of-the-art is an advanced and fluctuating position 
including such things as economic risk analysis, performance 
effectiveness measurement, multiple regression analysis, and 
unproven parametric models. These advances are not being 
ignored by the DSN, but rather will be incorporated as 
appropriate with the DSN environment. The DSN is using an 
accepted tool, adapted to local needs, to reduce project costs. 
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Table 1. Life-cycle cost categories 


Table 2. Weapon system life-cycle cost categories^ 


Major categories 

Suhcategories 

Research & Development 

Engineering 
Prototype Fabrication 
System Test & Evaluation 
Data 

Total Systems Management 
Training 

Producibility Engineering & Planning 

Investment, Nonrecurring 

Initial Production Facilities Effort 

Data 

Tooling 

Instruction Training 
Total Systems Management 

Investment, Recurring 

Hardware Cost 
Attrition 

Initial Spares & Repair Parts 
Support Equipment 
Au^dliary Equipment 
First Destination Transportation 
Initial Training 
Engineering Changes 
System Test & Evaluation 
Total Systems Management 

Operating 

Personnel 

Personnel, Replacement Training 
Consumption 

Integrated Logistics Support 
Depot Overhaul 
Transportation 

^Source: Thomas W Otto, Jr,, Life^Cycle Cost Model, U.S. Army 
Electronics Command Report ECOM-4338;JuIy 1975, p. 8. 


Major categories 

Subcategories 

RDT&E 

Research & Development 
’ Test Evaluation (RDT&E) 

Procurement 

Aircraft (Flyaway Cost) 

Support Equipment 

Training Equipment 

Technical Manuals 

SPO*Personnel 

Test & Evaluation (lOT&E) 

Ownership 

Spares: Initial, Replenishment, Engines, 
' WRM (HDWR) 

Maintenance: On-Equipment 
Maintenance: Off-Equipment 
Management Personnel: System/Ttem 
Managers 

Training. Tech Training 
Operations: Crew, CMD PERS 


— Military and Civilian 

— Base Only 


Base Operating Support 
Fuel 

Modifications (HDWR) 

Testing & Evaluation (OT&E) 

Item Transportation (2nd DSTNTN) 
Munitions and Missiles - Training 
Personnel PCS 
Attrition 
New Facilities 

Salvage & Disposal 

Disposal 

Modifications 

^Source: Marco R. Fiorello, Getting *‘Real*’ Data for Life-Cycle 
Costing, RAND Corporation, No. AD-AOlO 960. Jan. 1975, p.7. 


Table 3. LCC trade-off analysis of two competitive products for 
two years ^ 


Cost 

Vendor A 

Vendor B 

Product Price 

$200,000 

$170,000 

Installation 

3,000 

4,000 

Manning Labor (2 Yrs) 

46,720 

93,440 

Preventive Maintenance (2 Yrs) 

912 

1,632 

Corrective Maintenance (2 Yrs) 

2,800 

9,344 

Power Requirements (@ 0-025 /kW) 

1,168 

1,314 

Parts & Supplies Cost (@ 1% & 2% 



of Product Price Respectively) 

2,000 . 

3,400 

TOTAL 

$256,600 

S283.130 

^Source: Gary White and Phillip Ostwald,/* Life-cycle Costing,’' 
Management Accounting, Vol. 57, No. 7, Jan. 1976, p. 42. 
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A Brief Historical Introduction to 
Very Long Baseline Interferometry 

B. Benjauthrit 
TDA Engineering Office 


This article provides a short historical account of Very Long Baseline Interferometry, 
including the rationale, development, and experiments. 


I. Introduction 

The technique of Very Long Baseline Interferometry 
(VLBI) has been employed to achieve extremely high angular 
resolution in the study of radio sources. Unprecedented 
accuracies may be obtained for values of various geophysical 
quantities. As an example, it is now possible .to compare time 
scales on separate continents with an uncertainty of only one 
nanosecond (i0“^ sec). The main purpose of this article is to 
provide a short history of VLBI starting from the Michelson- 
Pease experiment of 1920 up to the present time. This 
includes the rationale, development, and various experiments 
behind the subject of VLBI. However, due to a large volume of 
references^ on the subject, only certain aspects are covered 
here. 

II. The Rationale Behind Radio 
Interferometry 

Until recently images of celestial objects obtained from 
radio telescopes lacked the detail of those obtained from 
optical telescopes. The reason is that the resolution of a 
telescope increases with the ratio of its aperture to the 
wavelength of the received signal, and radio wavelengths are 


^ All references in this article refer to “An Extensive Bibliography on 
Long Baseline Interferometry” following this article. 


about a million times longer than light wavelengths. Thus, 
radio telescopes are fundamentally confined to a poor angular 
resolution when compared to optical telescopes. However, in 
practice this is not the case for two reasons (Ref. .164): 

(1) The resolution of large optical telescopes is confined 
not by their size but by irregularities in the earth’s 
atmosphere. The limit is about 1 arc second. This is 
about 100 times better than the unaided human eye. 
Also, due to longer wavelength, the fluctuation of the 
incoming signal path length through the atmosphere at 
radio frequenaes is much smaller than that at light 
frequencies. 

(2) The signal must be coherent, or in phase, over the 
entire dimensions of the telescope. Coherent radio 
waves are much easier to obtain than coherent light 
signals, so that radio telescopes can operate much 
closer to the theoretical resolution limit than optical 
telescopes. 


III. Early Need for Radio interferometry 

Althou^ resolution increases as wavelength decreases, the 
performance begins to deteriorate significantly when the 
wavelength approaches the dimensions of the structural 
imperfections in the antenna. Further, the largest antennas 
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available were paraboloidal type; they had the least precise 
surfaces. Thus, the best resolution that had been obtained with 
paraboloidal antennas did not depend strongly on wavelength 
and was about 1 arc minute. Although it is possible to build a 
more precise large antenna, it would not provide better 
resolution than 0.1 arc minute (Ref. 177). This is the reason 
why radio astronomers turned to interferometry, where in 
effect two relatively small antennas function as opposite edges 
of a single huge radio telescope. 

The development of interferometry began as early as 1 890 
when Albert A. Michelson published a paper (Ref. 237) 
describing a technique of modifying an optical telescope to 
receive star light along two paths, and determining the angular 
diameter of a star by means of interference of its own light 
waves. However, it was not until 1920 that his technique was 
successfully demonstrated on the 100-inch telescope of 
Mt. Wilson, in collaboration with F. G. Pease, Figure 1 shows a 
sketch of the Michelson-Pease stellar optical interferometer. 
Here, a li^t wave is intercepted by two separated mirrors 
which reflect two beams to a common point where the beams 
are combined. If the path of one beam is made slightly longer 
or shorter than that of the other beam, the light waves in one 
beam will be out of phase with the waves in the other beam. 
When the beams are combined, one sees a pattern of 
alternating light and dark ‘'fringes.” 

Attempts to extend the optical interferometer to longer 
baselines in order to achieve higher resolution have not been 
very successful, mainly due to the difficulty in obtaining 
coherent light waves from the system and in maintaining the 
alignment of the mirrors within a fraction of the wavelength. 
It is much simpler to solve this problem by using radio waves. 
Two or more antennas can be used to synthesize large 
apertures and achieve high angular resolution. With such an 
instrument, one can determine not only the size and shape of 
discrete radio sources but also their precise positions in the 
sky. 

IV. Early Development of Radio 
Interferometry 

Early radio interferometers were essentially radio analogs of 
the above well-known optical device. However, in radio 
interferometers, the signals are combined and compared (or 
cross-correlated) electrically. This idea had evolved from the 
discovery of extraterrestrial radio-frequency radiation by 
Jansky^ in the early 1930's. Figure 2 depicts the radio 
equivalent first employed by Ryle and Vonberg for solar 


^Jansky, K. G., “Electrical Disturbances Apparently of Extraterrestrial 
Origin,*’ Proc Inst. Radio Engineers, N.Y., Vol. 21, p. 1387, 1933. 


observations in the late 1940’s.^ In the same period, a radio 
Lloyd’s mirror interferometer was set up in Australia. Mean- 
while, McCready and colleagues (Ref. 230) used a similar 
technique to study solar radiation and its relation to sunspots. 

Resolution obtained by radio instruments has been steadily 
upgraded by advanced technology in short wavelengths and 
increases in baseline lengths. Many radio sources are well 
resolved by an interferometer with a baseline of about one 
kilometer (km) and a resolution of the order of one arc 
minute. However, by the late 1950’s, it was realized that 
further increases in resolution would be required to study the 
structure of those sources in more detail and to resolve the 
smaller sources. 

V. Modifications to the Basic 
Interferometry 

The basic interferometer system can be modified in many 
ways. One way is to employ several intermediate frequencies, 
each generated by an independent local oscillator. Another 
way is to use complicated local oscillator chains for 
greater frequency versatility and phase-compensating systems. 
Detailed discussions on four modifications, namely, lobe 
rotation, spectral-line interferometry, Very Long Baseline 
Interferometry, and intensity interferometry were given by 
E. B. Fomalont and M. C. H, Wiight."^ We deal here solely with 
the subject of VLBL 

VI. Type of Radio Source 

It is now generally accepted that the radio emission of 
discrete sources such as radio galaxies and quasars is “synchro- 
tron” radiation from electrons moving m a weak cosmic mag- 
netic field at relativistic speeds, i.e., speeds close to the speed 
of light. Also, radio sources may' essentially be classified into 
two different types. One has a large angular extent and is 
strongest at the longer wavelengths. The other is relatively 
compact and is strongest at the shorter wavelengths. There is 
no simple relation between the angular extent of the radio 
emission and the optical emission from galaxies and quasars, 
however compact radio sources are not confined to quasars; 
many are identified with the nuclei of galaxies. 'Furthermore, 
many quasars are large extended radio sources. 

The large radio sources have a complex distribution of radio 
emission that typically extends over several hundred thousand 
li^t years of space, corresponding to angular dimensions of 


^Ryle, M., and D. D. Vonheig, Nature, pp. 158-339, 1946. 

^See Chapter 10 of Galactic and Extragalactic Radio Astronomy, 
edited by G. L. Verchuur and K. X. KeUemiann, 1974. 
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between a few arc seconds to a few arc minutes. On the other 
hand, the compact radio sources are so small and their particle 
densities are so great that, since relativistic particles absorb 
radiation as well as emit it, the source becomes opaque to its 
own radiation at long wavelengths. The self-absorption cutoff 
frequency depends only on the flux density, -the angular size, 
and weakly on the magnetic-field strength. The smaller the 
source, the shorter the wavelength at which it becomes 
opaque. 

Since the compact sources are relatively weak at the long 
wavelengths, where most of the early observations of extra- 
galactic radio sources were made, they remained essentially 
unnoticed for many years until sensitive receivers for short 
wavelengths became available. Further, William A. Dent (Ref. 
95) discovered in 1965 that some compact sources such as the 
quasar 3C273 are variable as observed from their radio 
emissions- 

VII. Development of VLBI 

Longer baselines for conventional interferometers were not 
feasible partly because of the increased cost of the cable 
between the interferometer elements and partly because of 
human and natural obstructions such as roads, rivers, moun- 
tains, and ultimately oceans. To overcome such problems, 
radio relay links were employed (see Fig. 3). The idea is to join 
the radio telescopes by a microwave-relay link. Here the 
local-oscillator signal is transmitted to the mixers to reduce the 
high frequency signal down to an intermediate level for 
returning for correlation. In order for the intermediate 
frequency signals from each antenna to be coherent, they are 
generally derived from a common source. This technique was 
first used by Australian and British radio astronomers to 
obtain baselines of more than 100 km and resolutions better 
than one second of arc. By the early 1960’s the efforts to 
extend the interferometer baseline were carried out success- 
fully at the University of Manchester (Ref. 274). Various 
interferometric techniques have also been applied to polariza- 
tion measurement^ and the mapping of source brightness 
distributions (Refs. 142 and 321). However, in order to resolve 
the compact variable sources with expected dimensions of 
about 0.001 arc second, baselines of size comparable to the 
dimensions of the earth were needed. Moreover, large improve- 
ments in the resolution of radio link interferometers were not 
practical because radio links are limited to line-of-sight 


^Ko, H. C., “Theory of Tensor Aperture Synthesis,” IEEE Trans. 
Antennas and Propagation, Vol. AP-15, pp. 188-190, Jan, 1967, 

Morns, D. V., et al,, “Preliminary Measurements on the Distribution 
of Linear Polarization Over Eight Radio Astrophys. J., Vol. 

139, pp. 560-569, Feb. 1964. 


operation, and installation of large numbers of repeater 
stations was both costly and technically complex. 

For many years radio astronomers had considered the 
possibility of completely eliminating the direct electrical 
connection between interferometer elements by separately 
recording the signals at each end on magnetic tape and later 
cross-correlating the two recorded signals. In order to utilize 
this technique successfully, the following requirements are 
necessary: 

(1) The recordings of the two tapes must be synchronized 
so that the time when a given wavefront is received at 
each station is precisely known. Also, the time tags on 
the recordings should permit reasonably short searches 
for fringes m the correlation process. 

(2) The actual radio signal frequencies are generally too 
high to be recorded directly on magnetic tape. Inde- 
pendent local oscillators must be employed to “hetero- 
dyne*’ the radio signal frequency (typically several 
gigahertz) to a much lower frequency so that it can be 
recorded on magnetic tape. If the two lower frequency 
signals are to be cross-correlated, then the oscillators 
must remain coherent over the observing time (or 
integration time). This implies that the relative phase 
change of the two oscillators must remain small during 
this period so that the change in frequency is less than 
the reciprocal of the recording time. For example, at a 
radio signal frequency of 10 GHz and a 100-second 
integration time, a frequency stability better than one 
part in one trillion (A f/f ~ 10“i^) is required, 

(3) Due to the time delay after heterodyning, phase rota- 
tion is also necessary. 

A typical tape-recording interferometer system and a more 
sophisticated one may be seen in Ref. 164. 

The idea of employing independent-local-oscillator tape 
recording interferometers was considered in the U.S.S.R. as 
early as 1961. However, stable frequency standards and wide- 
band tape recorders needed for high sensitivity were not 
generally available then. The frequency standards generally 
perform two functions (Ref. 340): 

(1) To generate sufficient phase stability that will ensure 
no significant loss of coherence in the two signal 
streams to be cross-correlated. 

(2) To keep time with sufficient accuracy to insure that 
the clock offset error remains constant over the interval 
required for the determination of the instantaneous 
baseline vector. 


148 



A tape-recording interferometer was first used in radio 
astronomy by a group at the University of Florida to study the 
dimensions of the radio storms on the planet Jupiter at a 
frequency of 18 MHz in 1965 (Ref. 40), Since the Jovian 
bursts are so intense, an integration time of considerably less 
than a second and a bandwidth of about one kilohertz (kHz) 
give an adequate sensitivity (requires Af/f = 10““^). Hence, a 
frequency stability of only one in a million and a time 
synchronization of about 1 millisecond from the National 
Bureau of Standards Station WWV were sufficient to maintain 
coherence at the two ends of the interferometer system. With 
this it was possible to determine on the surface of Jupiter the 
dimension of the radio emitting regions less than 1 arc second 
or about 320 km. This resolution is considerably better than 
the highest resolution in photographs made of Jupiter at 
optical wavelengths. 

Tape-recording interferometers for studying the much 
weaker radio emission of radio galaxies and quasars had to 
wait until stable atomic frequency standards and wideband 
tape recorders were commercially available. At that time, two 
systems for tape-recording interferometry were developed 
independently in the United States (Refs. 12 and 248) and in 
Canada (Refs. 26 and 27), 

A. Analog VLBI Development 

The joint Canadian team from the National Research Coun- 
cil and the University of Toronto built an analog recording and 
processing system employing television tape recorders. At each 
site data were recorded at 4 MHz on video recorders (down 
converted from 408 MHz received frequency) on a 90-minute 
reel tape. Then rubidium frequency standards and a special 
analog correlator were used to process the data to simul- 
taneously provide fringe patterns from a number of delays. 
These fringe patterns were then displayed on a multipen 
recorder. 

B, Digital VLBI Development 

A system built by two U.S. groups from the National Radio 
Astronomy Observatory and Cornell University (NRAO- 
Cornell) used a standard computer tape drive to record digital 
data in a 360 kHz band. Rubidium frequency standards were 
also employed here. The incoming signal (down converted 
from 610 MHz received frequency) was first bandlimited, 
infinitely clipped, sampled at 720 kilobits per second (kbps) 
and finally recorded on standard 7-track* computer tape at a 
per-track density of 315 bits per centimeter; hence, a full reel 
of tape held three minutes of data. This system preserves only 
the zero crossing information (i.e., the phase) of the signal. 
The normalized cross-correlation function can still be esti- 
mated from this infinitely clipped data by applying a correc- 


tion factor derived first by Van Vleck in 1943 (Ref. 382), A 
large general purpose computer was used to process the data. 

VIIL Further Developments and 
Experiments on VLBI 

The Canadian system has a greater bandwidth advantage 
but requires more complex special equipment to synchronize 
the tape on playback. In the U S. system digital data were 
simply fed into a large general purpose computer which stored 
and correlated the two data streams. In early experiments 
90 minutes of computer time was typically required to process 
a pair of tapes over 15 delays. Lately, only 50 seconds of 
computer time has been required to search over 7 delays. By 
using special-purpose digital correlators further reductions in 
computer time are possible. 

Both the Canadians and Americans had experienced diffi- 
culties in their early attempts. The Canadians obtained their 
fringes first, both on a baseline in Ontario between Shirley Bay 
and Algonquin Park (183 km) and on a much longer baseline 
between Algonquin Park and Penticton, B.C. (3074 km).'^ The 
NRAO-Cornell group obtained fringes at 610 MHz between 
Riverside, Md. and Green Bank, W. Va. (228 km), followed by 
a successful experiment at 1.665 MHz between Westford, 
Mass., and Green Bank, W. Va: (845 km), with the same digital 
equipment, by a group at the Massachusetts Institute of 
Technology. A hydrogen maser frequency standard was used 
at the Massachusetts site. 

The MIT group (Ref. 408), in cooperation with the North- 
east Radio Observatory Corporation, Haystack Observatory, 
Westford, Massachusetts, undertook a program to improve the 
instruments to allow accurate measurements of group delays 
(the frequency derivative of the phase delay) to make possible 
precise determinations of vector baseline (Ref. 138), radio 
source positions (Ref. 18), polar motion, universal time, and 
earth tides (Ref. 346). 

To measure the group delays with an accuracy greatly 
exceeding the inverse of the recorded bandwidth, the method 
of bandwidth synthesis was devised (Refs. 310 and 408) where 
the receiver passband is switched to sample signals over a 
wideband. This bandwidth synthesis technique is also 
employed in the VLBI development at JPL (Refs. 273, 
368-376). 

First long-baseline measurements with tape-recording were 
made in 1967 on the baseline across Canada and the U.S. at 


^Their experimental results were presented at the Ottawa meeting of 
the International Scientific Radio Union on May 23, 1967 (Refs. 26 
and 27). 
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radio wavelengths of 75, 50, and 18 centimeters* The results 
confirmed many small size galaxies and quasars. But many 
were still unresolved and higher resolution was needed. 
This was rapidly extended in a series of cooperative inter- 
continental experiments conducted in 1968 and 1969 by 
American, Swedish, and Australian radio astronomers. With a 
6-centimeter wavelength, a 0.001 -arc second resolution was 
obtained on the longest baselines. In spite of this, some small 
sources still remain unresolved. Hie 10^ km Calif ornia- 
Australia baseline was already more than 80% of the earth’s 
diameter, so further significant increases in the physical base- 
line were not feasible without expensive procedures of setting 
up stations in space or on the moon. 

A cheaper and simpler alternative was to observe at shorter 
wavelengths. Outside North America, however, only two radio 
telescopes were suitable for operation at short wavelengths and 
large enough to yield sufficient sensitivity for long-baseline 
interferometry, both of which were in the U.S.S.R. Due to 
distance, government regulations, lack of exchange in tech- 
nology, and other problems, arranging such joint experiments 
was both difficult and formidable. 

The first U.S.-U.S.S.R. experiment was completed late in 
1969, using the baseline between the 43 -meter radio telescope 
at NRAO and the 22-meter radio telescope on the shore of the 
Black Sea in the Crimea (Ref. 388). Later in Spring of 1971, a 
second experiment was conducted, involving, in addition to 
the above, the 64-meter radio telescope at- Golds tone, ^ Califor- 
nia and the “Haystack” telescope in northern Massachusetts.^ 
More than 20 investigators from eight institutions in both 
countries participated and it incl uded o bservations of inter- 
stellar clouds of water vapor as well as of radio galaxies and 
quasars. 

The data from the Goldstone-Crimea basehne obtained at a 
wavelength. of 3.5 cm (X-band) provided the highest resolution 
(about 0.0003 arc sec) achieved so far in the study of radio 
galaxies and quasars. The measurements of the water vapor 
clouds which were made at shorter wavelengths gave even 
hi^er resolution. 

From the various data taken, it was evident that some 
quasars appear to expand faster than the speed of light. 
Martin J. Rees (Ref. 297) of the University of Cambridge has 
suggested the “super light velocity” theory for interpreting 
such a phenomenon. If the radio source expands at a velocity 
close to the speed of light, then since it takes a finite time for 
radio signals to reach an observer, the signal arriving from the 


^Managed by JPL for NASA. 
^Managed by MIT. 


receding part of the source will have originated at an earlier 
time, when it was closer to the point of origin than the signals 
from the approaching parts. Under these conditions the 
apparent velocity of expansion may in fact exceed the speed 
of light. 

Evidence that this might be important in the quasar 3C279 
was first obtained in a series of transpacific observations made 
between 1968 and 1970 by a joint team from Australia (at 
Tidbinbilla) and the California Institute of Technology (at 
Goldstone). This group found that a component of the source 
that had first appeared in 1966 had reached a diameter of 
about 0.001 arc second by the end of 1969. This corresponds 
to a linear diameter of about 12 light-years. Thus, the apparent 
expansion velocity was indeed about twice the speed of light, 
as theorized by Rees. 

More detailed measurements were made of 3C279 in Octo- 
ber, 1970, by a group from MIT, NASA, and the University of 
Maryland employing a transcontinental baseline. These obser- 
vations were originally designed to measure the gravitational 
bending of the radio signals from 3C279 as it approached the 
sun on October 8, hut they showed clearly that the compact 
source in 3C279 was complex and appeared to have at least 
two components separated by 0.00155 arc second, or about 
20 li^t-years. This source was observed again in February 
1971 by the same group and also by investigators from the 
NRAO, Caltech, and Cornell group, using the same baseline 
and techniques. The source appeared to double in only four 
months, but the separation was greater by two light-years. 
Hence the two components seem to be receding from a com- 
mon point of origin with an apparent velocity about three 
times the speed of light. 

Although the sensitive long-baseline interferometer systems 
were initially developed to study the compact extragalactic 
radio source, the technique has also been used to study the 
radio emission from interstellar clouds of hydroxyl radicals 
(OH) and water vapor. Typically the clouds are dispersed over 
a volume of several li^t-years across, although the individual 
components are as small as one astronomical unit (the 
average distance between the sun and the earth). 

One limitation of tape-recording interferometry has been 
the lack of sufficient oscillator stability to determine the phase 
of the interference pattern; the improvement is to be expected 
from the use of hydrogen masers as frequency standards. Even 
with infinitely stable oscillators, there are still problems 
(Ref. 339) resulting from fluctuations in the path length 
through the atmosphere and ionosphere at the two widely 
separated observing sites. 
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Another limitation is the antenna flexure, but this type of 
limitation is often small and can be calibrated with a residual 
error of less than 1 mm. 


IX. Limitations on Baseline Length 

Extension of the baseline cannot always solve the resolu- 
tion problem. The present evidence indicates that only for 
wavelengths less than about 10 cm can baselines much 
greater than the diameter of the earth be effectively used. 
Hence, hydroxyl emission regions, which radiate at 18 cm, 
and pulsars, which radiate most strongly at meter wavelengths, 
are not likely to be targets for a space interferometer. 

There is an even more fundamental limit to the maximum 
baseline when dealing with the synchrotron sources. If the 
dimensions of the synchrotron system are below a critical size, 
the relativistic particles quickly lose all their energy by inverse 
compton scattering and so do not have time to radiate radio 
energy. This critical angular size is proportional to the wave- 
length at which the flux density is the highest. Since the reso- 
lution of a fixed length interferometer is inversely proportional 
to the wavelength of observation, the two effects cancel each 
other if the observations are made near the wavelength of 
maximum intensity (as they usually are to achieve the hipest 
sensitivity). Thus, the maximum baseline needed to resolve 
synchrotron sources is nearly independent of the wavelength 
of observation; in fact, it is comparable to the diameter of the 
earth for the stronger radio galaxies and quasars. Baselines in 
space or on the moon will therefore probably not be necessary 
to study even the smallest radio nuclei and quasars. 

Besides extending the baseline to achieve high resolution, 
intermediate baselines may be made very effective by linking 
together to form an ultra-high-resolution array (Refs. 74 
and 171); but this can be difficult. Moreover, the task of 
collecting all the tape recordings at a common location* and 
correlating all possible pairs of telescopes is formidable. An 
attractive alternative is to telemeter the data from each anten- 
na to a common site by means of , synchronously orbiting 
satellites, dispensing with the tape recording system complete- 
ly. This technique was successfully demonstrated via the 
Harmes satellite as reported in Ref. 422. 


X. Radio Interferometry for Earth 
Physics 

,The techniques of VLBI have also been applied to geophysi- 
cal studies at JPL. The study of various faults, earthquakes, 
and other, geophysical applications often requires a wide range 
of baseline lengths. A project that concentrates on this aspect 
of VLBI is refened to as the Astronomical Radio^ Interfero- 
metric Earth Survey (ARIES) Project (Ref. 218). The essence 
of the project is to use a pair of antennas, one fixed and the 
other portable, to provide the desired wide range of baseline 
lengths. The system characteristics, error sources, and other 
system details and requirements are given in Refs. 218, 272 
and 273. 

A mobile VLBI geodetic system is currently under 
development at JPL for NASA application to earth crustal 
dynamics studies, as well as for technology transfer to the 
National Geodetic Survey (NGS), the U.S. Geological Survey 
(USGS), and other geodetic user communities. This system 
consists of a 4-meter high mobility station and a 13 -meter 
transportable base station. The system characteristics are given 
in Table 1, The various error sources (both random and sys- 
tematic) for the mobile VLBI geodetic system are tabulated in 
Table.2. Certain applications of such a geodetic system are: 

(1) Environmental monitoring related to natural resource 
extraction such as subsidence resulted from geothermal 
field heat removal, oil drillings, and ground water 
pumping. 

(2) Global tectonics and regional crustal deformations for 
uses in the development of a fundamental scientific 
understanding of geodynamic processes and earthquake 
precursory model development in the pursuit of a reli- 

• able productive capability. 

Other possible VLBI applications to earth physics are global 
geodesy, tidal oscillations, crustal-block motions (including 
continental drift), polar motion, earth rotation, precession and 
nutation (including a test of general' relativity), obliquity of 
ecliptic, shape of sea surface geopotential and global time and 
frequency synchronization. More detailed discussions on these 
applications are given in Ref. 339. 


^^Renzetti, N A., “Radio Interferometric Geodetic System Applica- 
tions,” A JPL internal interoffice memo, dated June 1, 1978. 
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Table 2. Errors m mobile VLBI geodetic system 


Table 1. Mobile VLBI geodetic system characteristics 


Characteristic 

Mobile 

Transportable 

base 

Antenna diameters 

4 m 

13 m 

Antenna efficiency 

50% 

50% 

Receiver temperature (X-band) 

70 Kelvin 

30 Kelvin 

Band^vidth synthesized 

400 MHz 

400 MHz 

Radio source strength 

2 Jansky 

2 Jansky 

Digital data recording rate^ 

112 Mb/s 

112 Mb/s 

Observation time per source 

720 sec 

720 sec 

Note: Signal-to-noise ratio 

25*1 


Average delay 
precision 

22 psec (7 mm) 

Baseline measure- 
ment precision 
limit 

5 mm given 6 hr of data 
(single work shift) 

^3-channel bandwidth synthesis 


Errors 

Distance 

equivalent, 

mm 

Random Errors 


Delay precision 

7 

Time and frequency (H-masers, Af/f = 5 X 10“^^) 
Troposphere 

10 

Dry (surface meteorology) 

10 

Wet (water vapor radiometers) 

10 

Ionosphere (8.4 GHz, X-band) 

5 

Radio source catalog (0.003 arc sec) 

7 

Root sum square 

21 

Baseline vector precision given 12 hours of data 
Systematic Error Sources (Baseline <500 km) 

10 

UTl/polar motion (10-cm accuracy, supplied by 
external sources) 

8 

Radio source catalog (0.003 arc sec accuracy to be 
held essentially constant for monitoring missions) 

Troposphere 

4 

Dry (surface meteorology) 

10 

Wet (water vapor radiometry) 

7 

Ionosphere (X-band) 

5 

Antenna positioning relative to geodetic marks 

5 

Root sum square 

17 

Combined random and systematic 
error sources accuracy) 

20 
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Pig. 1. Michelson— Pease Optical Interferometer. (Light from a distant radio source is reflected 
from the outer mirrors to the inner mirrors and then optically combined at the projection screen. 
Moving one of the mirrors makes one light path longer than the other. The two out-of-phase 
beams, when combined at the screen, will interfere with each other to provide ‘‘fringes.” A 
photograph of actual double-slit interference fringes is shown on the right.) 


MIXER 



MIXER 


MIXER 




MIXER 


I. F. 


Fig. 2. Basic Radio Interferometer. (Two radio telescopes, sepa- 
rated by a few kilometers, are linked by cables. The high radio 
frequency (R.F.) signal from the radio source is mixed with a local 
oscillator signal to provide an intermediate frequency (LF.) signal 
which is suitable for further data processing.) 


Fig. 3. General Radio-link Interferometer Concept. (The local- 
oscillator signal is transmitted to the mixers and the I.F. signals are 
returned via microwave-relay iinks. The local-osciilator signal at 
each end is derived from a common source so as to provide 
coherent (in phase) I.F. signals from each antenna.) 
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An Extensive Bibliography on Long 
Baseiine Interferometry 

B. Benjauthrit 

Office of Tracking and Data Acquisition 


This article presents an extensive bibliography on the subject of long baseline inter- 
ferometry, starting from the time of Albert A. Michelson (1890) up to the present time. 
It contains over 400 references, including areas of long baseline interferometry 
applications. 


Over the past few decades, the subject of long baseline 
interferometry has grown from a simple concept of an optical 
device for detecting light fringes to a highly theoretical and 
sophisticated system assisted by computers and modern elec- 
tronic components. Interest in this subject is still increasing. 

areas of applications are rapidly widening. These include 
radio astronomy, -radio science, geodesy^, tectonics^, space 
navigation, and seismology. 


^ A blanch of applied mathematics that determines the exact positions 
of points and the figures and areas of large portions of the earth’s 
surface, the shape and size of the earth, and the variations of ter- 
restrial gravity and magnetism {Webster*s Collegiate Dictionary, 
1975 ). 

branch of geology concerned with structure, especially with folding 
and faulting (ibid.) 


With the above wide applications, it is not surprising that 
the number of references on the subject is almost inexhaust- 
ible. As one searches througli the literature, one often wonders 
which references one should select. The references selected 
here were considered most relevant to the subject of long 
baseline interferometry. However, certain papers and texts 
which have been found helpful in understanding the subject 
have also been included. 

Although there are over four hundred references listed in 
this bibliography, it is by no means complete. It may, how- 
ever, be considered one of the most extensive bibliographies 
on the subject. The bibliography is based on a computer 
data base search performed by the JPL Library. It is currently 
stored on the DACONICS computer system at the Laboratory. 
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In addition to searching through the appropriate references 
on the subject in the literature, the NASA/RECON computer 
printouts (the primary and alternate data banks and engi- 
neering index from 1962 to April 20, 1978), the INSPEC 
data base, and a Report Bibliography of the Defense Docu- 
ment Center were utilized. Several* references were obtained’ 
through private communications with researchers in the field. 

One can generally obtain a bibliography on the subject 
from computer data banks such as those mentioned above, 
with abstracts and other additional information. However, 
due to the volume of the data, one is often discouraged from 
looking through the printouts. With this present bibliography, 
one should find it much more expeditious and convenient to 
identify appropriate material. Tlie references are listed alpha- 


betically by the principal authors’ names. Papers by the same 
author are listed chronologically. If an article is published 
in more than one journal (with the same title and authors), 
the alternative sources are also noted at the end of the refer- 
ence. Moreover, those papers which can be ob tamed from the 
National Technical Information Service (NTIS) are indicated 
by the words, “Avail. NTIS.” Also, the Defense Document 
Center’s identification numbers should be found useful. 

According to the computer data banks mentioned above, 
over 1800 references have been published on VLSI and its 
related topics. Clearly, it was impossible to include all the 
available references in this bibliography. Many excellent 
papers on the subject have had to be omitted, as have the 
majority of the papers written on applications of VLBL 
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